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Abstract 

It  is  the  purpose  of  this  research  to  examine  the  engineering  process  required  to 
adequately  predict  the  total  resistance  of  unique  hullforms.  A  comparison  of  the  significant 
physical  effects  present  in  monohull  and  advanced  hullforms  testing  is  made.  Hie  historically 
accepted  method  of  scaling  using  Froude's  Hypothesis  is  studied  and  its  underlying 
assumptions  are  considered  in  detail.  .An  overall  design  philosophy  is  created  for  achieving 
resistance  minimization  through  a  range  of  operating  speeds.  This  discussion  is  applied  to  a 
specific  hullform  called  SLICE  which  is  an  adaptation  of  the  Small  Waterplane  Area  Twin 
Hull  (SWATH)  concept. 

To  further  develop  the  body  of  available  SLICE  resistance  data,  three  models  were 
built  and  tested  at  the  US  Naval  Academy's  380  feet  towing  tank  in  Annapolis  Maryland.  The 
results  of  these  tests  is  reported  on  and  a  video  demonstration  of  their  performance  is 
provided  as  part  of  this  report. 

Finally,  a  numerical  method  approach  to  the  solution  of  Mitchell's  integral  is  presented 
as  a  possible  means  of  calculating  SLICE  wave  making  resistance.  This  is  considered  as  an 
intermediate  step  tool  in  design  optimization,  useful  after  preliminary  model  tests  are 
completed  and  before  final  designs  are  decided  for  powered  model  testing. 
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Chapter  1  Introduction 


The  majority  of  the  world's  commerce  is  moved  by  monohull  displacement 
ships.  Ironically,  in  this  age  of  supersonic  air  travel,  guaranteed  overnight  mail 
delivery  and  fibre  optic  data  transfer  at  the  speed  of  light,  modern  displacement  ships 
are  capable  of  speeds  no  faster  than  those  built  fifty  years  ago.  The  apparent  lack  of 
progress  in  high  speed  ocean  travel  stems  from  the  physical  realities  of  the  interaction 
between  the  displacement  monohull  form  and  the  water  in  which  it  travels. 

The  total  resistance  curve  for  a  ship  is  generally  plotted  as  the  Total  Drag 
Coefficient,  CT  versus  Froude  Number.  These  are  expressed  as: 


^V/STgLT  0-]) 

*r-(ip  V2S)C,  (1-2) 

where  L  is  some  characteristic  length,  usually  length  of  the  waterline  of  the  body,  p  is 
the  density  of  the  water,  V  is  the  ship's  velocity  and  S  is  the  Wetted  Surface  Area. 
Total  resistance  is  generally  recognized  as  the  additive  effects  of  friction,  wave 
making  and  spray  generation  of  both  the  main  body  and  its  appendages.  The  extent  to 
which  each  of  these  components  effects  the  total  resistance  is,  however,  the  subject  of 
ongoing  debate  even  though  it  has  been  over  one  hundred  years  since  Froude  gave  us 
this  theory.  For  a  monohull  displacement  ship,  the  wave  making  resistance  portion  of 
the  total  resistance  curve  has  the  characteristic  shape  shown  in  Figure  (1-1).  Since 
resistance  is  obtained  by  multiplying  the  coefficient  shown  in  Figure  (1-1)  by  V2,  the 
total  resistance  curve  becomes  nearly  vertical  and  practically  insurmountable  by 
conventional  hull  forms  at  the  speeds  shown. 
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Figure  1-1:  Characteristic  Wave  Making  Resistance  Curve 
(from  LMSC  presentation) 


For  Naval  Architects,  the  search  is  ongoing  in  advanced  hull  form  design  for  a 
practical  method  to  "get  over  the  hump".  Most  advanced  hull  forms  do  not  resemble 
the  classic  displacement  monohull  shape  on  which  the  vast  majorin'  of  Naval 
.Architecture  experience  is  based.  The  purpose  of  this  research  is  to  examine  the 
underlying  assumptions  regarding  commonly  accepted  full  scale  resistance  prediction 
methods  in  an  effort  to  understand  the  important  features  of  design  optimization  and 
draw  some  conclusions  regarding  resistance  minimization  methods  .  The  Lockheed 
Missile  and  Space  Company  in  conjunction  with  Pacific  Marine  (LMSC/PM)  has 
proposed  a  hull  form  which  is  a  modification  of  Small  Waterplane  .Area  Twin  Hull 
(SWATH)  concepts.  Their  proposal,  called  SLICE,  will  be  presented  to  the  United 
States  Government  as  a  177  Lton  vessel  in  an  Advanced  Technology  Demonstration 
(ATD)  in  June  1996.  A  sketch  of  this  hullform  is  shown  in  Figure  1-2.  The  SLICE 
vessel  will  be  the  subject  hullform  for  the  research  conducted  as  part  of  this  paper. 
SLICE  is  protected  by  Lockheed  under  patent  law. 


Figure  1-2:  Artist's  Concept  of  SLICE  from  an  unpublished  LMSC, 
Pacific  Marine,  ONR  Article. 


1.1  The  SLICE  Concept 

High  speed  is  not  an  inherent  quality  of  SWATH  ships.  Although  there  may 
be  some  reduction  of  the  wave  making  resistance  due  to  smaller  waterplane  area,  there 
is  also  an  increase  in  the  total  wetted  surface  which  increases  frictional  resistance.  In 
SWATH  ships,  high  speed  is  a  design  feature.  When  compared  to  an  equivalent 
displacement  monohull.  SWATH  design  must  be  considered  as  a  trade-off  between  a 
definite  increase  in  frictional  resistance  and  the  potential  for  a  decrease  in  wave 
making  resistance.  TTie  design  question  becomes  how  to  ensure  that  the  decrease  in 
wave  making  resistance  is  greater  than  the  gain  in  frictional  resistance. 

LMSC's  SLICE  concept  differs  from  traditional  SWATH  design  in  that  the 
demi-hull  is  replaced  by  two  shorter  pods.  Each  pod  is  attached  to  the  main  box 
structure  by  a  single  strut.  It  is  through  this  design  feature  that  LMSC  intends  to 
minimize  the  impact  of  wave  making  resistance.  By  designing  the  vessel  with  an 
unusually  short  characteristic  length,  the  hump  of  the  wave  making  curve  will  in 


theory  be  driven  to  a  very  low  speed  regime.  Since  resistance  is  a  function  of  V2,  any 
left  shift  of  the  curve  could  have  a  substantial  minimizing  effect  on  the  total 
resistance  associated  with  the  hump.  It  is  further  theorized  that  in  the  high  speed 
regime,  total  resistance  will  be  dominated  by  frictional  effects  and  therefore  design 
efforts  focused  on  minimizing  the  wetted  surface  of  the  hullform  will  yield  the  best 
result.  LMSC  has  chosen  the  length  of  the  submerged  pod  as  the  characteristic  length 
critical  to  design  considerations.  It  should  be  emphasized  that  this  hullform  does  not 
eliminate  the  hump  effect  characteristic  of  the  wave  making  resistance  curve.  It  does 
not  specifically  cause  wave  cancellation  effects  nor  does  it  attempt  to  reduce  viscous 
drag  effects.  Rather,  the  apparent  thrust  of  the  LMSC  theory  is  to  stagger  the  peak 
effects  of  wave  and  skin  resistance,  thereby  reducing  the  overall  resistance  value  at 
specific  design  speeds. 

1.2  LMSC  Model  Test  Synopsis 

LMSC  has  conducted  two  different  sets  of  model  tests  at  the  US  Navy’s 
Surface  Warfare  Center,  Carderock  Division,  Bethesda,  Maryland  (NSWC).  In  the 
first  series  of  tests,  called  the  HM&E  tests,  several  different  configurations  of  struts, 
pods  and  foils  were  used  to  make  a  determination  of  the  overall  best  configuration  to 
use  as  the  baseline  for  the  production  vessel.  These  tests  were  conducted  in  the  Spring 
1994. 

The  second  set  of  tests,  also  conducted  at  NSWC,  were  used  to  obtain  detailed 
resistance  and  control  surface  performance  data  for  the  baseline  production  vessel. 
These  tests,  conducted  in  August  1994  and  continued  in  December  1994  were 
observed  by  several  MIT  students.  The  model  that  was  tested  was  a  l/8th  scale  of  the 
production  vessel.  It's  construction  was  fiberglass  bodies  mounted  to  an  aluminum 
box  frame.  For  all  the  tests  conducted,  the  model  was  fixed  in  pitch  and  heave  and 
data  was  collected  for  forces  and  moments  about  the  three  principle  axes.  There  were 
two  major  groupings  of  data  collection  runs  performed.  The  first  group  was  called  the 
dynamic  control  tests  and  was  primarily  used  to  judge  the  required  effectiveness  of 
control  surfaces  to  hold  the  ship  at  a  constant  draft  with  an  even  trim.  Each  run  had  a 
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slightly  different  model  configuration  with  canards,  stabilizers  and  rudders  tested  in 
different  locations  and  with  varying  angles  of  attack.  Each  configuration  was  tested  a 
single  time  unless  there  was  some  specific  problem  noted  in  the  conduct  of  the  test. 
The  second  major  grouping  of  tests  consisted  of  rather  major  modifications  to  the 
model  configuration  in  an  attempt  to  improve  the  design.  Haunches  on  top  of  the  struts 
were  removed  to  allow  for  deeper  strut  length,  wedges  located  at  the  top  of  the  struts 
were  added  and  removed  to  experiment  with  better  surface  piercing  shapes,  the  overall 
draft  of  the  model  w  as  adjusted  to  force  a  deeper  submergence  of  the  rather  bluff  pods 
and  finally,  a  new  model  was  tested  in  which  the  longitudinal  position  of  the  struts 
relative  to  the  pods  was  altered  to  modify  the  flow'  over  the  pods. 

Within  these  two  major  test  categories,  the  test  procedure  remained  consistent. 
The  model  remained  fixed  in  pitch  and  heave  for  all  tests.  Two  different  speed 
regimes  were  tested,  a  lower  speed  run  consisted  of  twelve  speed  increments  ranging 
from  3.2-  8.8  feet  per  second  (fps)  corresponding  to  full  scale  speeds  of  9-25  knots.  A 
high  speed  run  used  seven  increments  from  5.7-12.4  fps  which  corresponds  to  full 
scale  speeds  of  10-35  knots.  All  speed  increments  were  fit  into  a  single  pass  along  a 
tow  tank  length  of  approximately  1200  feet.  All  of  the  dynamic  control  tests  were 
conducted  only  in  the  lower  speed  regime,  the  major  model  modification  tests  were 
conducted  in  both  speed  regimes. 

1.3  Research  Motivation 

The  US  Navy  has  potential  applications  for  relatively  small  (200-700  Lton) 
high  endurance  craft  that  embody  the  seakeeping  attributes  of  SWATH.  The  Navy's 
interest  in  the  SLICE  program  is  demonstrated  by  the  funding  the  Navy'  has  committed 
to  the  ATD  that  LMSC  will  perform  in  June  1996.  It  was  LMSC/PM’s  entrepreneurial 
spirit  that  initially  brought  the  SLICE  concept  to  the  attention  of  the  Navy'  but  this 
same  entrepreneurial  spirit  has  the  potential  to  limit  the  usefulness  of  the  SLICE 
design  as  well.  LMSC/PM’s  allocated  timetable  from  concept  development  to  final 
product  delivery  is  just  over  two  years,  very  short  by  Navy  shipbuilding  standards.  As 
a  result  LMSC  must  focus  on  a  relatively  narrow  set  of  engineering  problems  that 
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specifically  address  the  needs  of  the  177  Lton  ATD  vessel.  Broader  based  questions, 
such  at  an  in-depth  understanding  of  the  effects  of  body  interactions,  how  these 
interactions  vary  with  modifications  to  the  basic  dimensions,  the  practicality  of  scaling 
the  SLICE  concept  into  a  different  size  regime  (1500  Lton,  5000  Lton,  etc.)  and  a 
relative  comparison  of  the  merits  of  a  SLICE  ship  as  compared  to  a  SWATH  or 
Monohull  design  for  a  given  mission  are  all  facets  of  the  development  of  this  concept 
that  are  not  specifically  addressed  by  LMSC  in  the  ATD. 

This  thesis  will  attempt  to  add  to  the  developing  body  of  SLICE  knowledge  by 
focusing  on  the  resistance  and  powering  aspects  of  the  SLICE  vessel  in  an 
academically  rigorous  manner.  In  so  doing,  three  benefits  will  be  realized: 

1)  Independent  research  into  this  hull  form  will  give  both  the  US  Navy  and 
LMSC  a  second  hydrodynamic  viewpoint  of  the  177  Lton  vessel  prior  to  the  ATD. 

2)  Development  of  a  sufficient  body  of  data  will  aid  optimization  of  the  177 
Lton  hull  form. 

3)  A  broader  understanding  of  the  hydrodynamic  features  of  the  hull  form 
will  facilitate  deciding  the  merits  of  pursuing  designs  in  other  displacements  more 
practical  to  the  US  Navy. 

1.4  Approach 

To  achieve  the  objectives  outlined  above,  there  are  several  questions  raised  by 
the  LMSC  efforts  that  require  answers.  These  questions  and  the  method  that  will  be 
pursued  for  answering  them  is  outlined  in  this  section. 

THEORY 

LMSC  theorizes  that  the  appropriate  characteristic  length  for  SLICE  is  the 
length  of  a  single  pod.  In  conventional  designs.  cFUs  based  on  the  length  of  the 
waterline  (LWL).  Representative  characteristic  lengths  for  SLICE  would  seem  to  be 
bounded  on  the  low  end  by  the  length  of  a  single  strut  and  on  the  high  end  by  the 
length  between  the  leading  edge  of  the  forward  strut  to  the  trailing  edge  of  the  aft 
stmt.  LMSC's  choice  of  an  underwater  body  for  the  characteristic  length  is  not  an 
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obvious  one.  Understanding  the  appropriate  characteristic  length  for  the  SLICE 
configuration  is  necessary  to  understand  the  magnitude  of  savings  that  can  be  expected 
over  competing  hull  forms. 

The  interaction  between  bodies  is  the  second  area  for  further  study.  There  are 
two  specific  situations  that  require  detailed  analysis.  The  first  is  the  effect  a  strut  and 
its  attached  pod  have  on  one  another,  the  second  is  the  effect  that  a  strut/pod 
combination  have  on  any  of  the  other  three  strut/pod  combinations  that  form  the 
underwater  body.  Understanding  these  effects  is  necessary  for  design  optimization  and 
also  deciding  the  practical  upper  limit  of  displacement  for  a  SLICE  ship. 

Finally,  the  LMSC  proposed  design  places  the  forward  struts  inboard  of  the  aft 
struts  as  illustrated  in  Figure  2.  This  arrangement  may  cause  some  significant 
outward  lateral  forces  on  the  struts  that  must  be  considered  in  the  structural  design. 

TESTING 

There  were  several  procedural  questions  raised  by  the  model  tests  conducted  by 
LMSC  in  August.  Keeping  the  model  fixed  in  pitch  and  heave  was  the  most 
significant  issue.  Although  it  was  intentionally  done  to  facilitate  the  dynamic  control 
analysis,  it  also  left  open  significant  questions  regarding  how  the  hull  form  behaves  at 
different  speeds.  Considering  the  recognized  sensitivity  that  SWATH  type  hulls  have 
in  tons  per  inch  immersion  (TPI)  and  the  moment  to  trim  one  inch  (MT1)  this  may  be 
significant. 

The  short  duration  of  each  test  run  speed  increment,  the  lack  of  repetition  of 
tests,  and  the  truncation  of  most  data  collection  runs  at  the  full  scale  equivalent  of  25 
knots  for  a  ship  required  to  achieve  30  knots  raises  questions  regarding  the 
thoroughness  of  the  results.  Although  this  is  understandable  considering  the  expense  of 
model  testing  and  the  limited  time  available,  more  testing  must  be  completed  before 
the  SLICE  body  of  data  can  be  utilized  to  make  design  decisions. 

OPTIMIZING 

There  are  several  striking  features  incorporated  in  the  LMSC  design.  The  most 


striking  is  the  bluffness  of  the  leading  sections  of  the  pods  and  struts,  especially  for  a 
high  speed  craft.  It  would  appear  that  minimizing  the  wetted  surface  to  displacement 
ratio  was  LMSC's  driving  concern  in  hullform  selection.  Although  this  approach 
minimizes  frictional  resistance,  it  also  introduces  significant  form  drag.  With  a  pod 
length  to  diameter  ratio  (L/D)  of  4.2,  SLICE  cannot  be  thought  of  as  a  slender  body. 
This  introduces  another  concern  in  evaluating  SLICE  performance.  The  vast  majority 
of  traditional  resistance  analysis  has  slender  body  theory  incorporated  to  some  extent. 
Any  evaluation  process  undertaken  must  be  carefully  scrutinized  to  ensure  that  the 
SLICE  hullform  does  not  violate  this  imbedded  assumption.  A  design  approach  which 
more  thoroughly  addresses  the  entire  speed  spectrum  is  preferable  even  with  a  specific 
mission  speed  identified.  All  speeds  from  zero  through  the  mission  speed  must  be 
attainable  so  some  design  consideration  must  be  given  to  the  entire  speed  spectrum. 
With  eight  component  bodies  (  4  struts  and  4  pods),  wave  cancellation  effects  through 
body  interaction  should  be  a  design  focus. 

RESEARCH  METHOD 

In  an  attempt  to  answer  these  questions  the  following  work  was  accomplished: 

1 )  A  literature  review  that  encompassed  one  hundred  years  of  ship  resistance 
research  was  performed.  This  was  accomplished  to  ensure  that  the  full  body  of 
assumptions  made  in  this  field  of  engineering  is  understood. 

2)  A  study  of  model  testing  techniques  was  made  to  develop  an  appreciation 
for  proper  test  methods. 

3)  Series  58  resistance  data  was  studied  for  potential  pod  optimization. 

4)  A  literature  survey  was  conducted  to  ascertain  the  state  of  the  science  in 
spray  drag  predictions. 

5)  Three  models  were  built  and  tested  to  expand  the  overall  data  bank. 

6)  A  numerical  calculation  method  for  wave  making  resistance  based  on  the 
solution  to  Mitchell’s  integral  was  attempted  in  an  effort  to  expedite  the  hullform 
optimization  process. 
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Chapter  2  Theory 


The  difficult)  in  predicting  full  scale  resistance  is  in  understanding  the  extent 
to  which  different  physical  phenomena  contribute  to  total  resistance.  Although  effects 
such  as  boundary  layer  development,  wave  making,  spray  generation  and  flow 
separation  are  known  to  be  present  as  a  body  moves  through  the  water,  methods  for 
quantitatively  evaluating  the  contribution  of  each  of  these  effects  to  total  resistance 
and  the  extent  to  which  these  phenomenon  may  interact  continues  to  be  the  subject  of 
ongoing  debate. 

This  chapter  studies  the  two  available  approaches  to  resistance  prediction, 
model  testing  and  analysis.  The  objective  is  to  identify  those  methods  that  will  be  best 
suited  for  SLICE  resistance  predictions.  Since  the  vast  majority  of  available  resistance 
data  was  developed  from  displacement  monohull  analysis,  care  must  be  taken  to 
ensure  that  underlying  assumptions  to  these  methods  are  not  violated. 

2.1  Model  Testing 

The  classical  approach  to  model  testing  is  well  documented.  Principles  of 
Naval  Architecture1  fPNA),  and  Rawson  and  Tupper2  are  probably  the  two  most 
recognized  authorities  outlining  this  process.  The  textbook  approach  is  definitively 
stated  and  easy  to  follow.  Nonetheless,  the  SLICE  experience  has  shown  that  specific 
model  testing  methodology  and  subsequent  data  reduction  can  be  hotly  contested.  The 
difficulty  lies  in  deciding  if  the  model  tested  accurately  replicates  the  physical 
phenomena  that  will  be  present  in  the  full  scale  ship  and  what  to  do  with  the  model 
scale  data  once  it  is  obtained. 

As  background,  the  model  testing/scaling  process  outlined  by  William  Froude 
is  presented.  In  this  discussion  the  subscript  m  refers  to  the  model  and  the  subscript  s 


1  Comstock.  J  P  .  Principles  of  Naval  Architecture.  SNAME,  New  York,  1967.  Chapter  7. 

2  Rawson.K.J.  &  Tupper.E.C.  Basic  Ship  Theory.  3rd  Edition,  Vol  2,  New  York.1983. 
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refers  to  the  full  scale  ship: 

1)  Build  a  model  that  is  geometrically  similar  to  the  full  scale  ship.  This  is 
accomplished  by  choosing  a  scale  factor  to  apply  to  all  linear  dimensions  of  the  ship 
The  scale  factor  X  is  then  defined  as: 


Lu. 

X=_jb?P_  (2-1 ) 

L 

^ model 

2)  Test  the  model  at  equivalent  speeds.  Through  dimensional  analysis,  it  is 
easily  shown  that  the  full  scale  and  model  scale  speeds  are  related  through  equation 
(1-1).  Specifically: 


Khip  ^ model  * 


(2-2) 


3)  Convert  the  total  model  resistance  measured  Rtm  to  a  coefficient  of  total 
resistance  C,m  through  relationship  (1-2). 

4)  Calculate  the  coefficient  of  frictional  resistance  Ctnr  The  method  for 
accomplishing  this  will  be  discussed  shortly. 

5)  A  quantity  called  the  coefficient  of  residuary  resistance  C^  is  obtained  by 
subtracting  Ctm  from  Ctm.  The  reason  Ctm  is  not  subjected  to  the  scaling  law  expressed 
in  (2-2)  will  be  discussed  shortly. 

6)  The  full  scale  coefficient  of  residuary  resistance  Crs  is  then  expressed  at  the 
same  Froude  Number  as  Cm.  Specifically, 

C  =C  (2-3) 

7)  The  full  scale  coefficient  of  frictional  resistance  Cfs  is  calculated  in  the 
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same  manner  as  in  step  4  and  is  added  to  Crs  to  obtain  the  full  scale  coefficient  of 
total  resistance  Cls. 

8)  Full  scale  total  resistance  Rls  is  then  found  using  (1-2). 

This  process  can  be  thought  of  as  a  system  of  one  equation  and  one  unknown. 
The  unknown  is  Cre;  the  equation  is  Q,  =  C,,  +  Cts.  This  method  of  scaling  data  is 
useful  only  if  the  following  conditions  exist: 

1)  A  geometrically  similar  model  can  be  built  and  tested  in  compliance 
with  the  equalities  set  forth  in  (2-1)  and  (2-2). 

2)  The  model  is  tested  in  a  condition  that  can  be  considered  equivalent 
to  the  normal  operating  condition  of  the  full  scale  ship. 

3)  Frictional  resistance  can  be  reasonably  calculated  at  both  the  model 
and  full  scale  size. 

4)  All  other  physical  phenomenon  present  must  either  obey  equality 
(2-3)  and  therefore  be  lumped  into  the  residuary  component  or  be  of  such  little 
concern  that  ignoring  them  is  inconsequential  to  the  final  outcome. 

If  the  last  condition  can  not  be  satisfied  for  a  particular  phenomenon,  then  its 
effect  must  be  calculated  at  both  the  model  and  full  scale  size  and  algebraically  treated 
the  same  as  the  friction  resistance  component. 

Froude's  experiments  were  conducted  on  slender  bodies,  generally  thought  of 
as  those  with  L/D  >  10.  Residuary  resistance  was  defined  as  the  combined  effect  of 
wave  making  resistance  and  eddy  resistance.  At  low  speeds,  eddies  appear  to  be  a  by 
product  of  the  boundary  layer  and  therefore  part  of  frictional  resistance. 

Experimentally,  Froude  demonstrated  that  eddy  resistance  did  not  change 
proportionally  with  skin  friction  as  parallel  midbody  was  added  to  a  given  hull  form. 

He  concluded  that  eddy  resistance  was  actually  related  to  the  form  of  the  hull.  Since 
wave  making  and  eddy  resistance  are  both  related  to  form  and  are  pressure  related 
phenomenon,  they  historically  have  remained  combined  as  the  residuary  resistance 
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component'.  For  slender  bodies  the  only  significant  physicaJ  phenomena  are  friction, 
waves  and  eddy  formation.  The  process  outlined  above  was  developed  with  slender 
bodies  and  only  these  three  phenomena  in  mind. 

2.2  Frictional  Resistance 

The  problem  of  calculating  the  frictional  resistance  of  a  ship  was  first  studied 
by  considering  flat  plates.  The  essential  issue  focuses  on  the  development  of  the 
boundary  layer  along  the  length  of  the  plate  and  its  transition  from  a  laminar  to  a 
turbulent  regime.  Frictional  resistance  is  usually  expressed  as  a  coefficient  of  friction 
in  the  form  of  (1-2)  and  is  generally  plotted  as  a  function  of  the  Reynolds  Number,  Jl, 


V 


(2-4) 


where  L  is  some  characteristic  length,  usually  length  of  the  waterline  of  the  body,  o 
is  the  kinematic  viscosity  of  the  water  and  V  is  the  ship's  velocity.  The  development 
of  a  laminar  flow  equation  was  the  work  of  Blasius  (1904),  while  the  turbulent 
solution  was  formulated  independently  by  Prandtl  and  von  Karman  in  the  1920's.  In 
the  range  of  Reynolds  Number  between  5x10'  and  9xl06,  neither  the  turbulent  nor  the 
laminar  expression  yields  accurate  results.  This  region  of  uncertainty  is  known  as  the 
transition  region.  Experimentally  it  has  been  shown  that  the  transition  from  laminar  to 
turbulent  flow  occurs  gradually  across  the  length  of  a  plate  as  its  speed  increases.  This 
gives  physical  meaning  to  the  rather  large  band  of  uncertainty  between  fully  laminar 
and  fully  turbulent  flow.  The  curves  of  laminar,  transition  and  turbulent  plate  flow  are 
shown  in  Figure  (2-1). 

Model  testing  through  the  1 940's  attempted  to  draw  parallels  between  monohull 


3  Gilmer, T.C  and  Johnson. B.  Introduction 
1982,  p.217. 


to  Naval  Architecture. 


Naval  Institute  Press,  Annapolis, 


forms  and  the  known  flat  plate  solution.  The  most  accepted  flat  plate  solution  to 
frictional  resistance  adapted  to  ship  form  use  is  the  Schoenherr  Formula.  This  flat 
plate  formulation  is  frequently  seen  in  the  literature  labeled  as  the  ATTC  Line.  In  an 
attempt  to  account  for  the  effect  of  hull  form  on  boundary  layer  formation  the 
International  Towing  Tank  Conference  (ITTC)  adopted  an  interim  solution  that  was 
derived  from  a  correlation  of  previously  obtained  model  tests. 


Figure  2-1:  Laminar  and  Turbulent  Flow  Skin  Friction  Lines4 


4  Principles  of  Naval  Architecture,  op.cit  ,  p295. 
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This  formulation,  known  as  the  1TTC  1957  line  is  given  as: 


0.075 

(log109t-2)2 


(2-5) 


It  should  be  stressed  that  this  line  was  developed  from  a  correlation  of  model  data,  it 
does  not  physically  represent  the  behavior  of  a  flat  plate  and  it  was  not  obtained  by 
extrapolating  flat  plate  data  to  ship  like  forms.  Forty  years  later,  the  ITTC  has  not 
replaced  this  interim  solution  with  a  better  alternative. 

At  the  same  time  that  the  ITTC  1957  formulation  was  published,  Hughes 
proposed  a  method  for  calculating  flat  plate  frictional  resistance  given  bv: 


0.066 

(log103t-2.03)2 


(2-6) 


To  the  basic  Hughes  line,  a  form  factor  is  then  added.  The  purpose  of  this  form  factor 
is  to  help  account  for  the  differences  in  the  laminar  to  turbulent  transition  seen  by  a 
flat  plate  and  a  fuller  body.  The  attractiveness  of  Hugh’s  method  is  that  it  attempts  to 
keep  friction  and  eddy  resistance  linked.  Physical  observation  of  a  body  moving 
through  the  water  would  seem  to  make  this  reasonable.  The  method  proposed  by 
Hughes  has  not  received  the  same  acceptance  within  the  Naval  Architecture 
community'  as  the  ITTC  1957  line.  The  earliest  criticism  was  the  low  value  that 
equation  (2-6)  yields  when  compared  to  flat  plate  testing.  Additionally,  the  use  of  a 
form  factor  is  also  questionable  as  outlined  by  Oosterveld5,  quoted  below: 

"  It  has  now  become  clear  that  in  many  cases  the  wave 
pattern  resistance  is  less  than  the  residuary  resistance  estimated 


Oosterveld.  M.W.C.  ed.,  "Report  of  the  Resistance  Committee",  15th  International  Towing  Tank 
Conference.  The  Hague,  1978,  p.21 


by  assuming  a  form  factor  independent  of  Froude  number  and 
Reynolds  number.  This  will  certainly  be  the  case  whenever  bow 
wave  breaking  occurs.  One  possible  reason  for  this  is  that  the 
viscous  resistance  may  vary'  significantly  with  Froude  number 
(which  effects  the  flow  at  the  edge  of  the  boundary  layer)  as 
well  as  with  Reynolds  number,  so  that  it  is  not  well  predicted  on 
a  constant  form  factor  basis.  However  it  is  conceivable  that  this 
variation  arises  at  least  in  part  from  an  interaction  of  the 
boundary  layer  and  wake  with  the  wave  system,  so  that  energy 
losses  which  originate  with  wave  making,  and  are  primarily 
Froude  number  dependent,  may  manifest  themselves  as  an 
apparent  increase  in  the  viscous  wake  losses." 

The  frictional  lines  discussed  in  this  section  are  shown  in  Figure  (2-2). 


Figure  2-2:  Standard  Skin  Friction  Lines”’ 


This  discussion  of  frictional  resistance  has  bearing  on  the  SLICE  project  in  two 


respects.  First  there  is  the  issue  of  turbulent  flow.  The  general  practice  in  model 
testing  is  to  artificially  stimulate  turbulence.  This  is  done  so  that  the  model  and  the 
full  scale  ship  operate  in  the  same  regime.  Without  turbulence  tripping,  most  models 
would  operate  at  least  partially  in  the  transition  regime  shown  in  Figure  (2-1),  making 
data  interpretation  difficult  since  C,  is  highly  variable  in  this  regime.  Hughes7 
describes  standard  methods  of  stimulation  (usually  pins,  wires  or  sand  strips), 
methods  of  accounting  for  the  resistance  caused  by  these  devices  and  the  uncertainty 
involved  in  insuring  that  turbulence  is  actually  achieved.  The  values  of  9?  at  which  the 
SLICE  models  and  full  scale  ship  operate  are  shown  in  Table  (2-1),  A  comparison  of 
these  values  to  Figure  (2-1)  would  appear  to  demand  that  some  form  of  turbulence 
tripping  be  used  in  testing. 


Full  Scale 

l/8th  Scale 

1/1 6th  Scale 

Low  End 

Reynolds 

Number 

2.0  x  107 

(Turbulent  in 
Fig  2-1) 

2.5  x  106 

(Transition  in 
Fig  2-1) 

1.2  x  106 

(Transition  in 
Fig  2-1) 

High  End 

Reynolds 

Number 

L - 

2.1  x  108 

(Turbulent  in 
Fig  2-1) 

3.0  x  107 

(Turbulent  in 
Fig  2-1) 

1.3  x  107 

(Transition  in 
Fig  2-1) 

Table  2-1:  Operating  Reynolds  Numbers  for  SLICE 
based  on  Strut  and  Pod  Length 


Figure  (2-1)  however  shows  the  laminar  to  turbulent  transition  for  a  flat  plate.  SLICE, 
with  a  nearly  hemispherical  nose  shape  and  an  overall  L/D  *  4.5  cannot  be  likened  to 
a  flat  plate.  The  laminar  to  turbulent  transition  for  a  sphere  is  known  to  occur  in  the 


'  Hughes, G.  and  Allan,!.  "Turbulence  Stimulation  on  Ship  Models.  SNAME  Transactions.  1951, 
p28 1-314 
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region  '3?  =  4.5  x  104  -  5  x  1CF  or  two  orders  of  magnitude  lower  than  that  of  a  flat 
plate.  The  actual  transition  for  SLICE  probably  lies  somewhere  between  these  two 
extreme  boundaries.  When  viewed  in  this  light,  turbulence  stimulation  is  probably  not 
worth  the  uncertainty  involved  in  mathematically  removing  the  resistance  effects  of  the 
trip  device.  For  the  purposes  of  testing  in  this  project,  at  the  model  scale,  turbulent 
flow  will  be  assumed  even  though  turbulence  stimulation  is  not  provided. 

The  second  issue  is  to  decide  the  most  appropriate  method  to  compute  the 
frictional  component  of  SLICE  resistance.  The  ATTC  line  should  be  excluded  since  it 
is  derived  strictly  from  flat  plate  results.  The  Hughes  line  with  form  factor  corrections 
may  initially  appear  attractive  since  it  incorporates  calculated  adjustments  for  form. 

The  fullness  of  the  SLICE  body  and  the  definite  appearance  of  a  breaking  bow  wave 
may  cause  the  inaccurate  accounting  discussed  by  Oosterveld.  The  ITTC  1957  line  is 
the  best  choice.  It  is  based  on  actual  ship  and  model  results.  Even  though  there  may 
be  uncertainty  in  describing  the  details  of  laminar  to  turbulent  transition,  and  the 
impact  the  fullness  of  a  ship  shape  has  on  that  transition,  the  ITTC  1957  line  must  be 
thought  of  as  correctly  including  the  phenomenon  since  it  is  based  solely  on  test 
results.  This  is  a  fallout  of  the  model  to  ship  correlation  method  used  to  derive  the 
line.  It  is  the  only  line  discussed  that  was  developed  in  this  fashion.  Further  work  in 
this  project  will  utilize  the  ITTC  1957  line  as  the  basis  for  calculating  frictional 
resistance. 

2.3  Wave  Making  Resistance 

Wave  making  resistance  has  two  parts.  In  an  obvious  way  it  accounts  for  the 
energy  expended  creating  the  visible  wave  system  seen  as  a  body  moves  through  the 
water.  It  also  accounts  for  the  energy  expended  overcoming  drag  due  to  the  pressure 
distribution  created  by  unseparated  inviscid  flow  around  a  body.  The  literature  is 
vague  on  the  second  part  of  this  definition  of  wave  making  resistance.  The  pressure 
distribution  around  the  body  is  directly  related  to  the  form  of  the  body  but  this  effect 
is  not  the  commonly  termed  "form  drag".  "Form  drag"  is  a  Reynolds  scaled  viscous 
effect  that  is  best  thought  of  as  the  correlation  between  flat  plate  frictional  resistance 
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prediction  and  actual  ship  friction  results.  The  form  part  of  wave  making  resistance 
should  not  be  confused  with  eddy  resistance  either  which  accounts  for  the  effects  of 
separated  flow.  Since  both  parts  of  the  definition  of  wave  making  resistance  are 
pressure  effects,  this  phenomenon  will  follow  Froude  scaling. 

One  of  the  reasons  that  model  testing  is  performed  is  because  the  ability  to 
calculate  wave  making  resistance  remains  an  elusive  goal.  This  section  reviews  the 
classical  work  that  has  been  accomplished  in  this  area.  The  objective  is  to  gain  some 
insights  into  the  physical  phenomenon  and  how  it  applies  to  the  SLICE  hullform. 

The  foundation  for  a  mathematical  approach  to  wave  resistance  was  set  in  place  at 
nearly  the  same  time  that  Froude  was  developing  his  model  testing  methods.  J.H. 
Mitchell  gave  an  analytic  solution  for  wave  resistance  of  the  form8: 


f sec30 1  f  f  aexp{(-^-)sec2d(y-ixcosd)}dxdy  1 2dd 
n  V  J  J  J  V2 


(2-7) 


Where  p  is  the  density  of  water,  V  is  velocity  and  c  is  the  potential  function  that 
represents  the  body  of  interest.  The  combined  complexity  of  writing  ct  to  comply  with 
the  necessary  boundary  conditions  at  the  free  surface  as  well  as  with  the  relatively 
complex  form  of  a  monohull  and  the  difficulty  in  actually  solving  this  integral  without 
the  advantage  of  numerical  computing  techniques,  were  serious  setbacks  to  the 
popularity  of  this  approach.  The  most  significant  pursuit  of  Mitchell's  work  was 
undertaken  by  T.  Havelock  and  W.C.  Wiglev.  Wigley  was  the  first  to  successfully 
decompose  the  waves  generated  by  a  hullform  into  components  and  thus  gain  an 
understanding  of  how  they  are  formed.  This  was  first  achieved  with  what  is  now 
called  the  Wigley  hull,  a  wedge  followed  by  parallel  midbodv  followed  by  an  identical 
wedge  pointed  opposite  to  the  bow.  Eventually  Wigley  refined  the  hull  to  a  smooth 


Newman.  J.N.,  Marine  Hydrodynamics.  The  MIT  Press,  Cambridge.  1989.  p.282. 
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torm.  He  discovered  that  there  are  actually  four  waves  produced  by  this  shape.  The 
first  is  produced  by  the  leading  edge  of  the  hull,  the  next  two  are  produced  by  the  two 
shoulders  formed  by  the  intersection  of  the  fore  and  aft  wedge  with  the  parallel 
midbody  and  the  final  one  is  produced  by  the  trailing  edge  of  the  aft  wedge.  This 
decomposition  held  true  for  hullforms  with  continuous  lines  as  well.  Figure  (2-3)  is 
reproduced  from  Wigley's  writings9  which  graphically  demonstrates  this  additive 
effect. 

This  discussion  is  useful  to  the  SLICE  program  because  it  gives  insights  to  an 
approach  at  hullform  optimization.  Since  the  amplitude  of  the  wave  system  generated 
is  related  to  the  form  of  the  body,  maintaining  lines  that  are  reasonably  fine  should 
help  to  minimize  the  hump  of  the  wave  making  curve.  As  a  approximation  without 
actual  calculations,  the  shape  of  the  struts  and  pods  of  the  SLICE  form  will  be 
adjusted  in  an  effort  to  find  a  reasonable  compromise  between  frictional  and  wave 
making  resistance.  This  is  different  from  the  LMSC  approach  which  appears  to  have 
concentrated  on  the  operating  range  of  >  0.9  with  effort  focused  on  minimizing 
frictional  resistance. 

The  SLICE  hullform  can  be  decomposed  into  eight  different  bodies,  four  struts 
and  four  pods.  Each  of  these  component  bodies  has  a  form  that  is  reasonably  similar 
to  the  hullforms  tested  by  Wiglev.  If  each  of  these  bodies  can  be  thought  of  as 
producing  four  waves  that  combine  into  a  single  wave  system  for  the  body,  then  an 
approach  to  optimizing  wave  making  resistance  can  be  created  by  adjusting  body 
spacing  once  form  issues  are  decided.  A  three  tier  system  of  looking  at  this  problem  is 
proposed: 

1)  First  consider  the  interaction  between  a  strut  and  its  attached  pod.  This  step 
can  be  thought  of  as  creating  an  oversized  bulbous  bow  for  the  strut. 

2)  Next  consider  the  effect  of  longitudinally  aligned  strut  pod 
combinations.  This  will  probably  have  a  diminished  effectiveness  compared  to  step  1 

£ 

WiglevAV.C.S..  "A  Comparison  of  Experiment  and  Calculated  Wave-Profiles  and  Wave- 
Resistances  for  a  form  having  Parabolic  Waterlines",  Proceedings  of  the  Royal  Society.  Series  A, 

London,  1934.  pp  144-159. 
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Figure  2-3:  Calculated  Wave  Profile  Analysis 


since  the  amplitude  of  the  forward  system  will  not  equal  that  of  the  aft  system  once 
they  interact. 

3)  Finally  consider  the  interactions  of  the  transversely  and  diagonally  aligned 
strut  pod  combinations.  Since  these  are  not  in  the  principle  direction  of  motion,  it  is 
expected  that  the  effectiveness  of  these  interactions  will  again  be  diminished  from  the 
previous  step. 

2.4  Eddy  Resistance 

As  previously  noted,  the  literature  is  noticeably  vague  in  defining  both  form 
and  eddy  resistance.  They  are  frequently  used  together  but  are  not  used 
interchangeably.  Articles  written  about  resistance  prediction  usually  account  for  eddy 
resistance  as  an  input  to  residuary  resistance.  The  research  for  this  thesis  uncovered  no 
direct  usage  of  form  resistance  in  prediction  schemes.  As  discussed  in  Chapter  2-2, 
this  thesis  considers  form  drag  to  be  the  increase  in  frictional  resistance  felt  by  a  ship 
shape  compared  to  an  equivalent  flat  plate.  The  1TTC  1957  line  directly  accounts  for 
form  drag  since  it  is  based  on  experimental  data  from  monohull  forms.  Other  schemes 
such  as  the  ATTC  line  and  Hughes  Line  require  some  form  compensation  for  accurate 
use  since  they  are  based  on  flat  plate  resistance. 

Eddy  resistance  on  the  other  hand  has  physical  origins.  In  this  thesis  it  is  taken 
to  mean  the  result  of  separation  effects  from  discontinuities  in  the  hull  form.  These 
can  result  from  actual  breaks  in  the  ships  lines  such  as  the  intersection  of  the  strut  and 
pod  at  both  the  leading  and  trailing  edge,  intersection  of  appendages  and  the  main 
body,  or  radii  of  curvature  which  are  too  hard  for  the  flow  to  follow.  As  discussed  in 
chapter  2-1,  eddy  resistance  is  a  pressure  effect  and  therefore  should  be  scaled  as  part 
of  the  residuary  resistance  component. 

Following  the  classic  approach,  the  minor  swirling  seen  along  the  sides  of  a 
ship  are  associated  with  viscous  effects  of  well  behaved  flow  and  are  accounted  for  in 
the  frictional  component10.  Significant  swirling  that  can  be  associated  with  flow 


10  Principles  of  Naval  Architecture,  op.cit.  p. 3 1 2. 


27 


separation,  such  as  that  seen  at  the  transom,  describes  the  type  of  effect  accounted  for 
by  eddy  resistance.  There  is  some  obvious  room  for  interpretation  of  cause  and  effects 
in  these  definitions  of  form  and  eddy  drag.  This  highlights  the  importance  of  observing 
physical  behavior  through  model  testing. 

2.5  Spray  Resistance 

The  formation  of  spray  and  the  impact  it  has  on  total  drag  has  historically  been 
assumed  insignificant  to  monohull  testing.  Simple  observation  of  full  and  model  scale 
monohulls  operating  at  speed  validates  this  assumption.  Although  spray  is  generated, 
the  region  of  generation  is  highly  localized  and  the  quantity  of  spray  generated  is 
usually  small.  With  the  increased  interest  in  advanced  hullforms  such  as  hydrofoils 
and  planing  craft,  neglecting  spray  was  recognized  as  a  poor  assumption  but  the  body 
of  available  knowledge  regarding  the  mechanism  of  sprav  generation  and  the 
appropriate  method  of  scaling  spray  results  is  relatively  small. 

When  spray  becomes  a  major  concern,  the  one  equation,  one  unknown  analogy 
used  in  section  2-1  becomes  one  equation  and  two  unknowns.  The  typical  solution  to 
this  problem  has  been  to  test  at  very  high  speeds  (SF >  3.0)  to  eliminate  wave  making. 
This  allows  spray  drag  to  be  analyzed  separately  but  is  only  a  partial  answer.  A  gap 
exists  at  the  lower,  more  practical  values  of  In  1971,  Chapman"  conducted  a 
significant  test  program  on  strut  forms  with  parabolic  leading  and  trailing  edges.  He 
concluded  the  following: 

1 .  Spray  drag  is  partially  dependent  on  both  strut  form  and  thickness.  Blunt 
leading  edges  produce  more  spray. 

2.  The  spray  mass  flow  rate  of  spray  depends  on  the  thickness  to  chord  ratio 
(t/c)  and  on  .  Blunt  bodies  send  the  spray  up  at  higher  angles  causing  greater 
wetting. 

3.  Skin  friction  due  to  wetting  is  the  primary  source  of  spray  drag. 


Chapman.  R.B  .  Spray  Drag  of  Surface  Piercing  Struts.  Naval  Undersea  Research  and 
Development  Center  Report  AD  730710.  1971. 
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4.  The  momentum  in  the  spray  sheet  is  sufficient  to  create  a  large  pressure 
drag  on  any  object  it  strikes. 

Chapman's  work  was  expanded  by  Zhu  Bing-quan  in  198612.  A  total  of  twelve 
struts  were  tested,  all  with  parabolic  trailing  and  leading  edges.  The  objective  was  to 
correlate  strut  performance  to  the  parameters  c  (chord  length),  l  (  maximum  thickness) 
and  x  (  longitudinal  position  of  the  maximum  thickness).  The  relationship  derived  is: 


R=pV2x2C. 


(2-8) 


and  Cs  is  given  as: 


C =0.068-  -0.555(-)2+0.696(-)3  +0.083(-)2^  (2-9) 

t  t  t  t 

The  research  also  showed  that  the  spray  drag  was  minimized  when  /  .y  =  0.313. 

This  has  applications  to  the  SLICE  analysis.  The  stmt  form  selected  by  LMSC 
is  contrary  to  that  recommended  by  the  published  literature  for  minimizing  spray  drag. 
The  apparent  design  goal  of  minimizing  frictional  resistance  fails  to  consider  spray 
effects.  In  deciding  the  extent  to  which  the  spray  effects  will  scale  from  model  size  to 
lull  scale,  the  dependence  of  (2-9)  on  dZ" implies  that  spray  drag  will  follow  Froude 
scaling.  The  conclusion  reached  by  Chapman  and  Zhu  Bing-quan  that  the  initiation  of 
spray  is  highly  form  dependent  also  supports  this  assertion. 

The  following  is  recommended  when  dealing  with  spray  effects.  Choose  stmt 
shapes  that  are  known  to  minimize  spray  generation.  Unless  there  is  a  significant 
structural  or  interior  arrangements  reason  to  do  otherwise,  stmts  should  have  parabolic 
leading  and  trailing  edges  and  t  x  should  be  reasonably  close  to  0.313.  For  the  ATD, 


^  Zhu  Bing-quan,  Ge  Wei-zhen.  "Experimental  Study  on  Spray  Performance  of  Surface  Piercing 
Struts,  Journal  of  Hydrodynamics,  ser  2  no  4,  1990,  pp  91-98. 
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this  will  help  minimize  spray  effects  and  allow  the  research  efforts  to  focus  on  the 
area  of  principle  interest,  the  wave  making  resistance.  Any  spray  that  is  observed  in 
testing  should  be  treated  as  part  of  the  residuary  resistance.  If  the  tentative  conclusion 
that  spray  follows  Froude  scaling  is  wrong,  then  this  process  will  over  predict  the  total 
resistance.  Considering  the  uniqueness  of  the  SLICE  hullform,  a  conservative  error  of 
this  nature  is  well  justified  for  the  first  full  scale  ship.  Design  work  in  this  project  will 
follow  these  conclusions. 

2.6  Appendage  Resistance 

The  effect  of  appendages  on  total  drag  has  historically  been  assumed  small  in 
monohull  testing.  This  assumption  is  based  on  the  fact  that  the  overall  size  of  the 
appendages  is  usually  one  or  two  orders  of  magnitude  less  than  the  main  hull.  For 
comparative  testing  of  hulls,  it  is  reasonable  to  ignore  the  appendages.  For  obtaining 
exact  resistance  predictions,  appendages  are  usually  included  in  the  final  testing  and 
stripping  tests  have  been  devised  to  help  account  for  their  effects.  In  stripping  tests, 
the  model  is  towed  first  with  all  the  appendages  on.  In  subsequent  tests,  the 
appendages  are  removed  one  at  a  time.  In  this  manner  the  impact  of  each  appendage 
can  be  assessed.  This  approach  assumes  that  the  interaction  effects  of  the  appendages 
with  each  other  and  with  the  main  hull  are  negligible. 

For  SLICE,  appendages  occupy  a  larger  percentage  of  total  volume  compared 
to  monohulls.  The  flow  around  SLICE  is  more  complex  and  therefore  interactions  may 
be  more  pronounced.  A  full  stripping  test  is  warranted  in  the  case  of  SLICE.  Since 
this  project  is  focused  a  comparison  of  alternate  hull  forms,  appendages  will  be 
excluded,  leaving  the  stripping  test  to  future  work. 

2.7  The  SLICE  Data  Scaling  Method 

Based  on  the  above  discussion,  the  traditional  method  of  data  scaling  outlined 
in  section  2-1  is  recommended  for  the  SLICE  project  and  will  be  used  in  this  research. 
Calculating  the  frictional  resistance  using  the  ITTC  1957  line  will  yield  a  moderately 
conservative  answer.  Although  the  ITTC  1957  line  was  derived  from  monohull  tests,  it 
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is  the  only  method  that  naturally  accounts  for  the  Froude  Number  dependence  of 
boundary  layer  development.  All  other  phenomena  (wave  making,  eddy  formation  and 
spray  development)  will  be  treated  as  part  of  residuary  resistance  and  scaled 
accordingly. 

It  is  important  to  recognize  is  that  this  technique  of  data  scaling  is  considered 
reliable  due  to  the  enormous  volume  of  model  testing  that  has  been  performed  on 
numerous  series  of  monohull  models.  The  validation  of  this  process  for  SLICE  will 
only  come  when  the  full  scale  ship  is  built  and  tested.  This  is  the  real  value  of  the 
ATD.  Building  and  properly  testing  the  full  scale  ship  will  significantly  contribute  to 
the  Naval  Architecture  community's  understanding  of  the  robustness  of  current  testing 
and  scaling  techniques. 

2.8  A  Numerical  Approach 

Model  testing  is  a  time  consuming,  detail  intensive  process  that  requires 
significant  capital  expenditure.  For  a  new  hull  form,  some  degree  of  model  testing  is 
mandatory  in  the  early  stages  of  design  so  observations  of  physical  behavior  can  be 
made.  This  is  critical  to  validating  assumptions,  but  from  a  cost  savings  perspective, 
the  ability  to  numerically  model  the  ship's  behavior  is  attractive.  Small  changes  to  the 
hull  can  be  made  and  tested  without  the  expense  of  building  a  new  model.  If  the 
numerical  model  can  be  kept  mathematically  simple,  comparative  results  can  be 
obtained  quickly  and  cheaply.  There  must  be  a  degree  of  confidence  however  that  the 
computer  model  fully  represents  actual  resistance. 

The  physical  phenomena  of  interest  are  friction,  wave,  spray  and  eddy 
resistance.  Friction  and  spray  resistance  can  be  calculated  using  (2-5),  (2-8)  and  (2-9). 
If  there  is  no  significant  separation  or  vortex  phenomena,  it  is  reasonable  to  neglect 
eddy  resistance  (only  observation  can  answer  this)  and  a  prediction  for  full  scale 
resistance  can  be  obtained  if  a  method  of  solving  (2-7)  can  be  found.  This  is 
significant  to  the  SLICE  discussion  since  it  may  provide  a  simple  and  fundamental 
tool  for  predicting  and  minimizing  the  wave  making  hump  magnitude. 
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A  solution  to  Mitchell's  integral  was  outlined  by  Lunde1,1  in  1951.  The  process 
was  further  applied  to  a  ship  called  the  Wave  Cancellation  Multi-hull  Ship  Concept  by 
Wilson  and  Hsu4.  They  express  the  total  wave  making  resistance  as  the  summation  of 
body  and  interaction  terms  as  shown: 


K 


2  N  N  N  N 


(2-10) 


w  here  N  represents  the  number  of  bodies  the  overall  hull  can  be  divided  into  and: 


TL 
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o  o 
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Lunde.  J.K..  “On  the  Linearized  Theory  of  Wave  Resistance  for  Displacement  Ships  in  Steady 
and  Accelerated  Motion".  SNAME  Transactions.  Vol  59,  1951,  pp.  25-76. 

14  Wilson.  MB.  and  Hsu.  C.C.,  "Wave  Cancellation  Multihull  Ship  Concept",  Presented  at  the 
Intersociety  High  Performance  Marine  Vehicle  Conference.  Arlington,  Va..  24-27  June  1992. 


3  2 


K  =  X 
V2 


(2-1  Od) 


In  these  expressions,  the  like  product  terms  are  the  body  terms  and  the  cross  product 
terms  are  the  interaction  terms. 

The  underlying  assumptions  to  Lunde's  approach  are  : 

1.  The  fluid  is  incompressible,  irrotational  and  inviscid. 

2.  The  wave  height  is  small  compared  to  wave  length. 

3.  The  ship  motion  can  be  reasonably  represented  by  a  fine  distribution  of 
sources  and  sinks  along  the  centerline  of  the  body. 

Prior  observation  of  SLICE  models  in  the  towing  tank  indicates  that  conditions 
2  and  o  may  not  be  satisfied  ideally.  Each  of  these  phenomenon  must  be  watched 
closely  in  subsequent  tests  and  considered  in  the  interpretation  of  results.  Havelock15 
demonstrated  good  correlation  using  this  technique  with  hull  forms  having  L/B  as 
small  as  10.6.  Data  presented  by  Newman16  indicates  that  there  is  better  correlation 
between  model  test  data  and  thin  ship  theory  at  higher  values  of  ZF.  SLICE  is  a  full 
form  with  L/D  =  4. 5-5.0  but  operates  at  &=  0.9-1. 2.  Although  SLICE  is  not  an  exact 
fit  for  thin  ship  theory',  the  approach  may  be  reasonable. 


L'  Havelock,  T  H  ,  "Wave  Resistance  Theory  and  its  Application  to  Ship  Problems",  SNAME 
Transactions.  1951,  pp.  13-24. 

!6  Newman.  J.N.  ,  op.cit..  p282. 
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Chapter  3  Model  Testing 


As  part  of  this  research,  three  models  were  designed,  built  and  tested.  Appendix  A 
provides  detailed  lessons  learned  from  this  testing  experience  which  may  be  useful  to 
other  students  attempting  to  undertake  similar  model  testing  projects. 

3.1  Model  Test  Objectives 

There  were  three  objectives  underlying  the  model  testing  phase  of  this  project: 

1.  Expand  the  quantity  of  available  model  scale  resistance  data  for  the 
baseline  SLICE  configuration. 

2.  Visually  observe  model  performance  in  sufficient  detail  to 
understand  the  physical  behavior  of  this  hullform.  As  discussed  in  chapter  2,  this  is 
important  for  justifying  the  scaling  technique  selected. 

3.  Test  operationally  equivalent  hullforms  to  help  assess  baseline 
SLICE  advantages  and  disadvantages. 

The  models  were  tested  at  the  United  States  Naval  Academy's  380'  tow  tank 
facility  during  the  week  of  6-10  March  1995. 

3.2  Model  Descriptions 

The  three  models  tested  as  part  of  this  research  were  designed  to  meet  the 
same  operational  constraints.  This  was  done  to  ensure  that  comparisons  made  at  the 
end  of  the  testing  could  be  made  on  the  basis  of  similar  ships.  All  models  when  scaled 
to  full  scale  have  a  displacement  of  177  Lton,  a  pod  diameter  of  7-8  feet  and  a  strut 
beam  of  3.25  feet.  These  dimensions  are  the  working  dimensions  of  the  ATD  vessel, 
driven  by  requirements  other  than  hydrodynamic  performance.  The  pod  diameter  is 
controlled  by  engine  size  requirements  and  the  strut  beam  is  controlled  by  access 
trunks  into  the  pods.  With  these  requirements  imposed  on  all  three  models,  they  are 
considered  equivalent  hulls. 

The  three  models  tested  were  designated  M-l.  M-2  and  M-3.  Detailed  tables 
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of  hull  offsets  and  design  calculations  for  each  model  are  contained  in  Appendices  B. 
C  and  D  respectively.  Individual  descriptions  of  each  model  follow. 

MODEL  M-l  BASELINE: 

This  model  utilizes  LMSC's  offsets  as  tested  in  December  1994  with  the 
exception  that  the  fillets  at  the  strut/pod  interface  and  the  haunch  structures  are 
excluded.  The  fillets  were  excluded  to  help  simplify  model  building.  This  modification 
is  not  considered  detrimental  to  the  testing  outcome.  The  combination  of  the  pod,  strut 
and  fillet  in  the  LMSC  model  creates  a  effective  L/D  of  3.9.  the  omission  of  the  fillets 
increases  this  ratio  to  a  more  favorable  4.2.  This  modification  should  result  in  slightly 
lower  resistance  values.  The  bow  wedges  installed  in  the  December  LMSC  model  are 
replicated  in  this  model.  The  purpose  of  this  model  is  twofold.  First,  it  provides  a 
means  of  correlating  the  data  obtained  in  these  tests  with  the  data  obtained  by  LMSC 
in  December  1994.  This  is  considered  important  since  the  tests  were  conducted  in  two 
different  facilities,  under  different  conditions  with  different  scale  models.  The  ability 
to  correlate  data  between  Model  M-l  and  the  LMSC  model  will  be  valuable  in 
demonstrating  repeatability  of  results.  The  second  use  of  Model  M-l  is  to  expand  the 
available  range  of  data  for  the  baseline  SLICE.  Model  tests  were  conducted  at  model 
equivalent  speeds  through  40  knots  to  ensure  that  the  entire  operating  range  was 
properly  represented  by  data. 

Table  3-1  provides  the  principle  dimensions  of  Model  M-l.  Profile  and  top 
view  sketches  of  the  forward  and  aft  strut/pod  combinations  are  in  Figure  3-1.  In  this 
drawing  the  alignment  of  the  struts  with  the  corresponding  pods  is  correctly 
represented.  The  spacing  betw'een  the  forward  and  aft  combination  is  not  preserved 
however  in  the  interest  of  fitting  the  drawing  onto  one  page.  This  model  was 
constructed  with  29  Ibf  of  reserve  buoyancy.  The  struts  are  manufactured  from  balsa 
wood  with  a  fiberglass  shell,  the  pods  are  molded  fiberglass. 
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Scale  Factor 


17.26:1 


Displacement 

75  Lbf 

Beam 

36.15  inches 

Draft 

8.34  inches 

Pod  Diameter 

5.5625  inches 

Fwd  Pod  Length 

22.9  inches 

Aft  Pod  Length 

25  inches 

Fwd  Pod  Separation  (CL  to  CL) 

23.4  inches 

Aft  Pod  Separation  (CL  to  CL) 

30.25  inches 

Strut  Length 

16.68  inches 

Stmt  Max  Beam 

2.26  inches 

Table  3-1:  Model  M-l  Baseline 
Principle  Dimensions 


Figure  3-1:  Model  M-l  Forward  and  Aft  Strut/Pod 

(Note:  In  this  drawing  the  fwd-aft  pod  spacing  is  not  preserved) 
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MODEL  M-2  RATIONAL  DESIGN: 


The  purpose  of  this  model  is  to  explore  the  effects  on  total  resistance  of  using 
finer  offsets  for  the  struts  and  pods  .  To  avoid  changing  too  many  variables  in  a  single 
test,  this  model  maintains  all  the  gross  parameters  of  the  Model  M-l  Baseline  (ie: 
Strut/Pod  alignment,  fore/aft  spacing  and  centerline  to  centerline  spacing).  The 
expectation  for  this  design  is  that  by  comparing  this  model's  results  with  that  of  M-l, 
some  conclusions  can  be  drawn  about  the  relative  magnitude  of  both  friction  and  wave 
making  resistance  through  the  range  of  speeds  tested.  The  anticipated  result  from  this 
model  is  that  the  amplitude  of  the  wave  making  hump  will  be  diminished. 

The  offsets  for  the  pods  are  adapted  from  Senes  58  data17.  This  series  was 
selected  since  it  represents  a  concerted  effort  to  optimize  submarine  hull  forms.  The 
model  selected  from  the  series  was  Model  4155.  This  model  was  selected  based  on  its 
length  to  diameter  (L/D)  ratio  of  5  which  is  similar  but  finer  than  the  baseline  SLICE 
L/D  of  4.2.  It  should  be  noted  that  within  the  Series  58  results.  Model  4155  falls 
outside  the  range  of  optimum  L  D  which  was  found  to  be  between  6  and  8.  Model 
4155  was  selected  despite  this  drawback  in  order  to  satisfy  the  operational  requirement 
of  creating  a  full  scale  177  Lton  vessel  with  pod  diameters  of  7-8  feet.  These  criteria 
could  not  be  satisfied  with  a  senes  58  hull  with  a  L/D  of  6-8. 

A  trait  of  the  Series  58  hullforms  is  that  they  have  no  parallel  midbody. 
Although  hydrodvnamically  this  is  superior,  it  would  significantly  complicate  the 
building  of  the  full  scale  ship.  The  Model  M-2  therefore  incorporates  a  modification  to 
Model  4155  by  replacing  some  of  the  midships  shape  with  parallel  midbody. 

The  strut  offsets  utilize  a  symmetric  parabolic  leading  and  trailing  edge.  The 
general  parameters  of  this  strut  are  outlined  in  Chapter  2-5  but  are  also  constrained  by 
the  design  requirement  that  the  full  scale  maximum  beam  must  be  at  least  3.25  feet. 
The  calculations  used  to  develop  the  stmt  and  pod  offsets  and  a  table  of  the  actual 

17  Gertler.  M.  "Resistance  Experiments  on  a  Systematic  Series  of  Streamlined  Bodies  of 
Revolution  for  Application  to  the  Design  of  High  Speed  Submarines",  David  Taylor  Model  Basin  Report 
C-297.  unclassified  and  approved  for  public  release  18  Dec  1972,  April  1950. 
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offsets  used  are  contained  in  Appendix  C. 

Model  M-2  represents  a  partial  optimization.  The  offsets  for  the  struts  and  pods 
are  based  on  research  conducted  in  optimizing  these  types  of  shapes.  The  optimized 
values  are  then  adjusted  to  conform  to  design  requirements.  Model  M-2  make  no 
attempt  to  optimize  the  placement  of  the  strut  pod  combinations  relative  to  one 
another.  Table  3-2  provides  the  principle  dimensions  of  Model  M-2.  Profile  and  top 
view  sketches  of  a  the  forward  and  aft  strut/pod  combinations  are  at  Figure  3-2.  This 
model  was  constructed  with  30  lbf  of  reserve  buoyancy.  The  struts  are  manufactured 
from  balsa  wood  with  a  fiberglass  shell,  the  pods  are  molded  fiberglass. 

MODEL  M-3  OPTIMIZED  HIGH  SPEED  SWATH: 

The  purpose  of  this  model  is  to  compare  SLICE  performance  to  a  more 
conventional  SWATH.  The  offsets  used  for  the  high  speed  SWATH  were  reported  by 
McGregor18.  The  demihull  diameter  was  forced  to  meet  the  7-8  foot  full  scale 
requirement,  the  strut  maximum  beam  was  forced  to  meet  the  3.25  foot  full  scale 
requirement  and  the  overall  displacement  was  forced  to  meet  the  177  Lton  full  scale 
requirement  established  earlier.  In  order  to  meet  these  requirements,  the  demihull  L/D 
was  established  at  12.86.  This  corresponds  to  a  full  scale  demihull  length  of  90  feet, 
about  15  feet  shorter  than  the  Length  Overall  (LOA)  of  the  proposed  ATD  SLICE. 

The  consequence  of  this  is  that  although  Model  M-3  is  hvdrodynamically  equivalent  to 
Model  M-l,  depending  on  the  extent  to  which  the  box  structure  is  cantilevered  forward 
and  aft  of  the  struts,  the  SWATH  will  end  up  with  as  much  as  600  square  feet  less 
arrangable  deck  space.  This  is  an  11%  decrease  from  the  SLICE  ATD  vessel. 

A  tandem  strut  arrangement  was  selected  in  order  to  minimize  the  displacement 
occupied  by  the  struts.  This  was  necessary  in  order  to  maximize  the  demihull 
displacement  and  maintain  a  reasonable  demihull  length.  As  with  model  M-2,  the 


18  McGregor.  R.C.  and  Chun,  H.H.,  "On  the  Potential  of  SWATH  Ships  for  Very  High  Speed 
Operation",  Fast  91.  pp.  491-506. 
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strut  offsets  are  based  on  the  guidance  provided  in  Chapter  2-5.  In  his  SWATH 
design,  McGregor  provided  guidance  for  the  optimum  placement  of  the  struts  in 
relation  to  the  pods.  This  placement  was  utilized  in  Model  M-3.  The  calculations  used 
to  develop  the  strut  and  pod  offsets  and  a  table  of  offsets  are  contained  in  Appendix 
D.  Table  3-3  provides  the  principle  dimensions  of  Model  M-3.  A  profile  and  top  view 
sketch  of  a  the  demihull  is  at  Figure  3-3.  The  struts  of  this  model  are  manufactured 
from  balsa  wood  with  a  fiberglass  shell,  the  demihulls  are  hollowed  PVC  rod,  shaped 
to  the  proper  offsets.  There  is  10  lbf  of  reserve  buoyancy. 

The  reason  the  three  models  were  not  all  built  to  the  same  scale  factor  stems 
from  the  actual  building  process.  There  is  a  slight  variation  in  the  full  scale  diameter 
of  the  three  underwater  hull  diameters.  In  full  scale.  M-l  is  8  feet.  M-2  is  7.74  feet 
and  M-3  is  7  feet.  PVC  pipe  was  used  to  mold  parallel  midbodv  sections.  This 
material  is  relatively  inexpensive  and  comes  in  a  standard  pipe  outer  diameter 
dimension  of  5.5625  inches.  This  dimension  became  the  driver  for  determining  the 
model  scale  factor.  Since  the  three  models  have  different  full  scale  diameters,  the  scale 
factors  vary  accordingly. 
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Scale  Factor 


16.68:1 


Displacement 


Pod  Diameter 


Fwd  Pod  Leneth 


Aft  Pod  Leneth 


Fwd  Pod  Separation  (CL  to  CL) 


Aft  Pod  Separation  (CL  to  CL) 


Strut  Leneth 


Strut  Max  Beam 


83.3  lbf 


>7.41  inches 


8.63  inches 


5.5625  inches 


27.8  inches 


27.8  inches 


24.17  inches 


;6.26  inches 


’.25  inches 


2.34  inches 


Table  3-2:  Model  M-2  Modified  Baseline 
Principle  Dimensions 


Figure  3-2:  Model  M-2  Forward  and  Aft  Strut/Pod 

(Note:  !n  this  drawing  the  fwd-aft  pod  spacing  is  not  preserved) 
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Scale  Factor 

15.10:1 

Displacement 

115.1  lbf 

Beam 

39.93  inches 

Draft 

11.13  inches 

Demihull  Diameter 

5.5625  inches 

Demihull  Length 

71.5  inches 

Strut  Length 

19.10  inches 

Strut  Max  Beam 

2.58  inches 

Table  3-3:  Model  M-3  McGregor's  High  Speed  SWATH 
Principle  Dimensions 


Figure  3-3:  Model  M-3  Demihull  and  Strut 
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3.3  Test  Description 

A  copy  of  the  test  plan  developed  for  the  Naval  Academy  Tow  Tank  tests  is 
included  as  Appendix  E.  This  enclosure  explains  the  planned  tests  in  detail. 

Two  basic  tests  were  planned  for  the  models.  The  first  measures  resistance  and 
lift  with  the  models  fixed  in  pitch  and  heave.  This  test  replicates  the  test  configuration 
used  by  LMSC  thereby  allowing  for  correlation  with  the  LMSC  data.  Data  collection 
runs  consisted  of  the  model  held  at  its  design  draft  with  no  trim.  At  the  slow  end  of 
the  speed  spectrum,  one  speed  change  was  allowed  during  data  collection.  At  higher 
speeds  only  one  speed  was  used  per  test  run.  Data  was  collected  at  full  scale 
equivalent  speeds  from  3-40  knots  with  the  following  model  configurations: 

1)  Model  M-l  Baseline  in  the  LMSC  configuration. 

2)  Model  M-l  Baseline  with  the  forward  and  aft  struts  in  line. 

3)  Model  M-l  Baseline  in  the  LMSC  configuration  at  varying  drafts. 

4)  Model  M-2  Modified  Baseline  in  the  LMSC  configuration. 

5)  Model  M-3  McGregor's  SWATH. 

The  second  test,  quasi  free  in  heave,  was  not  conducted.  The  intent  of  this  test 
was  to  replicate  the  effects  of  control  surfaces  through  a  counterweight  system  with 
the  hope  that  information  could  be  gained  regarding  the  necessary  sensitivity  of  the 
control  system.  Conceptually,  the  counterweight  would  provide  the  same  lift  as  the 
control  surfaces  without  added  model  building  complexity.  The  model  would  then  be 
free  to  respond  to  small  perturbations  and  its  stability  qualitatively  assessed.  This  test 
was  considered  a  second  priority  to  the  project  and  time  in  the  test  facility  ran  out 
before  it  could  be  accomplished.  It  is  considered  a  worthwhile  area  of  further  study. 

The  test  plan  indicates  that  a  fourth  model,  the  Midfoil,  was  tested.  This  model 
was  built  and  tested  at  the  request  of  ONR  but  is  not  considered  part  of  this  research. 

It  is  not  reported  in  this  thesis. 
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3.4  Physical  Observations 

The  importance  of  observing  the  behavior  of  the  models  can  not  be  overstated. 
A  videotape  of  the  model  tests  is  included  as  part  of  this  thesis.  Shown  in  the  tape  is  a 
port  beam  perspective,  both  above  and  underwater  as  well  as  a  stem  view  of  all  the 
models  tested.  The  video  is  organized  by  speed  increment,  showing  each  model  first  at 
the  peak  of  the  wave  making  hump  (10  and  12  knot  full  scale  equivalent)  then  in  the 
post  hump  valley  (16  knot  full  scale  equivalent)  and  finally  at  three  high  speed 
increments  (25,  30  and  35  knot  full  scale  equivalents).  The  following  qualitative 
observations  are  made  regarding  the  performance  of  the  models. 

PROPELLER  SUBMERGENCE 

The  ATD  vessel  will  have  propellers  mounted  on  the  aft  end  of  the  forward 
pods.  After  the  August  LMSC  model  tests,  there  was  significant  concern  that  the  wave 
system  created  in  the  hump  was  large  enough  to  expose  the  propeller.  This  evaluation 
was  made  by  observing  the  wake  from  above  the  free  surface  but  was  difficult  to 
assess  properly  due  to  spray  and  distortion.  With  M-l,  the  underwater  camera  used  in 
this  project  reveals  that  there  is  better  submergence  at  the  propeller  than  apparent 
when  viewed  from  above. 

HUMP  SPEED  PERFORMANCE 

It  was  expected  that  M-2  with  its  finer  lines  would  exhibit  less  wave  making 
than  the  very  bluff  M-l  in  the  hump  speed  range.  In  fact,  the  surprising  result  was 
that  the  wakes  produced  were  nearly  identical.  The  wake  of  M-3  at  the  same  speed 
was  much  smaller  than  either  M-l  or  M-2.  The  implications  of  this  observation  are 
unclear  without  examining  actual  data.  The  similarity  between  the  M-l  and  M-2  wakes 
would  suggest  that  the  form  of  the  bodies,  both  strut  and  pod,  is  not  a  controlling 
factor  in  hump  speed  resistance.  The  vast  difference  between  M-l /M-2  and  M-3 
suggests  that  the  underwater  body  plays  a  role  in  wave  cancellation. 

The  correct  conclusion  may  simply  be  that  this  is  a  demonstration  of  the  flow 
differences  of  bluff  bodies  (M-l  L/D  =  4.2,  M-2  L/D  =  5.0)  and  more  slender  bodies 
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(M-3  L/D  =  12.9)  or  there  may  be  some  more  complicated  mechanism  of  wave 
cancellation  that  needs  to  be  explored. 

HIGH  SPEED  PERFORMANCE 

In  the  full  scale  range  of  25-40  knots,  there  is  a  significant  difference  in  the 
performance  of  M-l  and  M-2.  At  these  speeds,  M-l  creates  a  very  large  vortex  at  the 
aft  intersection  of  the  strut  and  the  pod.  Above  the  surface,  this  is  manifest  as  a 
rooster  tail.  As  the  speed  increases,  the  inception  point  of  this  vortex  travels  further  aft 
along  the  pod  tail  and  at  the  full  scale  equivalent  of  40  knots  almost  reaches  the  end 
of  the  pod.  This  effect  is  absent  in  M-2  and  M-3. 

This  effect  does  not  appear  to  be  caused  by  the  relatively  small  scale  of  the 
models.  Reviewing  footage  of  the  LMSC  tests,  the  same  phenomenon  is  apparent.  It  is 
more  pronounced  in  the  August  tests  than  the  December  tests,  but  it  is  difficult  to  tell 
if  that  can  be  attributed  to  the  design  changes  LMSC  made.  There  is  not  enough 
quality  video  footage  of  the  December  tests  to  draw  this  conclusion.  It  is  possible  that 
the  August  and  December  tests  did  not  allow  a  full  steady  state  condition  to  be 
established.  In  both  the  August  and  December  tests,  speed  increments  were  maintained 
on  average  for  10  seconds  measured  from  the  start  of  a  speed  change  to  the  start  of 
the  subsequent  speed  change.  As  observed  in  the  March  tests,  this  amount  of  time  is 
barely  adequate  to  establish  steady  state  conditions.  When  accelerated  from  rest,  the 
March  data  shows  that  forces  did  not  stabilize  until  8-10  seconds  after  the  steady  state 
speed  was  achieved.  The  LMSC  testing  may  have  stabilized  sooner  since  the 
acceleration  was  only  between  speed  increments  and  not  from  rest.  Nonetheless,  the 
data  collection  period  was  extremely  short.  The  inception  speed  of  the  rooster  tail 
appears  to  be  comparable  in  both  the  LMSC  tests  and  this  project's  tests.  This 
indicates  that  the  effect  follows  Froude  scaling. 

After  the  August  LMSC  tests,  much  attention  was  focused  on  scaling  spray 
effects.  This  was  in  response  to  the  extreme  wake  created  at  high  speed.  LMSC's 
position  was  that  spray  effects  could  be  ignored  in  full  scale  predictions  and  designed 
out  using  spray  rails.  The  underwater  observation  of  M-l  demonstrates  that  much  of 
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what  was  called  spray  in  August  was  actually  produced  by  an  entirely  different 
mechanism.  The  vortex  that  creates  the  rooster  tail  is  a  separation  effect  created  as  the 
flow  unsuccessfully  attempts  to  conform  to  the  hullform.  When  M-l,  M-2  and  M-3  are 
compared  at  the  same  speed,  it  becomes  obvious  that  only  a  very  little  of  the  splashing 
seen  above  the  surface  in  M-l  can  be  attributed  to  spray  generated  by  the  struts.  The 
majority  of  the  splashing  effect  is  more  accurately  accounted  for  as  a  by-product  of 
eddy  resistance. 

WAKE  CHARACTERISTICS 

The  wake  produced  by  the  SLICE  hullform  is  compact  and  has  a  greater 
amplitude  than  that  of  the  SWATH.  The  SWATH  produces  a  relatively  flat  but 
significantly  longer  wake.  From  an  energy  in  the  water  standpoint,  comparison 
becomes  difficult  since  the  entire  volume  of  the  wake  must  be  considered. 

3.5  Test  Data 

Raw  data  collected  during  the  testing  is  presented  in  Appendix  F.  Test  1 
results  (fixed  in  pitch  and  heave  measurements  of  drag  and  lift)  are  presented 
graphically  in.  Drag  (positive  x)  is  taken  to  be  opposite  the  direction  of  forward 
motion.  Sinkage  (positive  z)  is  down  from  the  free  surface.  Several  speed  settings 
were  repeated  during  the  testing.  Repeated  speeds  are  recorded  in  these  figures  and  are 
designated  Dual  Point  (DP).  No  direct  comparison  of  performance  using  Figures 
(3-4)  -  (3-6)should  be  done  since  each  model  is  built  to  a  slightly  different  scale 
factor  . 

Figures  (3-7)  and  (3-8)  present  model  scale  data  at  the  same  scale  factor.  This 
required  scaling  M-2  and  M-3  data  to  the  M-l  scale.  Data  scaling  was  accomplished 
using  the  process  outlined  in  Chapter  2-7.  Scaling  between  models  was  accomplished 
using  the  composite  scale  factors  shown: 
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(3-1) 


^ mode!2  _  16.68 
^ modeli  ^7  ■  26 


=0.966 


Knoden  _  15.10 

^ modeli  T7  •  2 6 


=0.875 


(3-2) 


The  calculations  used  in  this  process  are  in  Appendix  G. 

Figure  (3-9)  shows  the  data  collected  for  Model  M-l  tested  in  both  the 
Baseline  configuration  and  an  in-line  configuration  of  the  forward  and  aft  strut/pod 
groups.  The  results  of  testing  Model  M-l  at  varying  depths  was  inconclusive. 
Insufficient  data  was  collected  to  develop  a  trend. 
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Figure  3-4:  Model  M-l  Fixed  in  Pitch  and  Heave 
Unsealed  Data 
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Figure  3-5:  Model  M-2  Fixed  in  Pitch  and  Heave 
Unsealed  Data 
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Figure  3-8:  Models  M-l/M-2  and  M-3  Fixed  in  Pitch  and  Heave 
Data  Presented  in  17.26:1  Scale 


Several  comments  can  be  made  about  this  data. 


REPEATABILITY  OF  RESULTS 

A  large  number  of  data  points  were  repeated  through  the  course  of  the  testing. 
Without  exception,  the  difference  in  recorded  data  is  statistically  insignificant.  Most 
repeated  data  was  recorded  in  successive  runs  down  the  tank.  A  notable  exception  to 
this  are  runs  70  and  71  in  Appendix  F.  In  this  case,  testing  of  Model  M-l  was 
completed.  Model  M-2  and  Model  M-3  were  tested,  then  Model  M-l  was  remounted 
and  realigned  on  the  carriage  and  five  data  points  collected.  There  was  excellent 
correlation  between  these  points  and  data  collected  two  days  earlier. 

HUMP  SPEED  MODEL  PERFORMANCE 

The  resistance  values  for  M-l  and  M-2  through  the  hump,  Figure  3-7,  support 
the  observation  made  earlier  that  the  wakes  are  very  similar.  The  similarity  in 
resistance  values  was  an  unanticipated  result.  M-2  represents  the  extent  to  which  L/D 
can  be  extended  and  still  retain  the  features  of  SLICE.  Since  this  variation  was 
insufficient  to  impact  the  amplitude  of  the  hump,  the  only  other  alternative  for 
diminishing  the  hump  magnitude  is  altering  the  relative  position  of  struts  and  pods. 

The  fact  that  the  hump  for  M-l,  M-2  and  M-3  all  occur  at  essentially  the  same 
speed  indicates  that  the  Strut  Length  is  the  appropriate  length  to  characterize  SLICE. 
This  conclusion  is  reached  by  considering  M-3  in  contrast  to  M-l  and  M-2.  The  hump 
is  expected  to  occur  at  &=  0.45.  For  the  demihull  of  M-3  this  corresponds  to 
V=  6.25  fps.  for  the  struts  it  corresponds  to  V=  3.2  fps.  The  hump  actually  occurs  at 
V=4.5, 

demonstrating  dominance  by  the  struts.  The  pods  of  M-l  and  M-2  are  less  than  half 
the  length  of  M-3's  demihulls  and  yet  the  hump  for  all  three  is  at  the  same  speed.  The 
struts  for  all  three  models  were  essentially  the  same  length.  The  underwater  bodies 
may  be  contributing  to  wave  making  through  cancellation  effects  but  the  position  of 
the  hump  is  driven  by  the  surface  piercing  bodies. 
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HIGH  SPEED  MODEL  PERFORMANCE 

In  the  range  6.5-10  fps  (full  scale  equivalent:  16-25  knots),  M-l  total  resistance 
is  up  to  20%  greater  than  that  of  M-2  even  though  M-l  has  less  wetted  surface  and 
therefore  less  frictional  resistance  than  M-2.  This  is  the  first  indication  that  model 
form  must  be  traded  off  with  friction  to  optimize  post  hump  resistance. 

For  speeds  greater  than  10  fps,  the  rapid  increase  in  M-l  resistance  can  be 
directly  correlated  with  the  appearance  of  the  trailing  edge  vortex  noted  in  the  physical 
observations. 

M-l/  LMSC  COMPARISON 

une  test  objective  for  M-l  was  to  compare  M-l  results  with  the  December 
LMSC  model.  Figure  3-10  presents  the  M-l  data  as  collected  as  well  as  a  sample  of 
LMSC  data  from  both  August  and  December  scaled  to  17.26:1.  Model  M-2  results  are 
also  presented  scaled  accordingly.  The  August  LMSC  curve  represents  LMSC  data 
runs  106-123,  The  December  LMSC  curve  represents  LMSC  data  runs  86-97.  These 
models  were  tested  with  appendages  and  haunch  structures  in  place. 

M-l  correlation  with  the  December  tests  is  good  below  8.25  fps  (full  scale:  20 
knots),  becomes  marginal  from  8.25-10  fps  (full  scale:20-25  knots),  and  does  not  exist 
above  this  value.  M-l  correlation  is  more  satisfactory  with  the  August  test  results, 
i  ne  lack  of  high  speed  correlation  requires  comment. 

Of  concern  is  the  vortex  created  at  10  fps  (full  scale:  25  knots).  The 
phenomenon  is  documented  as  existing  in  the  August  and  December  tests  as  well  as 
the  tests  conducteo  on  M-l .  The  question  is  deciding  if  the  lack  of  correlation 
indicates  that  effect  does  not  scale,  or  if  there  is  some  other  mechanism  that  is  causing 
the  difference.  The  possibility  that  the  LMSC  data  may  have  been  obtained  before 
steady  state  was  achieved  has  already  been  discussed.  Although  this  is  a  fault  of  the 
LMSC  test  program,  it  is  difficult  to  attribute  the  entire  difference  in  test  results  on 
this  shortcoming  alone. 
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Drag  Force  (lbf) 


0  5  10  15  20 

Velocity  in  FPS 


Figure  3-10:  Models  M-l  Performance  Compared  With  LMSC  Results 
Data  Presented  in  17.26:1  Scale 


In  the  tests  of  M-l.  there  is  direct  impingement  of  the  rooster  tail  on  the  model 
frame.  During  the  testing,  this  was  considered  acceptable  for  two  reasons:  1) 
Impingement  also  occurred  in  the  LMSC  model  tests  and  will  occur  in  the  full  scale 
ship  depending  on  the  construction  of  the  box  structure.  2)  During  the  testing,  the 
impingement  was  considered  a  minor  effect  since  the  profile  of  the  frame  structure 
was  small  (60  in:  total).  Reviewing  video  footage  of  these  tests,  the  direction  of 
rooster  tail  flow  was  perpendicular  to  the  frame,  thereby  maximizing  the  force  of 
impingement.  In  the  LMSC  tests  impingement  was  not  as  direct.  Appendix  H 
provides  a  calculation  that  demonstrates  that  impingement  could  account  for  up  to 
1 1  pounds  of  the  measured  force  at  15  fps. 

There  are  enough  differences  in  test  configurations  to  preclude  a  definitive 
explanation  for  the  lack  of  correlation.  The  rooster  tail  impingement  may  contribute  to 
the  lack  of  correlation.  The  conclusion  drawn  at  this  point  is  that  the  trailing  vortex  is 
real  and  will  scale  but  may  not  directly  contribute  to  overall  resistance  to  the  extent 
indicated  by  the  Model  M-l  results  in  Figure  (3-7).  Failure  to  achieve  steady  state  in 
the  LMSC  test  may  also  be  a  significant  factor  in  this  difference.  Re-testing  M-l  with 
the  frame  isolated  will  help  eliminate  this  ambiguity  but  ultimately  full  scale  test 
results  are  required. 

3.6  Full  Scale  Predictions  and  Discussion 

With  the  understanding  that  the  M-l  results  may  over  predict  full  scale 
resistance,  figures  (3-11)  and  (3-12)  present  full  scale  data  for  all  three  models.  Figure 
(3-13)  provides  power  requirements  based  on  these  values.  Data  scaling  was 
accomplished  using  the  process  outlined  in  Chapter  2-7.  The  calculations  used  in  this 
process  are  at  Appendix  G.  Based  on  figure  (3-11),  the  high  speed  SWATH  offers  the 
best  resistance  performance  through  the  entire  speed  range. 

Figures  (3-14)-  (3-19)  present  M-l,  M-2  and  M-3  data  in  model  and  full  scale, 
decomposed  into  the  frictional,  residuary  and  correlation  factor  components.  In  the 
post  hump  region  of  the  resistance  curve,  residuary  (ie  wave  making)  resistance 
accounts  for  at  least  60%  of  the  total  resistance  value.  This  further  supports  the 
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concept  that  the  hull  form  must  be  designed  as  a  compromise  between  frictional  and 
wave  making  resistance,  not  a  single  component  optimization. 
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Drag  force  (Ibt) 


EHP  (horsepower) 


Velocity  (knots) 


Figure  3-13:  Models  M-l/M-2  and  M-3  Fixed  in  Pitch  and  Heave 

Powering  Data 
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Figure  3-16:  Model  M-2  Resistance  Components  as  a  Percent  of  Total  Resistance 

Derived  from  Unsealed  Test  Data 


Percent  (as  fraction) 


Percent  (as  fraction) 


Velocity  ( knots) 

Figure  3-19:  Model  M-3  Resistance  Components  as  a  Percent  of  Total  Resistance 

Based  on  Full  Scale  Results 


Chapter  4  Numerical  Solution 


4.1  Objective  of  a  Numerical  Approach 

The  ability  to  predict  a  ship's  resistance  using  numerical  methods  has  many 
advantages.  If  an  appropriate  method  can  be  found,  the  expense  and  time  of  model 
testing  may  be  reduced.  Current  state  of  the  science  is  approaching  realistic  treatment 
of  the  free  surface  (ie:  non  linearized  theory)  but  these  methods  are  still  too  expensive 
to  use  as  an  iterative  design  tool.  Simpler  less  expensive  methods  may  not  be  accurate 
enough  to  exactly  predict  resistance  but  can  be  extremely  useful  for  comparative 
analysis  of  small  changes  to  the  same  hullform.  These  methods  require  tow  tank 
model  testing  as  verification  to  ensure  that  all  physical  phenomena  are  properly 
represented. 

This  thesis  attempts  to  validate  the  SLICE  applicability  of  a  numerical  method 
solution  of  Mitchell's  integral  as  an  approximation  for  wave  resistance.  As  cited  in 
Chapter  2.  the  approach,  first  outlined  by  Lunde.  was  used  by  Wilson  in  the  Wave 
Cancellation  Multihull  Ship  concept.  It  is  not  expected  that  this  thin  ship 
approximation  will  provide  exact  resistance  values.  The  goal  of  this  chapter  is  to 
develop  an  understanding  of  the  body  interactions  in  the  SLICE  hullform  to  assist 
resistance  optimization  through  body  placement.  Of  specific  interest  is  reducing  the 
magnitude  of  the  wave  making  hump.  The  thin  ship  approach  is  useful  since  it  allows 
wave  making  resistance  to  be  segregated  into  all  component  body  and  interaction 
terms.  This  is  important  while  studying  the  effect  of  small  geometry  changes  on  total 
resistance. 

Model  M-2  is  selected  as  the  subject  for  this  stud}-.  Physical  observation  during 
the  model  tests  demonstrate  that  this  model  has  a  reasonable  chance  of  satisfying  thin 
ship  flow  requirements. 
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4.2  Solution  Technique 

Restated  from  Chapter  2.8.  the  equation  to  be  solved  is: 


\  NN  NN 

Rw^=16it  PK2] (LtPiPJ+tTQ.QJ)seciQd6 
n  1  1  11 


(4-1) 


where: 


TL 

P=j jo(^^)<x^[K(xcosQ+ysinQ)soc2Q]eKsec2ddxdz  (4-la) 

0  0 


TL 

<2=j’J'o(x^)sin[K(!xccs6  +ysind)sec2d]<?  ’asec2e dxdz  (4-lb) 

oo 


2%  6x 


(4-lc) 


k=-£-  (4- Id) 

V2 

The  coordinate  system  that  will  be  used  has  its  origin  in  the  undisturbed  free  surface 
with  +X  in  the  direction  of  ship  motion,  +Z  vertically  up  from  the  free  surface  and 
+Y  in  accordance  with  the  right  hand  rule. 

DEVELOPMENT  OF  THE  BODY  AND  POTENTIAL  FUNCTIONS 

The  first  requirement  is  to  write  a  function  that  describes  the  strut  and  pod 
forms.  This  function  is  written  with  the  bodies  at  the  global  origin.  Later,  these  bodies 
will  be  translated  to  their  actual  relative  positions  through  the  X,  Y  and  Z  terms  in 
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Equation  (4-la)  and  (4-lb). 

For  M-2,  the  strut  form  was  developed  using  an  identical  parabolic  nose  and 
tail  shape  with  parallel  midbody  in  between.  The  equations  used  to  describe  the  strut 
shapes  are: 

NoseOfjseis= &  * [  1  (4-2) 


MidbodyOffseis = ~ 
2 


(4-3) 


TailOffsets=^  *[l-( 
2 


(4-4) 


where  for  this  specific  strut 

B  (the  maximum  beam)  =  0.195  ft 
Ln  (nose  length)  =  Lt  (tail  length)  =  0.51  ft 
Lpmb  (parallel  midbody  length)  =  0.42  ft 
n  =  o  05 

‘‘strut 

To  reduce  the  complexity'  of  the  strut  body  potential  function  as  applied  in  (4-la)  and 
(4-lb),  the  offsets  calculated  by  (4-2),  (4-3)  and  (4-4)  were  curve  fit  to  a  single  high 
order  polynomial  given  by: 


Strut =0.4893  ** 6  +2. 1 126  5  +3. 1328  **4 + 1.7248  **3  -0.0786  **2  -0.3890  ** +0.0008 

(4-5) 
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The  source  distribution  representing  the  strut  is  written  as: 


°sm=  *(2.936  **5 + 10.563  *x4 + 12.53 1  +5.174  *x2 -0. 1572  *x- -0.389) 


The  pod  shapes  are  created  in  a  similar  fashion.  Now  however,  the  nose  is 
represented  by  an  ellipse  and  the  tail  by  a  parabola.  The  pod  offsets  are  given  by  the 
following: 


PodNose=—*[ l-(^-f)  *”']'’«* 

2  Ln 


(4-7) 


D 


PodMidbody=— 


(4-8) 


„.„./)  X-L-L,  n*a 

PodTail=—*[l-( - .  *"*)  ] 

2  L  *  J 


(4-9) 


where  for  this  specific  pod. 

D  (the  maximum  diameter)  =  0.4635  ft 

Ln  (nose  length)  =  0.67  ft 

L,  (tail  length)  =  1.1  ft 

Lpmb  (parallel  midbodv  length)  =  0.60  ft 

IW  -  2.00 

n*i  =  2.40 

Again  a  single  high  order  polynomial  was  generated  to  represent  the  pod  and  is  given 
by: 
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Pod=-QM*x1-Q32*xt-l.l6*x5-2'21*xA-2.96*x*-221*x2-Q.99*x+Q.m 


The  source  distribution  representing  the  pod  is  written  as: 

apod=~~~  *(  -0.253  *x6  -1.93  *x  5  -5.82  *xA  -9.49  *x2  -8.88  *x2  -4.54  **-0.99 

271 

The  pods  are  treated  as  a  line  source  in  (4- la)  and  (4- lb)  by  fixing  the  z  term 
in  the  exponential  as  a  constant  equal  to  the  strut  depth  plus  the  maximum  radius  of 
the  pod.  This  follows  the  example  provided  by  Chapman19  for  dealing  with  submerged 
bodies  of  revolution. 

PAIRING  OF  TERMS 

Struts  and  pods  are  numbered  in  the  traditional  naval  manner.  The  scheme  is 
presented  in  Figure  (4-1).  There  are  a  total  of  36  possible  combinations  of  the  eight 
bodies  that  form  the  SLICE  hull.  Eight  combinations  are  "like"  or  body  terms  and  28 
combinations  are  different  or  "cross  product"  terms.  These  terms  are  presented  below 
grouped  by  relative  position  in  the  hull.  The  numbers  indicate  which  bodies  are 
interacting  in  the  term  calculated: 

1)  Body  Terms: 

1-1,  2-2,  3-3,  4-4,  5-5,  6-6,  7-7.  8-8 

2)  Same  Comer  Cross  Terms: 


19  Chapman.  R.B..  Hydrodynamic  Drag  of  Semisubmerged  Ships,  American  Society  of  Mechanical  Engineers, 
paper  number  72-Wa,  1972.  pp  879-884. 
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CODING 


MATLAB  was  used  for  the  required  computations.  This  choice  was  made  to 
utilize  the  installed  graphics  capabilities  of  the  program.  Plotting  results  as 
computations  were  made  was  useful  both  in  validating  geometry  assumptions  and 
debugging.  A  logic  flow  chart  of  the  overall  computation  scheme  is  at  Appendix  I. 

The  routines  used  are  also  presented  in  this  appendix.  The  pod  and  strut  geometries 
were  developed  in  the  files  PODGEO.m  and  STRUTGEO.m,  The  main  wave  making 
resistance  routine  is  contained  in  the  file  Q.m.  The  solution  to  equations  (4-1),  (4-la) 
and  (4- lb)  are  solved  in  called  subroutines  for  each  of  the  thirty  six  paired  terms  in 
files  named  for  the  bodies  that  are  paired.  Samples  of  a  strut  body,  pod  body’  and 
strut/pod  cross  term  are  contained  in  the  files  ONEONE.iru  SIXSIX.m  and 
ONESIX.m.  All  other  body  and  cross  product  subroutines  are  patterned  after  these 
files.  This  apparently  cumbersome  organization  had  one  significant  advantage.  Once 
the  files  were  written  and  assembles,  bodies  could  be  toggled  on  and  off  in  the 
computation  simply  by  adding  a  at  the  head  of  the  appropriate  subroutine  call 
line.  This  allows  the  ov  erall  SLICE  geometry  to  be  reduced  into  smaller  pieces  rather 
easily  and  helped  demonstrate  the  complex  interactions  that  occur  in  the  full 
multibodied  hull. 

4.3  Results 

Only  partial  results  were  obtained  for  the  numerical  solution  during  the 
allocated  period  of  study  for  this  thesis.  Complete  validation  of  results  was  not 
obtained  and  the  results  contained  in  this  section  require  further  analysis  before  design 
decisions  can  be  made.  Sufficient  progress  was  made  to  draw  some  preliminary 
conclusions.  Model  M-2  was  computed  in  the  speed  range  of  3-9  fps  (full  scale  11-22 
knots)  using  the  values  shown  in  Table  (4-1).  The  results  presented  in  Figure  (4-2) 
and  Table  (4-2)  are  based  on  using  fifty  subdivisions  in  the  x  integration,  12  divisions 
in  the  z  integration  and  10  divisions  in  the  0  integration.  Testing  was  conducted 
varying  the  number  of  divisions  and  this  was  found  to  be  a  minimum  acceptable  limit 
with  the  integrity  of  the  result  maintained. 
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Strut  Length 

1 .44  ft 

Strut  Beam  (max) 

0.195  ft 

Strut  Height 

0.36  ft 

Fwd  Strut  Set  Back  (x  axis  reference) 

0.0  ft 

Fwd  Strut  Offset  (y  axis  reference) 

±  0.979  ft 

Aft  Strut  Set  Back  (x  axis  reference) 

-3.48  ft 

Aft  Strut  Offset  (y  axis  reference) 

±  1.396  ft 

Strut  Submergence  (z  axis  reference) 

-0.18  ft 

Pod  Length 

2.321  ft 

Pod  Diameter  (max) 

0.4635  ft 

Fwd  Pod  Setback  (x  axis  reference) 

0.0  ft 

Fwd  Pod  Offset  (y  axis  reference) 

±  0.979  ft 

Aft  Pod  Set  Back  (x  axis  reference) 

-3.16  ft 

Aft  Pod  Offset  (y  axis  reference) 

±  1.396  ft 

Pod  Submergence  (z  axis  reference) 

-0.59  ft 

Table  4-1:  M-2  Gross  Parameters  Used  For  Calculation 

Figure  4-2  presents  only  the  wavemaking  portion  of  resistance  as  predicted  by 
the  numerical  scheme.  The  negative  overall  result  is  attributed  to  the  overall  sign 
convention  that  is  used  in  the  calculation  which  has  x  positive  in  the  direction  of 
motion.  An  unresolved  problem  with  this  conclusion  is  that  the  body  terms  (1-1 
through  8-8)  should  have  been  negative  on  their  own.  This  however  was  not  the  case. 
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R  (Ibf) 


Velocity  (fps) 

Figure  4-2:  Wave  Making  Resistance  of  SLICE  calculated  from  Thin  Ship  Theory’ 
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Speed 

(fps) 


1-1 

(Ibf) 


0.39E-3 


0.56E-3 


0.62E-3 


0.65E-3 


0.65E-3 


0.64E-3 


2-2 

(Ibf) 


0.39E-3 


0.56E-3 


0.62E-3 


0.65E-3 


0.65E-3 


0.64E-3 


0.63E-3  0.63E-3  0.63E-3 


3-3 

(Ibf) 


0.39E-3 


0.56E-3 


0.62E-3 


0.65E-3 


0.65E-3 


0.64E-3 


4-4 

(Ibf) 


0.39E-3 


0.56E-3 


0.62E-3 


0.65E-3 


0.65E-3 


0.64E-3 


0.63E-3 


5-5 

(Ibf) 


0.0004 


0.0020 


0.0038 


0.005- 


0.0061 


0.0067 


0.0070 


6-6 

(Ibf) 


0.0004 


0.0020 


0.0038 


0.005: 


0.0061 


0.0067 


0.0070 


Table  4-2a:  Wave  Making  Resistance  Component 

BODY  TERMS 


7-7 

(Ibf) 


0.0004 


0.0020 


0.0038 


0.0052 


0.0061 


0.0067 


0.0070 


erms 


8-8 

(Ibf) 


0.0004 


0.0020 


0.0038 


0.0052 


0.0061 


0.0067 


0.0070 


Speed 

(fps) 

2-6 

(Ibf) 

4-8 

(Ibf) 

3 

-0.556 

-0.556 

-0.181 

-0.181 

4 

-1.413 

-1.413 

-1.037 

-1.037 

5 

-2.043 

-2.043 

-1.773 

-1.773 

6 

-2.406 

-2.406 

-2.220 

-2.220 

7 

-2.598 

-2.598 

-2.466 

-2.466 

8 

-2.700 

-2.700 

mm 

9 

-2.751 

-2.675 

-2.675 

Tab 

le  4-2b:  Wave  Making  Resistance  Component  Terms 

SAME  CORNER  CROSS  TERMS 
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(lbf)  (lbf)  (lbf)  (lbf)  (lbf)  (lbf) 


0.276  -0.249  -0.245  0.001  0.421  -0.241  -0.436  0.0 


0.11  -0.547  -0.103  0.002  0.141  -1.054  -0.488  0.004 


0.356  j  -0.135  -0.528  0.0  0.350  0.383  -0.371  -0.003 


j  0.11  ! 

-0.547 

0.356 

-0.135 

-0.264 

1.128 

-0.412 

1.325 

-0.388 

1.142 

-0.103  ! 

1 

1 

0.682 

i 

BE 


M 


07  -0.422  2.035  1.709  -0.010 


2.059  -0.008 


1.106  1.527 


8  -0.882 


006  -0.903 


m 


0.302  0.799 


Table  4-2c:  Wave  Making  Resistance  Component  Terms 
LONGITUDINAL  CROSS  TERMS 


2-3 

2-7 

(lbf) 

(lbf) 

0.02E-3  i 

! 

-0.043 

-.1  IE-3 

-0.045 

-.09E-3 

-0.016 

-.30E-3 

0.870 

-.09E-3 

0.540 

-.31E-3 

0.986 

0.01E-3 

0.072 

-0.004 


0.0 


Table  4-2d:  Wave  Making  Resistance  Component  Terms 
DIAGONAL  CROSS  TERMS 
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Understanding  that  the  code  used  is  only  partially  validated,  the  following 
observations  are  made: 

STRUTS 

The  contribution  made  by  the  struts  is  small.  This  is  true  both  of  the  body 
terms  and  the  cross  terms  between  any  two  struts.  This  is  reasonable  for  the  strut 
geometry  used.  There  is  symmetry  about  midships,  the  nose  and  tail  sections  are 
slender  and  the  overall  L/D  is  7.4.  To  verify  that  the  low  values  obtained  were  not  a 
code  problem,  an  alternate  geometry  was  selected  that  did  not  possess  midships 
symmetry  but  retained  the  same  L/D.  A  significant  increase  in  resistance  values  was 
noted.  These  results  follow  the  conclusions  drawn  by  Chapman  and  Zhu  Bing-quan  on 
strut  optimization  cited  in  Chapter  2.5. 

CROSS  TERMS 

As  anticipated,  the  Same  Comer  Cross  terms  have  the  largest  impact  on  wave 
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making  resistance.  This  is  followed  in  importance  by  the  Longitudinal  Cross  Terms. 
The  Diagonal  and  Transverse  Cross  Terms  had  more  significant  impact  than  was 
anticipated.  Although  smaller  than  the  other  terms,  it  would  not  be  correct  to  eliminate 
the  Transverse  terms  from  optimizing  iterations.  Also  clear  from  a  review  of  Table 
(4-2)  is  the  complexity  of  the  body  interactions.  At  a  given  speed,  the  transverse  term 
for  a  specific  body  may  have  an  additive  effect  while  the  longitudinal  term  for  that 
same  body  has  a  subtracting  term  etc.  Additionally,  due  to  the  lack  of  fore  aft 
symmetry  in  the  baseline  SLICE  layout,  similar  body  groupings  do  not  behave  the 
same  at  a  given  speed.  Note  the  difference  between  3-8  and  7-4  in  the  Transverse 
Terms  Table  as  an  example  of  this.  This  decomposition  of  SLICE  wavemaking 
resistance  demonstrates  that  optimization  work  must  be  conducted  with  the  entire 
hullform  present.  Attempting  to  sub-optimize  a  component  body  (such  as  a  single  strut 
pod  combination)  without  observing  the  impact  on  the  entire  hullform  may  lead  to 
erroneous  conclusions. 
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Chapter  5  Conclusions 


5.1  Restatement  of  Thesis  Objectives 

This  thesis  attempted  to  add  to  the  developing  body  of  SLICE  knowledge  by 
focusing  on  the  resistance  and  powering  aspects  of  the  SLICE  vessel  in  an 
academically  rigorous  manner.  Three  specific  objectives  were  identified: 

1)  Independent  research  into  the  SLICE  hullform  to  give  both  the  US  Navy 
and  LMSC  a  second  hydrodynamic  viewpoint  of  the  177  Lton  vessel  prior  to  the 
ATT). 

2)  Development  of  a  larger  body  of  data  to  aid  in  optimizing  the  177  Lton 
hull  form. 

3)  A  broader  understanding  of  the  hydrodynamic  features  of  the  hull  to 
facilitate  deciding  the  merits  of  pursuing  designs  in  other  displacements  more  practical 
to  the  US  Navy. 

In  order  to  meet  these  objectives,  several  questions  required  answers: 

1)  LMSC  theorizes  that  the  appropriate  characteristic  length  for  SLICE  is  the 
length  of  a  single  pod.  Is  this  valid? 

2)  The  specific  interaction  between  bodies  has  not  been  defined  by  LMSC. 
What  is  the  relative  importance  of  form;  strut/pod  alignment;  and  the  effects  that  a 
strut/pod  combination  have  on  any  of  the  other  three  strut/pod  combinations  that  form 
the  underwater  body?  Understanding  these  effects  is  necessary'  for  design  optimization. 

3)  The  LMSC  proposed  design  places  the  forward  struts  inboard  of  the  aft 
struts.  Is  the  outward  lateral  force  (y  axis  lift)  created  by  an  effective  angle  of  attack 
on  the  aft  struts  significant  enough  to  create  structural  concerns? 

4)  The  LMSC  model  tests  were  all  conducted  with  models  fixed  in  pitch  and 
heave.  Is  this  an  accurate  enough  representation  of  reality  to  make  predictions? 

5)  The  LMSC  tests  maintained  speed  increments  for  an  extremely  short 
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duration  (usually  less  than  10  seconds)  and  frequently  truncated  test  runs  at  a  model 
scale  equivalent  of  25  knots  full  scale.  Was  this  sufficient  to  characterize  a  ship  with 
an  intended  speed  of  30+  knots? 

6)  Is  the  apparent  focus  on  skin  friction  reduction,  as  reflected  in  the  body 
forms  selected,  reasonable? 

7)  Does  SLICE  offer  the  stated  hydrodynamic  advantage  over  a  SWATH  ship. 

5.2  Conclusions 

CHARACTERISTIC  LENGTH 

The  appropriate  characteristic  length  for  SLICE  is  the  strut  length,  not  the  pod 
length.  This  is  illustrated  in  the  resistance  curves  for  M-l.M-2  and  M-3.  These  models 
had  significantly  different  underwater  bodies  but  very  similar  strut  lengths.  On  a 
resistance  vs  speed  plot,  the  position  of  the  wave  making  hump  for  all  three  hullforms 
is  nearly  identical.  In  order  to  align  these  with  S~  4.5.  strut  length  must  be  chosen  as 
the  characteristic  length  used  in  the  definition  S. 

THE  IMPORTANCE  OF  BODY  FORM  AND  BODY  INTERACTION 

Model  M-2  w  as  developed  with  the  expectation  that  finer  lines  would  help 
improve  the  hump  speed  performance  of  the  baseline  SLICE  hullform.  The  offsets 
selected  for  the  pods  increased  the  constructive  L/D  of  the  pods  from  4.1  in  the 
baseline  to  5.0  in  M-2.  This  was  the  maximum  improvement  that  could  be  made  to 
L/D  and  still  maintain  the  gross  characteristics  of  the  Baseline  SLICE.  This  change  in 
L/D  was  not  sufficient  to  impact  the  hump  speed  performance.  If  the  SLICE  hullform 
is  maintained,  improvements  in  hump  speed  resistance  must  come  from  body 
interaction.  Underwater  photography  showed  that  propeller  submergence  is  better 
through  the  hump  speed  than  previously  thought. 

The  finer  lines  of  M-2  did  make  a  significant  improvement  in  the  resistance 
measured  in  the  high  speed  range.  Underwater  photography  revealed  that  the  bluff 
lines  of  M-l  causes  the  shedding  of  a  significant  vortex  at  the  trailing  edge  of  the 
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struts  starting  at  model  scale  equivalents  of  25  knots.  These  vortices  were  incorrectly 
interpreted  as  spray  phenomenon  in  previous  testing.  In  the  full  scale  ATD  ship,  the 
flow  modification  caused  by  the  propellers  may  cause  these  vortices  to  be  ingested 
into  the  propellers  with  a  significant  loss  of  thrust. 

MODEL  TEST  CONDITIONS 

Attempting  to  test  the  SLICE  hullform  in  a  condition  other  than  fixed  in  pitch 
and  heave  introduces  a  significant  level  of  complexity  to  the  model  design.  The  added 
complexity  does  not  necessarily  make  the  results  more  useful.  The  difficulty  lies  in 
developing  a  system  that  connects  the  model  to  the  tow  carriage  through  the  line  of 
force.  With  four  underwater  bodies  all  towed  from  a  central  point,  this  is  extremely 
complex  to  manufacture.  Any  connection  point  somewhere  other  than  in  the  line  of 
force  is  easier  to  manufacture  but  creates  variable  applied  moment.  This  effect  cannot 
simply  be  subtracted  out  since  it  is  impossible  to  separate  it  from  the  effects  of  real 
hydrodynamic  forces.  The  remaining  option  then  is  to  lock  the  model  in  pitch  and 
heave.  For  the  purposes  of  conducting  comparative  testing,  this  simplification  is 
justified  and  was  a  reasonable  simplification  made  by  LMSC. 

The  tests  conducted  as  part  of  this  thesis  showed  that  there  was  a  significant 
transient  time  associated  with  acceleration  from  rest.  In  the  higher  speed  ranges,  forces 
required  4-6  seconds  to  achieve  steady  state  after  the  steady  state  speed  was  achieved. 
The  LMSC  tests  incorporated  more  speeds  in  a  single  run  with  fewer  accelerations 
from  rest.  LMSC  allowed  for  significantly  less  data  collection  time  at  each  speed 
increment  (less  than  10  seconds  compared  with  greater  than  20  seconds).  This  may  not 
have  been  sufficient  time  to  achieve  steady  state  and  may  contribute  to  the  differences 
in  LMSC's  test  results  and  those  conducted  as  part  of  this  thesis. 

The  numerical  analysis  conducted  in  Chapter  4  highlighted  the  extreme 
complexity  of  the  flow  around  SLICE.  Relationships  presumed  negligible  such  as  the 
transverse  and  diagonal  body  interactions  were  shown  to  contribute  noticibly  to  the 
overall  resistance  value.  This  has  an  implication  in  the  LMSC  test  results.  When  the 
August  control  surface  tests  were  conducted,  only  one  side  of  the  model  (starboard) 
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was  adjusted  for  various  angles  of  attack.  Based  on  the  results  of  the  wave  making 
decomposition  however,  drawing  conclusions  about  whole  body  behavior  based  on 
partial  bod}'  measurements  may  lead  to  incorrect  results. 

HULLFORM  OPTIMIZATION 

The  data  decomposition  presented  at  the  end  of  Chapter  3  demonstrates  that  it 
is  incorrect  to  optimize  the  hull  for  high  speed  operation  by  focusing  on  frictional 
resistance.  Powering  through  the  wave  making  hump  does  not  imply  that  wave  making 
resistance  becomes  negligible.  It  is  true  that  the  coefficient  of  wave  making  resistance 
Cw  becomes  a  small  number  but  resistance  continues  to  increase  as  velocity  squared. 
The  design  perspective  should  be  that  Froude  scalable  effects  will  account  for  at  least 
half  the  total  resistance  and  the  hull  should  be  optimized  accordingly. 

SLICE  APPLICABILITY  TO  THE  US  NAVY 

Figure  (5-1),  obtained  from  a  set  of  presentation  graphics  used  by  Dr.  R. 
Compton  at  the  U.S.  Naval  Academy,  depicts  a  predecessor  to  SLICE.  The  literature 
survey  conducted  as  part  of  this  thesis  did  not  uncover  any  design  work  accomplished 
on  this  hull. 

The  comparison  of  either  M-l  or  M-2  with  M-3  indicate  that  if  there  is  an 
advantage  to  the  SLICE  hullform  over  conventional  SWATH  huilforms,  it  must  be 
found  somewhere  other  than  in  powering  performance.  M-3  demonstrates  that  it  is 
possible  to  create  a  SWATH  to  the  design  constraints  of  SLICE  and  still  meet  or 
exceeded  SLICE  resistance  performance  through  the  entire  speed  range. 

SLICE  offers  a  possible  advantage  in  the  ratio  of  Deck  Area  to  Displacement. 
The  development  of  M-3  showed  that  the  177  Lton  SWATH  could  have  up  to  600  ft2 
less  deck  area  than  the  equivalent  SLICE.  The  extent  to  which  the  crossbox  is 
cantilevered  fore  and  aft  of  the  struts  becomes  the  controlling  factor.  For  low  density 
payloads,  such  as  ferry  traffic,  this  is  significant.  For  higher  density  payloads  typical 
of  military  craft  the  advantage  is  diminished. 
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The  following  recommendations  are  made  for  future  study  of  SLICE 
technology: 

1)  The  question  of  lateral  lift  on  the  aft  struts  created  by  the  flow  direction 
remains  unanswered.  This  question  should  be  addressed  before  the  ATD  vessel  is 
placed  in  service  to  ensure  a  structural  issue,  either  from  an  integrity  or  fatigue 
standpoint,  is  not  overlooked. 

2)  This  thesis  presented  a  reasoned  approach  to  scaling  model  data  to  full 
scale.  The  ATD  will  be  the  validating  test  of  assumptions  and  decisions  made  in 
presenting  the  case  for  using  traditional  methods  of  data  scaling.  If  conducted 
properly,  the  ATD  will  be  invaluable  in  demonstrating  the  robustness  of  Froude 
hypothesis.  During  the  ATD,  significant  effort  should  be  invested  in  photographically 
documenting  flow  behavior  both  above  and  below  the  surface.  Without  this 
observation  of  physical  phenomena,  specific  conclusions  about  scaling  issues  will  not 
be  possible. 

3)  The  numerical  approach  used  in  this  thesis  requires  further  development 
before  tuning  of  the  SLICE  hullform  can  be  fully  demonstrated.  Thin  ship  theory 
appears  to  be  reasonable  for  comparative  analysis  purposes  and  this  approach  should 
be  further  developed.  After  further  validation,  SLICE  component  positioning  can  be 
optimized  and  then  a  model  of  the  final  outcome  should  be  built  and  tested  in  a 
similar  scale  to  M-2  so  that  results  can  be  compared.  This  could  be  the  subject  of  a 
reasonable  follow  on  thesis. 

4)  The  SLICE  ATD  has  the  potential  to  answer  many  questions  that  have  been 
raised  regarding  the  classical  naval  architecture  approach  to  resistance  predictions. 

With  careful  planning,  the  ATD  will  be  much  more  than  a  demonstration  of  a  single 
ship's  performance.  It  has  the  potential  to  contribute  to  the  Naval  Architecture 
community's  understanding  of  the  robustness  of  current  testing  and  scaling  techniques. 
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Appendix  A 


Model  Testing  Lessons  Learned 
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The  Science  and  Apt  of  Model  Testing 


The  purpose  of  this  appendix  is  to  outline  the  thought  process  that  was  used  in 
designing  the  experiments  and  the  models  that  were  used  in  this  thesis.  The  tone  of 
this  section  will  be  conversational  and  it  represents  a  compilation  of  lessons  learned  as 
the  research  for  this  paper  progressed.  It  is  hoped  that  this  appendix  may  be  of  use  to 
other  students  who  may  choose  to  undertake  a  similar  type  of  research  effort. 

When  the  work  for  this  paper  began,  it  was  decided  that  in  support  of  a  Naval 
Construction  program,  there  would  be  a  wealth  of  experience  to  be  gained  from 
designing  and  building  a  ship  (albeit  a  small  one)  that  actually  had  to  work.  With  this 
decided,  I  embarked  on  a  challenging,  frustrating,  rewarding,  extremely  time 
consuming,  detail  intensive  process. 

The  title  to  this  appendix  is  not  meant  to  be  trite.  There  are  clearly 
scientific/engineering  principles  that  may  not  be  violated  if  reasonable  test  results  are 
expected.  Geometric  similitude  must  be  maintained,  tolerances  must  be  decided  on, 
load  and  material  strength  estimates  must  be  reasonably  predicted,  scaling  must  be 
calculated  properly  and  data  reduction  must  be  interpreted  accurately.  Beyond  this 
however,  I  discovered  an  entire  world  of  additional  issues  regarding  the  planning  of 
the  experiment  actual  construction  of  the  models  and  successful  time  and  resource 
management  during  the  critical  and  always  too  short  time  allocated  in  the  test  facility. 

\  ou  cannot  gain  an  appreciation  for  this  facet  of  experimentation  by  taking  a  course 
or  reading  a  book. 

RECOMMENDATION  #1:  Prior  to  doing  anything  on  your  own  experiment, 
observe  similar  projects  for  several  days.  In  the  case  of  this  research,  I  had  the 
opportunity  to  observe  the  SLICE  testing  conducted  at  NSWC  in  August  and 
December.  For  eight  days  that  I  was  able  to  watch  other  people  devise  solutions  and 
cope  with  frustrations  without  impacting  my  own  work.  Many  ideas  for  my  own 
project  were  bom  during  these  eight  days  and  most  importantly,  I  met  the  experts  in 
the  field  of  model  testing.  Later,  these  associations  were  invaluable  to  me. 
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If  this  period  of  observation  does  not  scare  you  away  from  proceeding  with 
}our  own  experiment  the  next  question  you  must  honestly  answer  is:  How  handy  are 
you  and  what  resources  do  you  have  at  your  disposal?  When  I  began  building  the 
models.  I  considered  myself  reasonably  proficient  with  tools  but  I  had  never 
undertaken  a  project  of  this  magnitude  before. 


RECOMMENDATION  #2:  As  a  means  of  determining  your  handiness,  I 
suggest  you  consider  the  type  tolerances  you  will  have  to  build  to  if  the  model  is  to  be 
considered  a  reasonable  geosim.  For  argument  lets  say  that  in  full  scale,  the  greatest 
building  error  that  will  be  tolerated  is  6  inches  (pretty  big  error!).  If  the  scale  that  you 
are  building  to  is  16:1,  inaccuracies  in  your  model  of  greater  than  3/8  inch  exceed  the 
building  tolerance  and  you  have  introduced  a  level  of  uncertainty  as  to  the  actual 
geometric  similarity  of  your  model  to  the  real  ship.  If  your  personality  is  such  that  you 
fail  to  recognize  the  difference  between  3/8  and  1.  2  inch  when  measuring  things,  then 
model  building  may  not  be  the  thing  for  you.  Scaling  errors  to  full  scale  to  determine 
reasonableness  was  a  process  that  slowly  evolved  during  mv  own  project 


At  this  stage,  some  decisions  need  to  be  made  regarding  the  specific  test 
objectives.  To  a  large  measure,  this  will  determine  the  complexity  of  the  model  to  be 
built.  The  speeds  to  be  modeled  must  be  matched  with  tow  tank  capabilities.  For  this 
project  full  scale  speeds  of  40  knots  were  desired.  Based  on  the  maximum  carriage 
speed  at  the  MIT  tow  tank,  this  would  have  required  building  the  models  in  a  100:1 
scale,  a  size  so  small  that  no  useful  data  would  be  collected.  The  alternatives  were  to 
either  slow  the  maximum  speed  or  find  another  facility.  The  United  States  Naval 
Academy  was  literally  the  closest  facility  with  the  necessary'  capability.  In  addition  to 
these  concerns,  there  were  also  the  questions  regarding  what  the  model  should  be 
capable  of  doing.  Early  on  in  the  project  it  was  decided  that  the  models  would  have 
no  control  surfaces.  This  removed  a  significant  level  of  complexity  from  the  model 
design.  It  was  also  decided  that  the  models  should  be  capable  of  being  modified  for 
self  propulsion  at  some  point  in  the  future.  Additionally,  since  the  general  test  scheme 
called  for  a  test  free  in  heave,  the  models  had  to  be  built  with  sufficient  reserve 
buoyancy  to  support  themselves,  the  mounting  frame,  the  towing  carriage  heave  post 
and  force  blocks  and  an  estimated  reserve  for  propulsion  equipment.  This  simple 
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exercise  placed  many  constraints  on  the  model: 

1)  The  models  would  have  to  be  particularly  rugged  to  enhance  the  chance  of 
survival  during  the  400+  mile  trip  to  the  test  facility  in  the  back  of  a  van.  I  eliminated 
the  conventional  foam  construction  in  deference  to  this  requirement. 

2)  By  understanding  what  the  models  had  to  do,  I  could  make  a  weight  budget 
to  build  to.  It  turned  out  that  the  heave  post  and  force  block  weighed  25  pounds 
collectively.  Properly  scaled.  I  could  estimate  the  displacement  of  the  models  (between 
75,  83  and  115  lbf  for  models  M-l,  M-2  and  M-3  respectively).  I  then  budgeted  10% 
for  reserve  buoyancy  and  allocated  15  lbf  for  frame  weight.  The  remainder  could  be 
used  in  materials  to  actually  build  the  models.  This  weight  budget  further  drove  the 
building  plan.  The  material  selected  had  to  be  strong,  extremely  light  and  simple 
enough  for  an  amateur  to  work  with. 


RECOMMENDATION  #3:  Design  your  models  only  for  the  capability  you 
plan  to  test.  Build  only  what  you  need.  Clearly  establish  exactly  what  the  models  will 
be  used  for  before  planning  their  construction.  Eliminate  multiple  functions  where 
possible  since  multiple  functions  will  most  likely  increase  the  model's  level  of 
complexity.  In  the  case  of  this  project,  significant  effort  was  expended  to  meet  an 
established  weight  budget.  The  weight  budget  was  based  on  the  assumption  that  the 
model  would  be  tested  in  a  free  to  heave  condition  but  this  test  was  never  performed. 
For  the  fixed  in  heav  e  tests,  the  models  would  have  been  acceptable  even  if  they  were 
negatively  buoyant. 


At  this  stage  there  is  almost  enough  preliminary  groundwork  set  to  begin 
building.  With  any  luck  you  will  be  able  to  enlist  the  help  of  others  in  the  actual 
building  process.  The  models  built  for  this  project  represent  the  collective  work  of  9 
individuals  with  a  total  of  roughly  700  manhours  invested  between  December  and 
March.  Admittedly  a  portion  of  this  time  was  expended  inefficiently,  devising 
solutions  to  emergent  problems,  obtaining  supplies  etc.  but  these  are  the  facts  of  model 
building  and  this  sort  of  time  expenditure  should  be  planned  for. 

RECOMMENDATION  #4  If  you  are  not  literate  in  a  CAD  system,  develop 
that  proficiency  well  in  advance  of  the  model  building.  As  early  as  possible,  make  full 
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scale  working  drawings.  This  is  critical  if  other  people  are  helping  you.  It  is  the  only 
way  to  know  that  you  helpers  understand  what  to  build.  If  the  drawings  are  done  soon 
enough  in  advance,  they  may  help  point  out  errors  in  the  design.  I  did  not  develop  a 
good  set  of  working  drawings  early  enough  in  the  process  and  it  hindered  progress  on 
several  occasions.  Once  the  drawing  process  was  brought  under  control,  I  found  it 
handy  to  cut  up  the  rather  disposable  drawings  for  templates. 

Originally  PVC  was  the  material  of  choice  for  the  model  construction.  It  is 
easy  to  shape  and  bond  and  is  a  relatively  inexpensive  alternative.  Additionally,  water 
absorption  is  not  a  concern.  There  are  two  drawbacks  though.  PVC  is  relatively  dense. 
At  81  lb/in3,  it  can  be  made  to  float  if  rather  thin  walled  sections  are  used.  This  is 
useful  for  long  stretches  of  parallel  midbodv  but  becomes  a  significant  drawback  for 
the  shaped  pieces  which  must  be  turned  on  a  lathe.  Hollowing  these  sections  out  to  an 
appropriate  wall  thickness  was  not  practical.  Additionally,  the  time  to  turn  a  single 
nose  cone  (inside  and  outside)  was  fund  to  be  about  10  hours  after  considerable 
practice.  A  total  of  16  nose  and  tail  cones  were  required  for  two  models.  160  hours 
was  considered  excessive  to  finish  this  single  piece  of  the  project.  The  second 
drawback  of  PVC  is  also  one  of  its  advantages.  PVC  is  relatively  easy  to  cut  because 
it  is  soft.  .Although  there  was  tittle  danger  of  it  deforming  in  a  gross  sense,  it  is  very 
easy  to  gouge  and  scratch.  Once  damaged  in  this  fashion,  it  is  difficult  to  repair 
because  although  it  is  easy  to  sand  and  it  bonds  well  to  itself  with  the  proper  adhesive, 
other  bonding  and  filling  agents  do  not  adhere  well.  In  short  the  probability  of 
achieving  a  satisfactory  repair  once  damaged  is  small.  In  tight  of  the  trip  the  models 
would  take  to  the  test  facility  and  the  potential  for  damage  during  this  trip,  this  was 
considered  a  major  drawback. 


RECOMMENDATION  #5  Treat  your  model  building  project  tike  an 
actual  shipbuilding  project  in  the  preliminary  phase.  Account  for  weight  and  buoyancy 
to  the  greatest  accuracy  possible.  I  developed  a  weight  budget  for  Structure  (the 
frame),  Propulsion  (the  tow  pad),  the  Payload  (the  heave  post),  a  building  margin  of 
10%  and  a  design  reserve  buoyancy  of  10%  (allocated  because  there  is  some 
intention  to  eventually  convert  my  models  to  self  propulsion  so  there  may  be  some 
growth).  By  calculating  the  displacement,  I  then  knew  how  much  each  component 
could  weigh  and  I  did  W’eigh  them  frequently  through  the  building  process.  As  work 
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proceeded,  building  margin  was  allowed  to  be  utilized  at  a  rate  commensurate  with  the 
stage  of  building.  I  incorporated  this  method  early  enough  in  the  project  to  abandon 
P VC  as  a  building  material  before  too  much  work  was  done,  but  not  soon  enough  to 
avoid  buying  a  quantity  of  material  that  was  not  useful. 

RECOMMENDATION  #6  Make  frequent  tests  of  methods  you  think  you 
will  use,  well  in  advance  of  when  you  need  them  in  the  building  project.  If  you  plan 
to  use  a  particular  machine,  check  it  out  to  ensure  all  the  pieces  are  there  and  it 
functions  properly.  If  you  plan  to  use  a  particular  adhesive  to  bond  two  different 
materials,  try  it  out  early  on  with  two  scrap  pieces,  even  if  the  adhesive  claims  to  be 
formulated  to  do  the  job.  I  had  significant  difficulty  with  this  on  several  occasions. 
Testing  gave  me  the  time  to  either  call  the  manufacturer  for  assistance  or  to  find  a 
different  product.  The  two  times  I  did  not  perform  tests  on  intended  procedures 
resulted  in  near  disaster.  The  one  time  I  did  not  perform  a  bonding  test  was  with  PVC 
and  PVC  cement,  items  I  had  worked  with  in  the  past  and  was  assured  by  the  PVC 
supplier  would  work.  For  a  reason  I  never  determined,  the  cement  did  not  bond.  This 
occurred  late  on  a  Saturday  immediately  proceeding  the  Sunday  I  was  to  depart  for 
Annapolis.  The  search  for  a  new  bonding  agent  was  frustrating  since  it  was  not 
supposed  to  happen,  difficult  because  of  the  hour,  and  nearly  cost  me  a  significant 
piece  of  the  test  program  (the  high  speed  SWATH  model  did  not  require  significant 
reserve  buoyancy  and  so  was  built  from  PVC  to  utilize  the  material  purchased.).  I 
intended  to  paint  the  models  bright  yellow  with  marine  paint.  Yellow  is  the  color  of 
choice  for  model  testing  because  it  enhances  contrast  and  is  useful  for  photography. 
The  marine  paint  was  difficult  to  apply,  did  not  cover  well  and  was  slow  to  dry.  This 
provided  for  extreme  frustration  at  the  already  frantic  end  of  the  building  phase.  In  the 
end  analysis,  any  spray  enamel  with  a  10-15  minute  dry'  time  would  have  performed 
equally  well  in  the  tow  tank  tests,  and  I  would  have  known  that  if  I  tested  the  paint 
systems  in  February  rather  than  waiting  until  the  end  of  the  project. 

For  the  SLICE  hullforms.  molded  fiberglass  was  found  to  be  the  best 
construction  method.  With  other  materials  such  as  shaped  foam,  PVC  rod  turned  on  a 
lathe  or  stacked  wood,  construction  is  a  100%  hands  on  process.  Only  one  piece  can 
be  manufactured  at  a  time.  Each  SLICE  model  has  the  following  components:  4  struts, 
4  nose  cones,  4  tail  cones  and  4  parallel  midbody  sections.  There  were  two  SLICE 
models,  therefore  requiring  the  manufacture  of  32  components.  Hand  crafting  each 
component  would  not  only  have  been  prohibitively  time  consuming  but  also  would 
have  introduced  slight  variations  between  components  that  should  have  been  identical 
(the  four  nose  cones  for  example).  Molded  fiberglass  construction  offered  the 
following  advantages: 
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1.  A  mold  insures  that  the  four  components  have  the  same  shape. 

2.  Once  the  molds  are  made,  component  fabrication  time  is  primarily 
epoxy  curing  time.  This  frees  up  time  for  the  builder  to  do  something  else. 

3.  By  laying  in  the  fiberglass  cloth  in  the  correct  fashion,  excellent 
strength  is  obtained  for  very  little  weight. 

4.  Using  the  proper  additives  in  the  epoxy'  provides  an  outside  surface 
that  requires  no  finish  work  once  the  piece  is  removed  from  the  mold. 

The  West  Epoxy'  System  was  the  material  of  choice  for  this  project.  It  is 
relatively  expensive  but  is  readily  available  at  Marine  Supply  stores,  very  easy  to  use, 
and  is  associated  with  an  excellent  consumer  product  hotline  that  is  responsive  and 
helpful.  Most  of  the  products  referred  to  in  the  building  process  are  West  System 
names.  The  West  system  can  be  used  with  either  a  fast  setting  hardener  or  a  regular 
hardener.  The  fast  hardener  cuts  the  working  time  and  also  the  curing  time  of  the 
mixture  in  half.  Additives  which  are  called  for  in  some  of  the  steps  also  reduce 
working  time.  Times  cited  below  are  based  on  using  the  fast  hardener.  The  West 
System  sells  a  set  of  proportioning  pumps  which  helps  to  ensure  that  the  correct  mix 
ol  resin  and  hardener  is  achieved.  The  pumps  are  a  worthwhile  investment. 

Note:  The  epoxy  will  not  adhere  well  to  plastic  so  working  on  a  plastic  drop 
cloth  will  help  cleanup.  Mixing  the  epoxy  in  plastic  pails  can  be  reused.  Just  let  the 
epoxy  harden  and  the  tap  the  pail. 

The  following  outlines  the  steps  used  constructing  the  pods: 

1.  A  CAD  drawing  was  made  of  the  pod.  Five  or  six  tangent  lines  were  drawn 
along  the  curvature.  .Angles  with  the  vertical  and  insets  from  the  leading  edge  were 
calculated.  This  is  shown  in  figure  A-l. 

2.  A  male  plug  was  cut  out  of  PVC  rod  on  a  lathe.  This  was  done  by  marking 
the  insets,  and  then  cutting  at  the  appropriate  angle  from  the  inset  mark  to  the  leading 
edge.  The  transitions  between  angles  was  fared  with  a  file  and  the  entire  plug  polished 
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with  consecutive  applications  of  200,  320  and  400  grit  wet  sandpaper  followed  by  an 
abrasive  cleanser  and  finally  brass  polish.  .411  these  steps  were  performed  on  the  lathe 
to  ensure  roundness. 

3.  A  female  mold  was  then  made  from  the  male  plug.  This  was  done  with  the 
West  System.  First  the  Plug  is  coated  with  a  layer  of  wax.  This  is  necessary  to  ensure 
the  epoxy  does  not  adhere  to  the  plug.  The  first  epoxy  coat  consists  of  the  appropriate 
mix  resin  and  hardener  with  sufficient  microfibers  added  to  give  the  mixture  the 
consistency  of  paste.  Since  the  additives  decrease  the  working  time.  I  suggest  you 
pump  the  resin  into  the  mixing  pail,  add  the  additive  and  mix  it  thoroughly  and  then 
as  a  last  step,  add  the  correct  amount  of  hardener.  This  method  will  save  a  few 
minutes  of  working  life.  This  mixture  is  painted  onto  the  plug  and  serves  as  veiy 
smooth  gel-coat  layer.  Brushes  can  be  conserved  by  washing  them  with  acetone  after 
use  but  I  recommend  always  doing  this  first  layer  with  a  new  brush.  This  application 
is  allowed  to  begin  to  harden  but  not  cure  (approx  2-3  hours).  Even  with  steep  angles, 
this  mixture  should  be  thick  enough  so  there  is  very  little  running  of  the  coat  while  it 
sets.  A  second  layer  is  added  which  is  also  a  paste  like  consistency.  This  time,  the 
additive  is  a  50°o-50%  mix  of  microfibers  and  colloidal  silica,  which  will  be  less 
smooth  but  has  greater  strength.  This  also  is  set  aside  to  begin  to  harden  but  not  cure 
(2-3  hours).  Finally,  the  real  strength  layer  is  added.  For  the  molds,  this  consisted  of 
squares  of  fiber  matting  soaked  in  mixed  epoxy'  (no  additives  in  the  mixture  this  time). 
The  matting  has  a  random  fiber  lay  and  tends  to  be  fairly  (0.125"  thick).  Do  not  let 
the  matting  get  too  wet  or  else  it  will  tend  to  slide  when  placed  on  the  plug.  The  most 
effective  way  to  do  this  step  is  with  your  hands.  Latex  gloves  are  an  absolute 
requirement.  This  step  is  just  like  working  with  paper  mache.  When  the  matting  has 
begun  to  harden  (2-3  hours)  a  layer  of  12  ounce  woven  fiberglass  cloth  is  added  as  a 
final  step.  Pre-cutting  the  cloth  into  smaller  patches  will  help  since  your  gloved  hands 
will  be  covered  in  epoxy.  Since  there  are  no  additives  in  the  epoxy  in  the  last  two 
steps,  use  the  epoxy  sparingly,  just  enough  to  barely  wet  the  cloth.  .Allow  the  entire 
mold  to  cure  for  the  appropriate  time  (  8  hours  with  fast  hardener).  The  mold  may 
have  to  be  beaten  pretty  hard  to  pull  it  free.  In  the  case  of  my  shapes,  I  found  that 
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removing  the  plug  actually  drew  a  significant  vacuum  because  of  the  closeness  of  the 
fit.  The  mold  is  strong  enough  to  be  beaten  pretty  hard  with  a  rubber  mallet.  A  small 
hole  drilled  into  the  mold  and  later  patched  may  also  help. 

4.  Manufacturing  pieces  is  now  essentially  a  repeat  of  step  3  with  a  few 
exceptions.  I  chose  to  leave  out  the  fiberglass  matting  in  the  final  pieces.  There  was 
sufficient  strength  with  the  two  initial  coats  and  the  12  ounce  cloth  and  I  wanted  to 
reduce  weight.  Take  Note :  The  epoxy  cures  using  an  exothermic  process.  In  the 
steps  used  up  till  now,  this  was  insignificant.  Since  the  female  mold  is  a  cup,  as  you 
manufacture  the  component  pieces,  any  excess  resin  will  pool  in  the  bottom  of  the 
mold.  On  two  occasions  I  cracked  my  molds  because  of  the  heat  released  by  the 
excess  pool.  I  recommend  that  you  set  the  mold  into  a  pot  of  cold  water  at  least 
during  the  final  step  (glass  cloth  and  epoxy).  This  seems  to  be  adequate  to  control  the 
temperature.  Take  Note:  You  need  to  provide  a  method  of  pulling  the  piece  out  of 
the  mold.  I  did  this  by  setting  two  12"  pieces  of  nylon  string  into  the  wet  epoxy  of  the 
final  coat  and  covering  the  end  with  a  square  of  fiberglass  cloth.  This  gave  a  secure 
method  of  pulling  the  piece  out. 

5.  After  the  components  are  pulled  from  the  mold,  they  should  be  washed 
with  acetone  to  remove  any  residual  wax.  If  this  is  not  done  later  epoxying  during 
assembly  may  not  bond  well. 

6.  When  all  the  components  are  manufactured,  they  can  easily  be  assembled  by 
epoxying  them  together.  I  cut  rings  out  of  balsa  wood  that  were  fitted  into  the  joints  to 
increase  the  gluing  surface.  The  timing  of  this  assembly  is  dependent  on  how  the 
struts  get  attached.  For  the  SLICE  models,  the  struts  are  bolted  to  the  pods.  In  order  to 
facilitate  this,  the  tails  were  joined  to  the  parallel  midbody  sections.  The  struts  were 
then  attached  to  the  sub  assembly  and  finally  the  nose  cones  sealed  the  assembly. 


The  following  process  was  used  to  build  the  struts. 

1.  A  full  scale  CAD  drawing  of  the  strut  cross  section  was  made.  This  was  cut 
out  and  glued  to  a  piece  of  thin  plywood.  The  cross  section  shape  was  cut  out  of  the 
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wood  so  that  a  "slip  over  template  was  created. 

2.  Blocks  of  balsa  wood  were  cut  to  rough  dimensions  of  the  strut  and  glued 
together  to  create  the  correct  strut  length. 

3.  Once  dry.  the  balsa  wood  was  sanded  to  shape  so  that  the  template  made  a 
loose  fit  over  the  balsa  wood. 

4.  The  balsa  wood  was  sealed  using  a  coating  of  epoxy  painted  onto  it.  This 
was  left  to  cure. 

5.  The  strut  was  sanded  to  remove  large  defects.  (The  balsa  wood  will  cause  a 
fair  amount  of  bubbling  in  the  epoxy')-  It  is  then  covered  in  12  ounce  fiberglass  cloth. 
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Appendix  B 

Model  M-l  Baseline 
Table  of  Offsets 
and 

Design  Worksheet 


97 


Model  M-1  Table  of  Offsets 

Foward  Pod  in  Inches  Aft  Pod  in  Inches 


Struts  in  Inches 
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These  are  the  full  scale  SLICE  offsets  taken  from  the  December  1994  LMSC 
variant  used  to  generate  Model  M-l  offsets. 
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Model  M-1  Baseline  LMSC  SLICE 

This  is  based  on  the  drawings  provided  by  LMSC  for  the  December  baseline  ship. 

la.  Units  Defined 

ft=lL  IhslM  sec  =  iT  Base  Units 
nm  =6076  ft  Nautical  Mile 

kt  =  I  —  Knot 

lir 

Iton  =  22401b  Long  Ton 

1b.  Foward  Pod  Offsets  in  Full  Scale 

l  =(.)..  is 

Xfp.  -READ(mlpfx)  yf  =READ(mlptv) 

*  *  i 

Aft  Pod  Offsets  in  Full  Scale 

xap  =  RE  AD  (ml  pax)  y  =  READ(mlpay) 
i  n 

Strut  Offsets  in  Full  Scale 

j  =  u..  :T 

xs  --READ(mlsx)  v  =  READ(mlsv) 

U  'j 

Model  Scale  Factor 

>•  = — —  -  >.  =  17.258 

5.5625  in 
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Model  Offsets 


^fpvxipy  ~xfp./- 


.  - 1 


Y  fpv>'  ip j  =>'  fp.  ;;  '• 12  Foward  Pod 


Xapvxap/  “xap'/-  12  Y  ap\>  apj  “ >'  ap.  A  12  Aft  Pod 


x  'j  =x.  ;;‘  i2 

S\  z>  I  .v 

J 


Ysl>s  =ys.;.-'i2 


Struts 
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1c.  Model  Scale  Weight  Balance 

Fixed  Variables 

l  :=  I “.258 

2 

P  svv  =  1.9912  lb  — 

a4 
...2 

vsw  13343  I  O'5  — 
sec 


Density  of  Sait  Water  at  56°  F  An  assumed  average 
operating  temperature. 


Kinematic  Viscosity  of  Salt  Water  at  56°  F. 


pf\v  :  9367  ib  -j±_ 


density  of  r,esh  Water  at  68°  F.  An  assumed  average  tow 
:a.nK  temperature.  Use  this  if  preciction  is  for  a  model  test. 


v 


f\v 


1  0804  U)'5-  — 
sec 


Kinematic  Viscosity  of  Fresn  Water  at  68°  F. 


g  —  32. 1 744ft*sec 


Acceleration  due  to  gravity 


V  =0  kt.  1  kt..  35  kt 


Modify  this  to  reflect  the  velocity  range  of  interest. 


Indicator  Settings 

Water 


SW  ■  o 


M1"  inoicates  the  snip  is  operating  in  Salt  Water 
”0"  indicates  the  snip  is  operating  in  Fresn  Water 


P  =iiiSw=i,Psvv,Pfw; 


v  =  iltSWs  I ,  v  sw,v 
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Hull  Displacement  Definition 

...  _  170-lton  This  defines  the  the  reauireded  (d)  displacement  of  the  vessel 

>/'  note:  170  H°n  is  usea  here  sines  that  is  the  displacement  of  the  "filletless'4  hull. 

vd  ~^d'—  v  d  =  1  -it''  The  requireded  (d)  displaced  volume 

pg 

Bouyant  Force 

1%  =  V  j-p-g  Fb  =  74.084* lh 

Weight  Account 

FRAME 

W  --  l  x-  lb 

Force  Block  +  Heave  Post 

Wb,  =  25  lb 

Other  Model  Supplies 
Woth  -  2- lb 

Desired  Reserve  Bouyancy 

F  brd  -  h-06  F  b 
Fbrd  =4'445*ib 


Allowed  Weight  of  a  single  strut  pod 


W 

sp 


Fb-  W,--  Wb[-  '.V  0(h  -  Fbrd 


W  =6. 16*  lb 
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Appendix  C 

Model  M-2  Modified  Baseline 
Table  of  Offsets 
and 

Design  Worksheet 
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Model  M-2  Table  of  Offsets 


Pods  (both  fwd  and  aft  have  the  same  offsets) 
in  inches 


Struts 
in  inches 


xm(  x) 

i 

in 

ym(v) 

0 

0 

0.556 

0.8 

1.1 12 

1.145 

1.669 

1.41  1 

2.225 

1.632 

2.781 

1.82 

3.337 

1.983 

3.894 

2.1241 

4.45 

2.2461 

|5. 006  | 

2.351! 

5.562 

2.441 1 

6.1 19 

2.5181 

6.675 

2.5821 

7.231 

2.635 

7.788 

2.678 

8.344 

2.7 13  j 

8.9 

2.74  I 

9.456 

2.758! 

10  012 

2.771 

10.569 

2.779  j 

1  1.125 

2.781  i 

1  1.681 

2.78  | 

12.238 

2.7731 

12.794| 

2.763! 

13.35 

2.75  | 

13.906 

2.7351 

\m(-x).  -in 


1 4.463  i 

15.019 

15.575 

16.131 

16.687 

17.244 

17.8 

18.356 

118.912! 

1 9.469 1 

120.025 

20.581 

21.137 

21.694 

22.25 

22.806 

23.362 

23.919 

24.475 

25.031 

25.588 

26. 1 44 

26.7 

27.256 

27.812 

ynny)  -  in' 


2.716 

2.6941 

2.6691 

2.642 

2.61 

2-575  j 

2.535  j 

2.4911 

2.4411 

2.5X5! 

2.25! 

2.172 

2.083 

1.984 

1.873 

1.75 

1.614 

1 .463 

1.296 

1.113 

0.9091 

0.684| 

0.429! 

0 

xscale(xl  )•-  1 


0 

0.72 

1.44 

2.16 

2.87 

3.59 

4.31 

5.03 

p.75  | 

6.47  i 

7.19 

7.91 

8.62 

9.34 

10.06 

10.78 

1  1.5 

12.22 

12.94 

13.65 

14.37 

15.09 

15.81  j 

16.531 

17.25  j 

yscale(xl ) 
ft 

0 

0.287 

0.529 

0.729 

0.889 

1.009 

1.093 

1.144 

1.166 

1.168 

1.168 

1.168 

1.168 

1.168 

1.168 

1.168 

1.166 

1.144 

1.093 

1.009 

0.889 

0.729 

0.529 

0.287 

0 


Model  M-2  The  Better  Slice 


This  variant  uses  the  same  gross  dimmensions  as  the  baseline  SLICE  but  varies  the 
body  shapes  to  follow  more  traditional  forms.  The  struts  use  parabolic  leading  and 
trailing  edges  and  the  pods  are  based  on  series  58  shapes. 


la.  Units  Defined 

ftsIL  lbslM  secs  IT  Base  Units 

nm  =  6076  ft  Nautical  Mile 

kt  =  1  —  Knot 

hr 

lton  =  2240  lb  Long  Ton 

lb.  Fixed  Variables 


Psw 


=  1.99121b 


sec 


Density  of  Salt  Water  at  56°  F.  An  assumed  average 
operating  temperature. 


Vsw  -  I  3343' I0'5  — 
sec 


Kinematic  Viscosity  of  Salt  Water  at  56°  F. 


pfw 


=  i.  9367-lb- 


Density  of  Fresh  Water  at  68°  F.  An  assumed  average  tow 
tank  temperature.  Use  this  if  prediction  is  for  a  model  test. 


Vftv  =  1.0804  10‘5'  — 
sec 


Kinematic  Viscosity  of  Fresh  Water  at  68°  F. 


g  -  32. 174*ft*sec  2  Acceleration  due  to  gravity 


V  :=  0* kt,  1  kt..  35- kt  Modify  this  to  reflect  the  velocity  range  of  interest. 
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1c.  Indicator  Settings 

Water 

$$$  SW  -  I  *1*  indicates  the  ship  is  operating  in  Salt  Water 

*0"  indicates  the  ship  is  operating  tn  Fresh  Water 

p  -  if(SW=  I  .psw,Pf\v) 
v  .  -  if(S W—  I ,  v  sw>  v  faj 
Body  Segmentation 

$$$  SCg  :  =  1  •  ft 

2a.  Hull  Component  and  Required  Displacement  Definition 

Aj  =  1 77  lton  This  defines  the  the  requireded  (q)  displacement  of  the  vessel 

vd  V  d  =  6 1 89*ft3  The  requireded  (q)  displaced  volume 

Pg 

Strut  Height 

H  =  6-  ft  Ammount  of  Strut  in  the  Water  Defined  per  the  Baseline 
Strut  Beam 

D  j  =  3.25-ft  Maximum  Strut  Beam  Required  Defind  per  the  Baseline 
Considered  a  required  parameter. 

Overall  Strut  Length 

...  Li  524  ft  Overall  Strut  Length 


Length  to  Beam  Ratio 


LD1 


LD1  =7.385  Length  to  Maximum  Diameter  Ratio 
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Length  Breakdown 

Lpmbi  =7*ft  Length  of  Parallel  Midbody 

Lnl  ~a5'(L  1  '  Lpmbl)  Ln,  =8.5-ft  Strut  Nose  Length  (2.4  xd  is  optimum) 

L  tj  -  L  |  -  Lpmbi  -  L  ni  Ljj=8.5*ft  Strut  Tail  Length  (3.6xD  isoptimum) 

Nose  and  Tail  Shape  Exponent 

nnl  =2.25  ***  n  tl  =2.25 

Calculation  of  Nose  and  Midbody  Shape 

xl  :=0  ft,-seg..- fL-i  *  L^ui  ■  ,  ^ 

ni  pmoi  defines  the  x  values  of  the  subdivisions 

defines  the  y  values  of  the  nose  based 
on  a  parabola 

D1 

ypmb  1  ( x  1 )  -  —  defines  the  y  values  of  the  midsection 

2 

A I  ( xl )  =  iff  jxl  |  _L  n  j  ,ynl(xl  ),ypmbl(x  1 )  ] 

Calculation  of  Tail  Shape 

xl  0-ft,-  seg..  -  ^Lnj  -r-Lp^j  +Ltj^  defines  the  x  values  of  the  subdivisions 

i"  til 

j  defines  the  y  values  based 

I  j  on  a  parabola 

J  j 

bodvi(xl)  :=  iff  |  x  1 1  iLn|  +  Lpmb| ,  A 1  (x  1 ) ,  >t;  1  ( x  1 ))  Function  that  defines  the 

body  geometry 


vti(xl)  = — L-i  l 

2  I 


xl  -  (L, 


ynl(xl )  :=  — -• 
2 


/L  nl-  !xli\nnl; 


Strut  Section 


Body  I 


Wetted  Surface  and  Volume  Calculation 


Perimeter 

Nose 


1 - ;  D]!  .‘4  Ln,  -  ,/D  ,2- I6-Ln,2', 

N„  =05w'D|2Tl6Ln|V--i-lnl  1  1 - 1 

8'Lnl  '  u  1 


Np  =  17.55-ft 


Tail 


r 


I - - - -  D,2  4  Ltl  -r  ,,'D  ,2-  1 6  L t , 2 

Tp  -0.5-^D  ,2.  i 6  L t , 2 - —  Ini - ■t.?- •  1 - 1±- 


8L 


tl 


D 


Tp  =  1 7.55*  ft 
Midsection 

,C*p  “^Ppmbl 
Sp  =  !4-ft 
Total  Perimeter 

P  =Npr.Sp.Tp 
P  =  49.  l*ft 

Wetted  Surface  of  the  Strut 

WS$  ^  P  H  WS  s  =  294.602-l\2 


for  one  strut 


Areas 


nose  area 


Na  =f'Ln,-D 


tail  area 


Ta  =-.L„D| 

mid  section  area 
Sa  -D  ]  Lpmj,] 


N  a  =  1 8.4 17*  ft2 

Ta  =  18.41 7-ft2 
S  a  =22.75*ft2 


Total  Area 

Area  =Na+Sa-Ta  Area  =  59.583-ft2 


Volume  per  Strut 

vol  =  Area*  H  vol  =  357.5*ft3 

Total  Strut  Volume 

VOL  =4- vol  VOL  =  1430ft3 

Pod  Volume  Required 


v  pods  *  v  d  '  V0L  v  pods  =  4759'ft"'  for  four  pods 

V  d 

V pod  = vpod  =  ll90‘ft3  per  pod 

Pod  Definition  based  on  the  Model  41 55  of  Series  58 

^  1  =  4.3  the  scale  factor  used  to  get  from  the  data  on  page  49  of  senes  53  (Lpod=  9 

ft)  to  the  size  needed  per  baseline  lpod=  38.7  ft)  is  4.3 

Lp  ~  9- ft  X  ,  D  p  =  1.8  ft  X  j  D  p  =  7.74*ft 

i  =0, 1 ..  53 

x.  :=READ(betx)  y.  =  READ(betv)  non  dimmensional 

xfs(x)  =x-Lp  yfs(y)  ry-Dp 

WS  58  =  39.75  ft2 


1  10 


full  scale  dimmensions 


for  one  pod 


WS  fs  =  WS  58-/.  , 2  WS  fs  =  734.977  *  ft2 


Vol  -  1 4.89- ft-*  As  printed  in  the  Series  58  data 


Vol  fs  =  Vol  5g-/.  Vol  fs  =  1 1 83. 859* ft3 


VOLUME  -  4  Vol  ('s  -  4  vol 


VOLUME  =  6165.437* 


V  ^  —  6  1  88.72  /  *ft  Compares  calculated  volume 
to  needed  volume. 


A  =  VOLUME-pg 


A  -  176.334  ‘lion 


Now  Scaie  to  Mode!  Size 
5.5625-in 


'« 1 


D  , 


;.o  =0.06 


—  =  16.698 
7  i 


i  =0.1. .25  j  =26.27..  50 

xrrnx)  -\fs(\)  /0  vm(y)  =  yfs(yK  n 


Better  Slice  Strut  Offsets  a  0  =0.06 


xriA/.n 


xscaiet  x  1 ) 


yscale(xl)  =  body  1  ( x  1  )•  1 2  7 
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Weight  /  Volume  Balance 


Density  of  PVC 


D  r  --6.1875  111  lT  -  7.9375  m 
Dr  2 

V..  -j'i —  L  V  =0.1 38- ft'' 

*  0  5  1 


.7-lb 


ppvo  -  v 


i)  j\v  g  =n2.3i  Mb* it 


Ppvc  =S4.70S-:bMt  " 


diameter  and  length  of  rod 
weighed 

volume  of  1  foot  of  rod 

density  of  the  PVC 

density  of  FW  at  68°F. 


Displaced  Volume 


V 


model 


VOLUME 


V  model  Ej24M\ 


Bouyant  Force 


‘  b  v  mouei  f)  t\v  ^  F -r  ><2.523 •  lb 

Weight  Account 

FRAME 


Wf  =  15-lb 

Force  Block  +  Heave  Post 
Wb3  -  25-lb 
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Parallel  Midbody  X  4  pods 


„  _  5.5625  .  2  5.5625  -  0.5  .  2 '  . 

’  pvcpmb  _  “  in  ~  T  "  in:  ^pmb  4 

V  pvcpmb  =0-065*n5 

pmb  "  ^  pvcpmb  P  pvc  ^  pmb  “  -  482-lb 

Pod  Nose  ana  Tail  x  4  poas 

5  5o°5 

•  pcxir.t  :  v  pods'-  -  4L  pmb  ^ 1  in  :  volume  ofthe  nose 

and  tail 

;XHint  =0  645- ft3 


\  r 

’  nhorel 

1  (  2.25  in  )“•  1 .25  in  2 

V  nborel  =  39  761  •'n"' 

’  nbore2 

-  t(  2.0  in'r- 1  in-2 

Vnbore2  =25133  'in'' 

V  - 
nnore.i 

-  i  (  1 .5  in  r  2- in- 2 

^  nborej  ~  28  2  4  ‘,11 

v  nbore4 

-  t(  1.25-in  )“■  1 .0  in- 2 

^  nbore4  =  9  817  -in' 

r  :;bore5 

-  ti  0.75- in  )“•  1 .0-in-2 

V  nbore5  =  3  334  *in  ' 

^  tborel 

=  ::■(  2.25- in  f -2- in-2 

V  tborel  =63  617‘in'' 

^  tbore2 

=  *t(2.0-in j2  3 -in- 2 

V  tbore2  =  75398 ’in'' 

^  tbore3 

=  7t-(  1 .5-  in  r  2  in-2 

Vtbore3  =  28274 

X,/ :  bore  4 

-  n-(  1 .0  -in  T  2  in  2 

Vtbore4  =  '2.566  -in'' 

3  > 

v  :  bore 5 

*i  (  0.5  in  r  1  in  2 

vthore5  =  L57!'in'' 

11? 


V 


bore 


^  nborel  *  ^  nbore2  ^  nbore3  ^  ^  nbore4  ^  nboreo  ^  ^  tborel 
^  ’bore2  “  ^  lbore3  ~  ^  lbore4  T  ^  tbore5 


V  bore 


=  287.947  •in3 


V  pvepod  V  podnt  V  bore 

V  pvepod  ”^-479*ft 


W, 


podnt  ^  pvepod’ Ppvc 


W  podnt 88  40.547- ib 


4  Struts 

V  ^  -  Area  4  }.  0“-6  in 


V  ,  *0.427*n“' 


V  shore  KHn- l.X-in-6-in  t 

^  shore 

V  sbore2  4  ( 2. 75- in  0.7 -in- 6- in )•  2 

^  sbore2  ~ ( F  0 5  3  *  1 1 

^  pvcstr  ^  s  “  ^  shore  ~  ^  sbore2 

v  pvcstr  =  0.1 24*  ft3 

W  =  V  n 

stmts  pvcstr  pvc 

^  stmts  =  10.499Mb 
Other  Model  Supplies 
Woth  ^5  !b 

Checking  for  a  Balance 

F  ^r(j  =  1 0  lb  Desired  Reserve  Bouyancy 

W  model  -  W  f --  W  ^1  i-  W  pm'D  -  W  p0Cjnl  -  W  struls  *  W 
W  model  =  I -527Mb  Actual  Model  Weight 


oth 


Fbrn  =Fb-w 


model  Actual  Reserve  Bouyancy  (Should  =  Rbd) 


Fbra  =  -19  005  *lh 

i‘  bra  ''  brd  =  ~29.1K)5  •  lb  Must  be  0  or  greater. 


114 


Areas  to  Help  the  Balance 


reserve  =  areas  with  excess  weight  included  that  can  be  trimmed  down  by  the 
indicated  ammount 

r f  =  2-lb  frame  rbj  =  0  lb  force  block  rQlh  =  Mb  other  model  items 
reserve  =  rrrbrroth 

axed  =  W  fr-  W  bj  -  VV  ot|1  -  reserve  Weights  that  are  unchangable 

excess  =  reserve  bouyancy  that  is  desired  but  not  needed 
excess  -'Mb 

f  hr  •  F  xrd  -  excess 

Calculate  the  upper  limits  for  strut/nose  and  tail  cone  weights 


u  modelmax  ^  b  r  br 

u  modelmax  -  hxed 

yV - 

The  Max  Allowed  weight  for  a 

spa  ^ 

single  strut/pod  combination 

W  spa  =  7  63 1  *lb 

w  w  ,  w 

stmts  podnt  pnin 

4  4  4 

Summation  of  Allowed  Weights 


Wr  rt  =  13-lb 

Frame  Weight 

Wb|-rh|=25-!b 

Force  Block  and  Heave  Post 

Wo.h-  roth=4*lb 

Miscalleanous  Weight  Allowed 

W  spa  =  7-63  Mb 

The  Max  Allowed  weight  for  a 
single  strut/pod  combination 

Calculated  Weight  of  a  single 
strut/pod  combination 


Model  Dimmensions 


i 


16.68 

Lstrutm  =  Lr*m 

Lstrutm  =  1.439‘ft 

W.Sstrut  m  =WS  s  ;.m2 

WSstrut  m  =  1 .059* ft2 

LPod  m  r  Lp  m 

Lpod  m  =  2 . 3  2*  It 

wsPodm  =  wsfa.;.mJ 

WSpod  n)  =  2.642*  ft2 

116 


Appendix  D 

Model  M-3  Optimized  SWATH 
Table  of  Offsets 
and 

Design  Worksheet 
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Model  M-3  Table  of  Offsets 


Demihuli  Offsets 

in  inches  Strut  Offsets 

in  inches 


xh!2  ft* 1  -  ;20FF(xhi2';  ff 1  !2 


0 

0 

xl  ft- 1  - 

0.794 

0.667 

0 

1.589 

1.212 

0.794 

2.383 

1.641 

1.589 

3.178 

1.959 

2.383 

3.972 

2.173 

3.178 

4.766 

2.294 

3.972 

5.561 

2.335 

4.766 

6.355 

2.303 

5.561 

7.15 

|2.285 

6.355 

7.944 

i  2.268 

•  15  1 

8.738 

j  2.25 

7.944 

9.533 

2.232 

8.738 

10.327 

2.214 

9.533 

1 1.122 

2. 1 96 

10.327 

1  1.916 

2.178 

1  1.122 

12.71 

2.16 

1 1.916 

13.505 

2.142 

12.71 

14.299 

]2. 1 24 

13.505 

15.094 

j2. 107 

14.299 

15.888 

'2.089 

15.094 

16.682 

|2.07  1 

15.888 

17.477 

12.053 

16.682 

18.271 

12.035 

17.477 

19.066 

2.017 

! 

18.271 

19.86 

2.313 

1 9.066 

u 

21.449 

| 

2.814 

2.78 

2.78 

24.626 

2.78 

2.78 

26.215 

2.78 

27.01 

2.78 

27.804 

2.78 

28.598 

•  2.78 

29.393 

:  2.78 

bodyi(xl ) ■  ft' 1  •  12 

__ 

0.26 

0.488 

0.685 

0.852 

0.99 

1.099 

1.182 

1.24 

1.275 

1.29 

1.291 

1.291 

1.291 

1.29 

1.275 

1.24 

1.182 

1.099 

0.99 

0.852 

0.685 

0.488 

0.26 
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Model  M-3  The  High  Speed  Swath 


This  variant  uses  the  same  gross  dimmensions  as  the  baseline  SLICE  but  uses  a 
shape  characteristic  of  a  high  speed  swath.  Tandem  strut  design  was  selected  to 
minimize  strut  displacement. 

la.  Units  Defined 

ftslL  lb=lM  sec=IT  Base  Units 

nm  -  6076 Tt  Nautical  Mile 

kt  - 1—  Knot 

hr 

Iton  -  22401b  Long  Ton 

lb.  Fixed  Variables 

k  = 0  0662 

Density  of  Salt  Water  at  56°  F.  An  assumed  average 
operating  temperature. 

v  SW  2  1334310  -  Kinematic  Viscosity  of  Salt  Water  at  56°  F. 

sec 

Density  of  Fresh  Water  at  66°  F.  An  assumed  average  tow 
tank  temperature.  Use  this  if  prediction  is  for  a  model  test. 

V  f\V  =  '  ^804-10  -  Kinematic  Viscosity  of  Fresh  Water  at  68°  F. 

sec 

~2 

g  =  32. 1 74Mt*sec  Acceleration  due  to  gravity 

V  =  O  kt,  1  kt..  35  kt  Modify  this  to  reflect  the  velocity  range  of  interest. 


p  f.v  -  ;  9jO/  ib- - 

it4 


Psw  -  i  9912  lb 


1c.  Indicator  Settings 


Water 

$$$  SW  -0  "1”  indicates  the  ship  is  operating  in  Salt  Water 

"0"  indicates  the  snip  is  operating  in  Fresh  Water 


p.  =  ifiSW=l.Psw,pftvy 
V  =  if:  SWs  | .  y  sw,  v  ( 

Body  Segmentation 

$$$  Mig  =  lit). 

2a.  Hull  Component  and  Required  Displacement  Definition 

-d  =  1 77- Iton  /.'  This  defines  the  the  reauireded  (d)  displacement  of  the  vessel 

v  d  '  Ad'~  v  d  =  The  requireded  (d)  displaced  volume 

Pg 


2a. i.  Strut  Creation 

Strut  Height 

H  -  tS-  ill  Desired  to  allow  demi  hull  submergence  of  at  least  one  diameter. 

H  =  6.355  *m 

Strut  Beam 


Bj  :  3.25  ft;.  Maximum  Strut  Beam  Required  Defind  per  the  Baseline 
Considered  a  required  parameter. 

B  j  -  2.582  *in 

Overall  Strut  Length 

***  L  j  =  24- ft  /.  Overall  Strut  Length 
L  !  =  19.066 -in 
Length  to  Beam  Ratio 
L1 

LB1  = —  LB1  -7.3S5  Length  to  Maximum  Beam  Ratio 
B  i 


Length  Breakdown 


***  Lpm^j  -  3 ■  ft- X  Length  of  Parallel  Midbody 

***  Lnl  -0.5  :L  |  -  Lpmbl  Lnj  =  0.695*  ftStrut  Nose  Length  (2.4  x  d  is  optimum) 

Ljj  -  L  j  -  L  pm^  j  -  L  n  j  Ltj  =  0.695*ft  Strut  Tail  Length  (3.6  x  d  is  optimum) 


Strut  Nose  and  Tail  Shape  Exponent 


nns.2.25 


n,.  =  2.2o 


Calculation  or  Nose  and  Midbody  ShaDe 
xl  =  U  li.  sea..  L„|-L  hi 


defines  the  x  values  of  the  subdivisions 


vn  1  ( \  I )  -  — —  i  _ 


3,  Lnj-  xl 


'  nl 


vpmb  1  ( \  I )  -  • 


defines  the  y  values  of  the  nose  based 
on  a  parabola 


defines  the  y  values  of  the  midsection 


A 1  ( x  1 )  -  if  x  1 ;  l  L  n  j ,  yn  1  ( x  1 ) ,  ypmb  1  ( x  1 ) ) 
Calculation  of  Tail  Shape 


xl  =  Oft.-  sec..  Lnj  -  Lpm^j  -  Ltj  defines  the  x  values  of  the  subdivisions 


13  , 

ytl(xl)  --  i 


^  n  i  1'pmbl 
.  Ui 


defines  the  y  values  based 
on  a  parabola 


body  1  ( x  l’i  -  if  'xl!  _Lnj  ~  Lrmhl ,  Al(xl  ),ytl(xi )  Function  that  defines  the 

body  geometry 
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2.a.ii.  Wetted  Surface  and  Volume  Calculation  For  the  Struts 


Perimeter 

Nose 


N_  - ' )  5-  B 

r 

N  =  I  42  !*tt 


16  L 


3 


ill 


8L 


-  In  i 


4-Lnr.Br 


nl 


Tail 


B,2  '  4-L 

— 1 — Ini - 


tl 


SL, 


Tp  =  1.421-rt 
Midsection 

e  -  ->  t 

‘  P  ~  ~  ^  pmb  1 
sp  =0.397-ft 

Total  Perimeter 

p  =np--vtp 

P  =3.24-11 

Wetted  Surface  of  the  Strut 

WS  -PH  W.S  c  =  l.716*fi‘  for  one  strut 


Areas 


nose  area 


Na  =fL„l-B 


Na  =0.1 -It2 


tail  area 


T  a  =f  LtrB  1 


Ta  =0  i-ir 


mid  section  area 

a  "  ^  1  '^pmbi 


sa  =0.04.1 -ft2 


Total  Area 


\rea  \  •  S  a 


Aren  =0  242*11" 


Volume  per  Strut 
voi  -Area-H 


ml  =0.128*11-' 


Total  Strut  Volume 
VOL  s  =  4-vol 


VOL  s  =0.51 3*fV 


2.a.iii  Demi-Hull  Creation 


Demi  Hull  Volume  Required 


vdh  =VJ- V°LS  vJh  =  ! 


V  dh 


=  0.7  -IV' 


for  two  hulls 


per  demihull 


Demi  Hull  Shape 


D  |,  =  7ft/.  Max  diameter  of  the  Hull  D  ^  =  5.56 1  -in 


D  hi  =0-84-Dh 

0^1=  4.67  1  *in 

Section  #1-2  interface 

D  h2  =  73  0  h 

D  =  4.059  *in 

Section  #2-3  interface 

D  h3  ' D  h 

D  ^  =  5.56 1  *in 

Section  #3-4  interface 

D  h4  - D  h 

D  j  ^  =  5.56 1  *m 

Section  #4-5  interface 

Dhs  -;l73Dh 

D  =  4.059 -in 

Section  #5-6  interface 

Dh(,  s.:  D, 

I)  Kr,  =4  -::i 

Section  #6-7  interface 

=  90-rt- /. 

Overall  Length  of  the  hull. 

L  h  =5.958*11 

L  hi  =«-08-Lh 

Lhi  =5-72,in 

Section  #1 

L  h2  z  !9  Lh 

L  |l2  =  13.584  -in 

Section  #2 

L  h3  =il05Lh 

L  h3  =  3.5“5  *in 

Section  #3 

L  h4  =0  25  Lh 

L  =  1 7.87 4  *m 

Section  #4 

Lh5  =!)'>Lh 

L  =  10.724  *m 

Section  #5 

L  h6  :;)09  Lh 

Lh6  =!r435,in 

Section  #6 

L  Ii7  -‘!)|l>Lh 

L  =  1  3.584  *m 

Section  #7 

I-h 

1.2  D  L2D  =  12.857 

19  h 

Segment  1  &  2 

Calculation  of  Nose  Cone  Shape 

n  i  c  o  25 
n  1 

Nose  cone  shape  exponent 

\h!2  =  0-ft.-seg..-  LhrLh2J 


defines  the  x  values  of 
the  subdivisions 
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the  subdivisions 


yh  i  (\h 1 2) 


yh2(\hl2 j 


phl  j,  Lh|-  xhi2.\"nl 
2  Lhl 

D  hi  D  h2“  ^  hi  ,  r  ,  ,  . 


defines  the  y  values  of  the 
nose  based  on  an  ellipse 

defines  the  y  values  of 
see  ]  j  section  2 


Al(xhl2)  =  if  I xii  1 2 i  < L h , , y h I f xli 1 2 ) , yh2 ( xii ! 2 / ' 


offsets  for  the  first  two  body 
segments. 


Adding  in  Segment  3 


\bl2  =  O  ft,- seg..  -  Lh,  -  Lh2-L!l3/ 

D  to  D  -  S  i,-i 

vh3  f  xh  1 2 ;  -  — -  -  xh  1 2 

■>  ->  T 


defines  the  x  values  of  the 
subdivisions 

to  -  see  :  defines  the  y  values 
of  section  3 


BK\hi2)  -  if  xh  1 2  i  <Lhl  -  Lh2,AI(\hl2),yh3i  \hi2) 


offsets  for  the  first  three 
body  segments. 


Adding  in  Segment  4 

defines  the  x  values  of  the 
subdivisions 

defines  the  y  values  of 
section  4 


xii  i  2  ufi  ,-sec..  Li 


yh4(\hl 2) 


13 


h3 


C 1  ( xh  1 2 )  =  if  xhi2!<Lh]  -Lh2-Lh3lBlfxhi2\yh4(\lil2Vi 

offsets  for  the  first  three  body 
segments. 

Adding  in  Segment  5 


\h!2  -(Mt, -see..  I' hi  -  L  h2  “  L  h3  *•  L1i4  -  L 


,  -  ,  ,  ~ ^  h4  ^  ho  li4  ,  - 

yho(\hl2)  - - - -  xii  1 2  L  u 

2  ^  L  t  - 

“  L  ho 


defines  the  x  values 
of  the  subdivisions 

L  h2  -  L  h3  T  L  h4  ~  se? 


defines  the  y  values  of 
section  5 


0FF(xhl2) 


DI(xhl2)  =  if  ;.\]il2;<Lhl  ^Lh2^Lh3-rLh4,Cl(xhl2),yh5(\lil2)j 


offsets  for  the  first  three  body 
segments. 


Adding  in  Segment  6 


xhl2  -Olt.-seg..  -  Lhl  T  Lh2  ~  Lh3  ~ ‘'h4*  Lh5  *  Lh6/  defines  the  x  values  of  the 

subdivisions 


yh6(xhl23 


^  h5  ^  I16  ^  h5  i  - 


I  ^irl  2  -r  L  u  |  -  L  uo  -  L  -  L 


'  li6 


h2~  Lh3~Lh4  +  Lh5  ;  ~  se§] 


defines  the  y  values  of 
section  6 


El(.\hl2)  \hi2  <L h j  -  L h2  -  L h3  -  L  h4  -r  L h5,DI( xh!2),yh6(xhl2) ) 


offsets  for  the  first  three  body 
segments. 


Adding  in  Segment  7 


vhn  T  r  t  i  i  t  t  defines  the  x  values 

mu 2  u  n.-aes..-  L hi  -Lh2*Lh3^Lh4"Lh5-Lh6~  Lh7i 


of  the  subdivisions 


yh7 ( \h i 2 ) 


D  h6  ‘  D  h6 

2  i  - 

-  L  h  / 


\h  1 2  -  L  h  j  -  Lh2-  Lh3-Lh4>Lh5tLh6  -  segj  j 

defines  the  y  values  of 
section  7 

OFF( xh  12'.  - :!'.  xhl2!  <L  hl  -  L h2  -  L  h3  T  L  h4  T  L h5  *  L  h6,El(xhl2).yh7(xhl2)] 


offsets  for  the  first  three 
body  segments. 


0.4 
0.2 
o 

-0.2 

”6  -5  -4  -3  -2  -1  0 

\h!2 


126 


Wetted  Surface  of  TWO  demihulls 


r  Lh 

WSh  -  2 •  ■  ■  2-wOFF(xlil2)  dxli 1 2 :  WS  h  *  13.2- ft2 

-  0  ft 


Volume  of  TWO  demihulls 


VOLh  -2-:  it-OFF(xhl2)2  dxli  1 2 VOL  h  =  !  ,3-ft3 

'0  ft  ! 


Deck  Area  Lost  compared  to  SLICE 

C  =3  ft  This  is  the  combined  ammount  the  bovv/stern  are 
cantelivered  over  the  ends  of  the  hulls. 

Lh 

105  ft — —  -  C  :  55  ft  =  660* It"  (negative  implies  a  gain  in  deck  area) 

/. 
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Model  Building  Balance 


Bouyant  Force 

Fb  avail  Fb  =  112.496-lb 

Weight  Account 

FRAME 

Wf  =  18  lb 

Force  Block  +  Heave  Post 

W  hl  -  25  lb 

Other  Model  Supplies 

W0.h  =  2-lb 

Desired  Reserve  Bouyancy 

Fbrd  =00S-Fb 
Fb|.d=  5.625- lb 


Allowed  Weight  of  a  demihull  +  2  Strut  combo 


W*P 


Fb- Wf-Wbl-Woth-Fbrd 


W  =30.936- lb 
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SLICE  Project  Test  Plan  for  the  Week  of  6-10  March  1995 


Anticipated  Number  of  Models  to  Test:  4 

Desired  Data  to  Obtain:  Resistance  (x)  for  all  4  models.  Dynamic  Sinkage  (z)  for  2 

models. 

Model  Descriptions:  Model  M-l  is  the  LMSC  baseline  design. 

Model  M-2  anticipates  improvements  to  the  LMSC  baseline  by 
improving  the  shape  of  the  bodies  but  retaining  the  gross 
geometry  of  the  baseline  ship. 

Model  M-3  is  an  optimized  SWATH  based  on  the  results 
presented  by  McGregor 
Model  M-4  is  the  Mid  Foil  Concept. 

Desired  Testing  Schedule: 


6  March  95  Set-Up  M-l;  M-l  Test  1  20runs 

7  March  95  M-l  Test  2;  Set-Up  M-2;  M-2  Test  1  30  runs 

8  March  95  M-2  Test  l(fmish)  ;  Set-Up  M-3;  M-3  Test  1  30runs 

9  March  95  M-3  Test  1  (finish);  M-3  Test  2  30  runs 

10  March  95  Set-Up  M-4;  M-4  Test  1  20  runs 


Test  Description: 

Test  1:  Resistance  Locked  in  Heave  and  Pitch. 

Objective:  To  measure  forces  in  the  X  and  Z  directions  over  the 
operating  range  of  speeds. 

Discussion:  Previous  testing  has  indicated  that  a  significant  Z  force  is 
generated  and  that  the  magnitude  of  the  force  is  geometry  dependent.  This  method  of 
testing  replicates  the  LMSC  methodology  and  therefore  allows  for  a  basis  of  model 
comparison  with  LMSC  results.  The  LMSC  justification  for  this  methodology  is  that 
control  surfaces  will  be  employed  to  fix  the  operational  draft. 

Description:  With  the  model  rigidly  fixed  to  the  carriage,  the  following 
sequence  of  speeds  will  be  used  to  collect  data  for  both  X  and  Z  forces  :  Model  Scale 
Equivalent  (MSE)  of  3,6,9,12,16,20,25,30,35,40  kts  for  basic  data.  Each  speed  will 
accomplished  with  one  run.  Then  MSE  of  20,25,30,35,40  kts  to  demonstrate 
repeatability.  If  time  in  the  day  permits  then  additional  runs  will  be  made  at  MSE  of 
10,11,13,14  and  22  kts  to  fill  in  hump  data  and  the  high  speed  knuckle.  This 
represents  a  total  of  20  runs  for  the  test  completion. 


Test  2:  Free  to  Heave. 


Objective:  To  quantify  the  difference  in  resistance  results  for  a  model 
rigidly  fixed  and  one  that  is  allowed  to  oscillate  about  a  desired  waterline. 

Discussion:  Although  control  surfaces  may  hold  the  ship  to  a  fixed 
depth  in  a  gross  sense  of  the  word,  it  is  unlikely  that  the  system's  sensitivity  will 
exactly  replicate  a  rigid  mounting.  This  test  will  use  a  counterbalance  to  replicate  the 
control  surface  but  allow  the  model  to  oscillate  about  this  point. 

Description:  This  test  will  use  the  Z  force  results  from  Test  l  to 
determine  the  counterbalance  weight  required.  The  model  will  be  attached  to  the 
carnage  free  to  heave  but  locked  in  all  other  respects.  The  necessary  counterbalance 
weight  be  added  and  the  model  towed  at  the  given  speed.  (Note,  due  to  the  small  TPI 
associated  with  SWA  TH  vessels ,  a  method  for  adding  the  weight  as  the  carriage 
comes  to  speed  may  be  necessary.  This  is  a  Discussion  Point  for  the  3  Feb  meeting.) 
The  test  ng  will  require  a  positive  stop  in  the  negative  z  direction  with  a  method  of 
determining  that  the  stop  has  been  reached.  This  stop  must  be  set  to  no  more  than  3 
inches  above  the  design  waterline  of  the  model. 

The  following  sequence  of  speeds  will  be  used  to  collect  data  for  X 
force  only  with  note  being  taken  of  the  counterbalance  weight  added.  Model  Scale 
Equivalent  (MSE)  of  9,12,16,20,25,30,35,40  kts  for  basic  data.  Each  speed  will 
accomplished  with  one  run.  Then  MSE  of  20,25,30,35,40  kts  to  demonstrate 
repeatability'.  If  time  in  the  day  permits  then  additional  runs  will  be  made  at  MSE  of 
3,  6,10,1 1,13,14  and  22  kts  to  fill  in  hump  data  and  the  high  speed  knuckle.  This 
represents  a  total  of  20  runs  for  the  test  completion. 
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Appendix  F 

Model  Test  Results 
Raw  Data 
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Index  of  Test  Runs  Conducted  6-10  March  1995 


Run  #  Model/Test 

Full  Scale 
Equivalent  \, 
_ (Kts) 

Test 

Velocity 

(EPS) 

Drag 

(LW) 

Lift 

(Lbf) 

1  I  M-1  in  Baseline  Confiq. 

10 

4.08i 

6.46 

1.27 

2  i  M-1  in  Baseline  Confiq. 

3; 

1.21  ! 

0.151 

0.191 

M-1  in  Baseline  Confiq. 

6  i 

2.41  : 

0.689 

0.562 

M-1  in  Baseline  Confiq. 

9: 

3.64 

4.24 

1.09 

3  i  M-1  in  Baseline  Confiq. 

12 

4.85! 

9.08 

2.42 

M-1  in  Baseline  Confiq. 

13.5 

5.47 

6.16 

3.33 

6 !  M-1  in  Baseline  Confiq. 

16! 

6.47 ! 

5.61 

5.18 

7 ;  M-1  in  Baseline  Confiq. 

20 

8.08! 

7.19 

7.73 

R7  M-1  in  Baseline  Confiq. 

gran 

8 1  M-1  in  Baseline  Confiq. 

■Bjuifl 

9  i  M-1  in  Baseline  Confiq. 

30. 

RS  M-1  in  Baseline  Confiq. 

30 

12.13! 

19.39 

15.46 

10 !  M-1  in  Baseline  Confiq. 

35 

14.14 

28.82 

11  M-1  in  Baseline  Confiq. 

<10 

16.18 

35.32 

26.42 

12  :  M-1  in  Baseline  Confiq. 

10 

4.05! 

6.21 

1.02 

M-1  in  Baseline  Config. 

11 

4.45 

8.91 

1.64 

14 1  M-1  in  Baseline  Confiq. 

13 

5.25 

6.86 

2.74 

M-1  in  Baseline  Confiq. 

14- 

5.66 

5.65 

3.57 

16  i  M-1  in  Baseline  Confiq. 

15 

6.02! 

5.48 

4.25 

17 '  M-1  in  Baseline  Confiq. 

40 ; 

16.18! 

35.38 

26.65 

181  M-1  with  Spray  Rail  Aft 

30 

12.13! 

19.45 

15.48 

191  M-1  with  Spray  Rail  Aft 

35 

14.14; 

27.95 

20.61 

20  i  M-1  with  Spray  Rail  Aft 

40 

16.18l 

36.87 

26.87 

24 1  M-2  in  Baseline  Confiq. 

10 

4.11  : 

6.83 

0.443 

25 !  M-2  in  Baseline  Confiq. 

3, 

1.22! 

0.172 

0.247 

M-2  in  Baseline  Confiq. 

6 

2.46 

0.667 

0.845 

27 1  M-2  in  Baseline  Confiq. 

9: 

3.68 

4.52 

0.416 

M-2  in  Baseline  Confiq. 

12 

4.93 

9.92 

2.37 

28  1  M-2  in  Baseline  Config. 

6 

2.46 

0.679 

0.724 

29 1  M-2  in  Baseline  Confiq. 

12: 

4.93 

9.76 

2.51 

30  i  M-2  in  Baseline  Confia. 

16 

6.58 

6.28 

5.95 

31 !  M-2  in  Baseline  Confiq. 

16:' 

6.58- 

6.3 

5.89 

32 1  M-2  in  Baseline  Confiq. 

20 1 

8.23! 

7.25 

9.35 

R32  M-2  in  Baseline  Confiq. 

20 

8.23 

7.26 

9.29 

34 1  M-2  in  Baseline  Confiq. 

25 

10.28 

10.23 

12.97 

R34  M-2  in  Baseline  Confiq. 

25 

10.28 : 

10.19 

13.03 

36 1  M-2  in  Baseline  Confiq. 

30 

12.35  ; 

14.54 

17.52 

30 

12.35! 

14.53 

17.48 

35! 

14.4; 

20.17 

22.33 

»— T1  ■  II II  ra— "Bara 

35 

14.4 

20.19 

22.3 

gnraBj 

40 

16.45 : 

26.42 

28.25 

40 

16.45 

26.36 

28.35 

42 ;  M-2  in  Baseline  Confiq. 

11 

4.47! 

9.81 

0.77 

11.5 

4.72! 

10.78 

1.82 

13.5 

5.47' 

6.68 

3.25 

45 ;  M-2  in  Baseline  Confiq. 

12.5 

5.22 

7.85 

2.86 

46  i  M-2  in  Baseline  Config. 

14.5 

5.97 

6.05 

4.27 
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M-1  In  Line  Confi 
M-1  In  Line  Confi 


M-1  In  Line  Confi 
M-1  In  Line  Confi 
M-1  in  Line  Confi 
M-1  In  Line  Confi 
M-1  in  Line  Confi 
M-1  In  Line  Confi 
M-1  In  Li 
M-1  In  Line  Co 


11.48 

22.29 


35.17 

44.51 


5.36 

7.54i 

10.26: 

14.22: 


19.67 

26.45 


TIME(SEC) 


TIME(SEC) 


TIME(SEC) 


TIME(SEC) 
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TIME(SEC) 


TIME(SEC) 


152 


TIME(SEC) 


TIME(SEC) 


154 


VELOOO  012(FPS)-ALL  T229/S0 


TIME(SEC) 


Csl 
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TIME(SEC) 
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TIME(SEC) 
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TIME(SEC) 


TIME(SEC) 


TIME(SEC) 


161 


TIME(SEC) 
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TIME(SEC) 


TIME(SEC) 


IME(SEC) 
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TIME(SEC) 


166 


IME(SEC) 


167 


TIME(SEC) 


168 


TIME(SEC) 


169 


TIME(SEC) 


TEST: Model  M-2  lest  I  RUN  N():036 
05/30/95  13:45:37 


170 


TIME(SEC) 


171 


TIME(SEC) 


CN 
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TIME(SEC) 


TIME(SEC) 
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TIME(SEC) 


175 


TIME(SEC) 


TIME(SEC) 


177 


TIME(SEC) 


TIME(SEC) 


179 


(oas)Tiwi.i. 


TIME(SEC) 


VELOOO  048(FPS)-ALL  T229/S0 


182 


TIME(SEC) 


VELX)00.053(FPS)-ALL  T229/S0 


TIME(SEC) 


i-ALL  T229/S0 


VELOOO  057(FPS)-ALL  T229/S0 


187 


VEL000.058(FPS)-ALL  T229/S0 
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TIME(SEC) 
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TIME(SEC) 


190 


TIME(SEC) 
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TIME(SEC) 


TIME(SEC) 
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TIME(SEC) 
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TIME(SEC) 
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TIME(SEC) 
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TIME(SEC) 


Appendix  G 


Data  Scaling  Worksheets 
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This  Calculation  Scales  the  Test  Data  obtained  for  Model  M-1  to  Full  Scale . 


NOTE:  The  System  of  Subscripts  in  this  Routine  can  become  rather  bulky...  *R\  *C".  “M",  ‘S’  are  Resistance,  Coefficient 
of  Resistance,  Measured  and  Ship  respectively...  where  measured  refers  to  the  measured  model  (Model  M-2  Data 
unsealed)  and  Ship  refers  to  the  full  scale  calculation.  Also,  "s",  "p",  "r",  T,  "t"  as  subscripts  are  strut,  pod,  residual, 
frictional  and  total  respectively.  Other  uses  of  symbols  are  defined  as  they  arise. 

Input  Data  Section 

1  z0..  16  Sets  the  data  range 

RM  ^  j  • =  READ(m  1  tl  fx)  Total  measured  X  component  force  for  test  #1 :  fixed  in  pitch  and  heave, 
i 


VM.  :=READ(mltlv) 


\  =  17.26 


VMlcts.  =VM, 


1 


1  1.688 


VSkts.  =  ^kts  -A 


Speed  increments  used  in  the  test,  (inputted  in  fps) 


Scaling  Factor  for  Lengths 


Measured  Velocity  in  Knots  (converts  the  inputs  which  are  in  fps) 


Scaled  Velocity  in  Knots 


VS.  :=  VU.-Jl 

1  »  V 

:=  1.2260- 10' 5 
osw  :=  1.1794- 10’ 5 
Pf\v  =  1-9384 
psw:=  1.9890 


Equivalent  Scaled  Velocity  in  FPS 


ftA2/sec  Kinematic  Viscosity  of  Fresh  Water  at  59  F  (Test  Temp) 


ftA2/sec  Kinematic  Viscosity  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


lb  secA2/ftA4  Density  of  Fresh  Water  at  59  F  (Test  Temp) 


lb  secA2/ftA4  Density  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 
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Geometry  Section 


LM  s  :=  1.39 

SM  s  =0.715 

LMp  =2.14 
SMp  =2.01 

SMt  :  =  4  SMs+4  SMp 


feet  Length  of  one  Start 


ftA2  Wetted  surface  of  one  Strut 


feet  Length  of  one  Pod 

ftA2  Wetted  Surface  of  one  Pod 


ftA2  Wetted  surface  of  the  entire  underwater  body 


Resistance  Calculations 

NOTE:  In  this  section  Rn  is  calculated  for  the  pod  and  the  strut  separately  since  their  lengths 
are  different.  Coefficients  of  friction  are  then  calculated  separately  for  each  based  on  the  1 957 
ITTC  Une.  These  are  converted  to  Resistance  Values  and  then  subtracted  from  the  Total 
Resistance  value.  Finally  this  is  converted  to  a  Coefficient  of  residuary  resistance  which  is  equal 
for  the  model  and  the  ship. 


VM.LM  s 

RnM  ~  - - - -  operating  Rn  of  the  Strut 

J  Ufw 


RnM  p 


VM.LM  n 
■  P 

uf\v 


operating  Rn  of  the  Pod 


CM 


fs. 


0.075 


^log^RnM 


Coefficient  of  Frictional  Resistance  based  on  Strut  Length 
and  using  the  1957  ITTC  line. 


CM 


0.075 

^ log /RnM  -  2j2 


Coefficient  of  Frictional  Resistance  based  on  Pod  Length 
and  using  the  1957  ITTC  line. 
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Frictional  Resistance  of  four  Struts 


RMfsi:=CMfs.~Pft,(VMi)J4SMs 


RM^J  i-CM^ — •  p  ^VMJ2-4- SM  p  Frictional  Resistance  of  four  Pods 


RMr  :=RMU  -RMfg  -RMf 

i  i  i 


Residual  Resistance  for  the  Model  as  measured. 


CM, 


RM, 


j"Pfw  (VM.)2SMi 


RnM , 


RnM 


Pi 


RM 


tl. 


RMfs  RM 


fPi 


RM, 


CM 


fs. 


CM  +•  CM, 


lfPi 


137186 

273238 

412692 

459176 

462577 

504527 

530604 

549878 

595228 

620171 

641713 

733548 

916085 

1 146240 

1375261 

1603148 

1834437 

211207 

420669 

635367 

706933 

712170 

776754 

816900 

846574 

916395 

954796 

987961 

1129347 

1410375 

1764715 

2117308 

2468157 

2824241 

0.151 

0.689 

4.24 

6.29 

6.46 

8.91 

9.43 

9.08 

6.86 

6.16 

5.65 

5.61 

7.18 

10.645 

19.465 

28.82 

35.35 

0.031 

0.102 

0.211 

0.254 

0.258 

0.3 

0.328 

0.35 

0.402 

0.432 

0.459 

0.582 

0.865 

1.29 

1.786 

2.351 

2.994 

0.077 

0.258 

0.535 

0.647 

0.655 

0.765 

0.836 

0.891 

1.026 

1.104 

1.173 

1.489 

2.216 

3.312 

4.595 

6.057 

7.723 

0.043 

0.328 

3.494 

5.389 

5.547 

7.845 

8.265 

7.839 

5.432 

4.624 

4.017 

3.538 

4.099 

6.043 

13.084 

20.413 

24.633 

0.008 

0.006 

0.006 

0.006 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 

0.007 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 

0.004 

0.004 

0.003 

0.005 

0.025 

0.031 

0.032 

0.038 

0.036 

0.032 

0.019 

0.015 

0.012 

0.008 

0.006 

0.006 

0.008 

0.01 

0.009 

Note:  these  values  for  Rn  are  Low  for  ITTC  1957  use.  Should  be  on  order  of  10*7 
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Scaled  Resistance  Calculations 


CS  r  :  -  CM  r  Residual  Resistance  Coefficient  remains  the  same  in  the  saling. 


VS.-LM 

i 


operating  Rn  of  the  Strut  in  the  new  scale 


VS.-LM  n-Jt 
i  P 


operating  Rn  of  the  Pod  in  the  new  scale 


l  logj  RnS  \  -  2 


V  I 


Coefficient  of  Frictional  Resistance  based  on  scaled  Strut  Length 
and  using  the  1957  ITTC  line. 


(log/RnS  \-  2) 


Coefficient  of  Frictional  Resistance  based  on  scaled  Pod  Length 
and  using  the  1957  ITTC  line. 


CSQ  :  =  0.0004 

a 


Correlation  Allowance 


j  2  2 

RS  =  CS  —  *p  sw  f  VS.  j  -4-  SM  Frictional  Resistance  of  the  Struts 

i  i  2  ”  '  ‘5 


j  2  2 

RS  =  CS  —  -p  svv  ^VSj)  -4-SMp-X  Frictional  Resistance  of  the  Pods 


RS  r  -  CS  r  — 'P  sw‘  (VS.^2-SM  |-X2  Scaled  Residual  Resistance 


RS  a  -  CS  a  —  p  svv*  (VS.)  SM|  A  Correlation 


RSt.;=RSfs.-RSfp.  +  RSr  +  RSa. 

i  i  r  i  i  i 
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Scaled 

Scaled 

Total  M-2 

Speed 

Speed 

Resistance 

(fps) 

(knots) 

at  Full 

Scale 

VS. 

VSkts. 

I 

RS  t 

i 

Residual  PodFrictional  StrutFrictional  Reynolds  #  Reynolds  # 

Component  Component  Component  of  the  Strut  of  the  Pod 


2.978 

5.932 

8.959 

9.968 

10.042 

10.952 

11.518 

11.937 

12.921 

13.463 

13.93 

15.924 

19.887 

24.883 

29.854 

34.801 

39.822 


488.89 

2680.46 

20480.51 
30929.77 
31798.85 
44370.38 

46882.96 
44860.57 

32717.51 

28777.95 

25868.95 
24672.03 
30738.84 
45766.63 
88587.27 

133802.96 
163618.18 


225.11 

1732.1 

18434.19 

28431.8 

29266.18 

41391.37 

43609.44 

41361.02 

28658.3 

24394.36 

21195.81 
18667.49 

21627.81 
31883.98 
69030.12 
107699.01 
129965.39 


167.14 

592.83 

1268.39 

1544.88 

1566.09 

1838.72 

2018.35 

2156.07 

2496.66 

2693.8 

2869.63 

3676.58 

5551.69 
8417.37 

11809.69 
15708.47 
20189.01 


63.99 
226.03 
482.53 
587.39 
595.44 
698.78 
766.85 
819.03 
948.04 
1022.69 
1089.26 
1394.66 

2103.76 
3186.42 

4466.99 

5937.76 
7627.04 


10225869 

20367227 

30762119 

34227083 

34480617 

37607536 

39551296 

40987988 

44368441 

46227690 

47833405 

54678821 

68285143 

85440940 

102512226 

119499001 

136739309 


p. 

ri 

15743425 

31356738 

47360385 

52694933 

53085266 

57899371 

60891923 

63103809 

68308247 

71170688 

73642796 

84181783 

10512964 

13154216 

15782457 

18397687 

21051951 


WRJTE(mlfulfx)  -RSj  WRITE(mlfulv)  -  VS  Writes  the  scaled  data  to  a  file  for  future  use. 


RSft.  ;RSfp.  +  RSfs. 
1  r  1  1 


1  RS, 


—  Percent  of  total  resistance 

1  [  attributable  to  friction/wave. 


CA. - percent  of  total  resistance  that  is  attributable  to  correlation  allowance 

1  RSf 


WRITE(mlfulfp)  :  =  FF.  WRITE(mlfiilrp)  :  =  RF.  WRITE(mlfulcap)  =  CA. 
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This  Calculation  Scales  the  Test  Data  obtained  for  Model  M-2  to  the  scale  that  Model  M-1 
is  built  to.  This  enables  a  direct  comparison  of  the  Data  obtained  from  the  two  models. 


NOTE:  The  System  of  Subscripts  in  this  Routine  can  become  rather  bulky...  "R",  "C",  "M",  "S"  are  Resistance,  Coefficient 
of  Resistance,  Measured  and  Scaled  respectively...  where  measured  refers  to  the  measured  model  (Model  M-2  Data 
unsealed)  and  Scaled  refers  to  the  Model  M-2  data  scaled  to  the  Model  M-1  scale  factor.  Also,  "s",  "p",  "r“,  T,  T  as 
subscripts  are  strut,  pod,  residual,  frictional  and  total  respectively.  Other  uses  of  symbols  are  defined  as  they  arise. 

Input  Data  Section 

i  =0..  15  Sets  the  data  range 

RM  ^  j  :  =  READ (  m2 1 1  ix)  Total  measured  X  component  force  for  test  #1 :  fixed  in  pitch  and  heave, 

i 


Speed  increments  used  in  the  test,  (inputted  in  fps) 


Scaling  Factor  for  Lengths  M2/M1  scale  factor. 


Measured  Velocity  in  Knots  (converts  the  inputs  which  are  in  fps) 


Scaled  Velocity  in  Knots 


Equivalent  Scaled  Velocity  in  FPS 


^  fw  ~  1  ‘2260'  10  ^  ftA2/sec  Kinematic  Viscosity  of  Fresh  Water  at  59  F  (Test  Temp) 

0  sw  =  1-17 94  *  1 0  5  ftA2/sec  Kinematic  Viscosity  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


P  f\v  : ”  1  -9384  lb  secA2/ftA4  Density  of  Fresh  Water  at  59  F  (Test  Temp) 


p  sw  "  1  9890  lb  secA2/ftA4  Density  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


VM.  :  =  READ(m2tlv) 


16.68 


17.26 


VM  kt3,  r  VM.- 


1.688 


^kts.  VM  kts.'  V- 


VS. 

1  I  V 
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Geometry  Section 


LM  s  :=  1.4375 

SM  s  :=  1.059 

LMp  =2.32 
SMp  =2.642 

SMt  :  =  4-SMs-4  SMp 


feet  Length  of  one  Strut 


ftA2  Wetted  surface  of  one  Strut 


feet  Length  of  one  Pod 


ftA2  Wetted  Surface  of  one  Pod 


ftA2  Wetted  surface  of  the  entire  underwater  body 


Resistance  Calculations 


NOTE:  In  this  section  Rn  is  calculated  for  the  pod  and  the  strut  separately  since  their  lengths 
are  different.  Coefficients  of  friction  are  then  calculated  separately  for  each  based  on  the  1957 
ITTC  Line.  These  are  converted  to  Resistance  Values  and  then  subtracted  from  the  Total 
Resistance  value.  Finally  this  is  converted  to  a  Coefficient  of  residuary  resistance  which  is  equal 
for  the  model  and  the  ship. 


VM.-LM  s 

RnM  s  = -  operating  Rn  of  the  Strut 

j  ufw 


VM.-LM  _ 
RnM  = - - - L 

!  uf\v 


operating  Rn  of  the  Pod 


CM 


fs. 


0.075 


log/RnM  < 


Coefficient  of  Frictional  Resistance  based  on  Strut  Length 
and  using  the  1957  ITTC  line. 


Coefficient  of  Frictional  Resistance  based  on  Pod  Length 
and  using  the  1957  ITTC  line. 
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RMfs.  =CMfs.|pfv(VMi)2-4SMs 


Frictional  Resistance  of  four  Struts 


RM 


RMr 


CMr 


RnM  s 

i 

143046 

288438 

431485 

481903 

524113 

553426 

578049 

612051 

641364 

699990 

771513 

964978 

1205343 

1448053 

1688418 

1928783 


,  =CMrpi|pfw(VMi)2  4  smp 


Frictional  Resistance  of  four  Pods 


=  RM tj  -RMfs  -RMf 

i  i  *i 


Residual  Resistance  for  the  Model  as  measured. 


RM, 


jip^VUfSM, 


RnM  p 


230865 

465514 

696378 

777749 

845873 

893181 

932920 

987798 

1035106 

1129723 

1245155 

1557390 

1945318 

2337031 

2724959 

3112887 

RMtl  RMfs  RMf  RMr 

1  i  i  *1  i 


0.101 

0.025 

0.346 

0.172 

0.704 

3.5 

0.857 

5.589 

0.995 

8.37 

1.096 

9.194 

1.185 

8.126 

1.312 

5.953 

1.426 

4.619 

1.667 

3.641 

1.983 

3.425 

2.96 

2.981 

4.413 

3.843 

6.139 

5.638 

8.097 

8.51 

10.297 

11.556 

0.046 

0.156 

0.316 

0.383 

0.445 

0.49 

0.529 

0.585 

0.636 

0.742 

0.882 

1.314 

1.954 

2.713 

3.573 

4.537 

0.172 

0.673 

4.52 

6.83 

9.81 

10.78 

9.84 

7.85 

6.68 

6.05 

6.29 

7.255 

10.21 

14.49 

20.18 

26.39 

CMfs  CMf  CMr 

1  rl  1 


0.001 

0.002 

0.018 

0.023 

0.029 

0.029 

0.023 

0.015 

0.011 

0.007 

0.006 

0.003 

0.003 

0.003 

0.003 

0.003 

0.007 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.008 

0.006 

0.006 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 
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Note:  these  values  for  Rn  are  Low  for  ITTC  1957  use.  Should  be  on  order  of  10A7 


Scaled  Resistance  Calculations 


CS  r  -  CM  r  Residual  Resistance  Coefficient  remains  the  same  in  the  saling. 

i  i 


RaS  s. 

I 


VS.-LM  ,-X 
vfw 


RnS 

r  1 


VS.-LM n  X 
«  P 


uf\v 


operating  Rn  of  the  Strut  in  the  new  scale 


operating  Rn  of  the  Pod  in  the  new  scale 


CS 


0.075 


Coefficient  of  Frictional  Resistance  based  on  scaled  Strut  Length 
and  using  the  1957  ITTC  line. 


cs 


fP; 


_ 0075 _ 

(log(RnSpj  j  -  2j2 


Coefficient  of  Frictional  Resistance  based  on  scaled  Pod  Length 
and  using  the  1957  ITTC  line. 


CS„  =0 

RSfs.  =CSf4'P|\v-  VS,)24  SMs.)2 

1  I  / 

RSfPl  =  CS fp ~'Pftv  ;VS.)2.4.SMp  X2 

Rsr.  =csr.  VPftv  (VSi)2  SM  l  *-2 

RS  t  =RS  +  RS  f  -  RS  r 

1  1  r  1  1 


Correlation  Allowance 


Frictional  Resistance  of  the  Struts 


Frictional  Resistance  of  the  Pods 


Scaled  Residual  Resistance 
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Scaled  Scaled  Total  M-2 

Speed  Speed  Resistance 

(fps)  (knots)  at  17.26:1 

Scale 

VSkts.  RSt. 

1  ! 

0.711 
1.433 
2.143 
2.394 
2.603 
2.749 
2.871 
3.04 
3.186 
3.477 
3.832 
4.793 
5.987 
7.192 
8.386 
9.58 


0.16 

0.61 

4.09 

6.18 

8.87 

9.75 

8.9 

7.1 

6.05 

5.48 

5.71 

6.59 

9.28 

13.16 

18.32 

23.95 

VS. 


1.2 

2.42 

3.62 

4.04 

4.39 

4.64 

4.85 

5.13 

5.38 

5.87 

6.47 

8.09 

10.11 

12.14 

14.16 

16.17 

Residual  PodFrictional  StrutFrictional  Reynolds  #  Reynolds  # 

Component  Component  Component  of  the  Strut  of  the  Pod 


RSr. 


0.02 

0.16 

3.16 

5.04 

7.55 

00 

7.33 

5.37 

4.17 

3.29 

3.09 

2.69 

3.47 

5.09 

7.68 

10.43 

RSfp.  RSfs. 

ri  1 


0.09 

0.04 

0.32 

0.14 

0.64 

0.29 

0.78 

0.35 

0.91 

0.41 

1 

0.45 

1.08 

0.48 

1.2 

0.53 

1.3 

0.58 

1.52 

0.68 

1.81 

0.81 

2.7 

1.2 

4.02 

1.78 

5.6 

2.47 

7.38 

3.26 

9.39 

4.14 

RnSs  RnSp 

i  M 


135897 

274022 

409919 

457817 

497918 

525766 

549158 

581461 

609309 

665005 

732953 

916748 

1 145100 

1375679 

1604031 

1832383 

219326 

442248 

661574 

738877 

803596 

848540 

886293 

938428 

983372 

1073259 

1182922 

1479552 

1848092 

2220227 

2588767 

2957306 

WRITE(m2sclfx)  =  RSt_  WRITE(m2sclv)  :=  VS. 


Writes  the  scaled  data  to  a  file  for  future  use. 


i 


RS 


ft. 


FF.  =  — 
’  RS, 


RF.  =  1  -  FF.  Percent  of  total  resistance 
attributable  to  friction/wave. 


WRITE(m2ff)  :=FF. 


WRlTE(m2rf)  =RF. 
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This  Calculation  Scales  the  Test  Data  obtained  for  Model  M-2  to  Full  Scale, 


NOTE:  The  System  of  Subscripts  in  this  Routine  can  become  rather  bulky...  "R",  "C",  "M",  "S*  are  Resistance,  Coefficient 
of  Resistance,  Measured  and  Ship  respectively...  where  measured  refers  to  the  measured  model  (Model  M-2  Data 
unsealed)  and  Ship  refers  to  the  full  scale  calculation.  Also,  "s\  "p",  V,  T,  T  as  subscripts  are  strut,  pod,  residual, 
frictional  and  total  respectively.  Other  uses  of  symbols  are  defined  as  they  arise. 

Input  Data  Section 

i  ;  =  0 ..  15  Sets  the  data  range 


RM  1 1  -  READ(m2tlfx)  Total  measured  X  component  force  for  test  #1:  fixed  in  pitch  and  heave, 

i 


VM.  :  =  READ(m2tlv) 


X  ;=  16.68 


VM 


kts. 

i 


VM.- — - — 
1  1.688 


!  1 


Speed  increments  used  in  the  test,  (inputted  in  fps) 


Scaling  Factor  for  Lengths  M2/M1  scale  factor. 


Measured  Velocity  in  Knots  (converts  the  inputs  which  are  in  fps) 


Scaled  Velocity  in  Knots 


VS.  =  VM.-Jli 
1  1  v 


Equivalent  Scaled  Velocity  in  FPS 


UfSv  =1.2260- 10' 5 
Usw  =  1.1794  10  5 
PfVv  =  1-9384 
Psw  =1.9890 


ftA2/sec  Kinematic  Viscosity  of  Fresh  Water  at  59  F  (Test  Temp) 


ftA2/sec  Kinematic  Viscosity  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


lb  secA2/ftA4  Density  of  Fresh  Water  at  59  F  (T est  Temp) 


lb  secA2/ftA4  Density  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 
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Geometry  Section 


LM  s  :  =  1.4375 

feet 

Length  of  one  Strut 

SM  s  :=  1.059 

ftA2 

Wetted  surface  of  one  Strut 

LM  p  =2.32 

feet 

Length  of  one  Pod 

SMp  =2.642 

ftA2 

Wetted  Surface  of  one  Pod 

SM  j  :  =  4  •  SM  s  +  4  •  SM  p  ftA2  Wetted  surface  of  the  entire  underwater  body 


Resistance  Calculations 


NOTE:  In  this  section  Rn  is  calculated  for  the  pod  and  the  strut  separately  since  their  lengths 
are  different.  Coefficients  of  friction  are  then  calculated  separately  for  each  based  on  the  1957 
ITTC  Line.  These  are  converted  to  Resistance  Values  and  then  subtracted  from  the  Total 
Resistance  value.  Finally  this  is  converted  to  a  Coefficient  of  residuary  resistance  which  is  equal 
for  the  model  and  the  ship. 


VM.-LM  s 

RnM  „  = -  operating  Rn  of  the  Strut 

1  ufw 

VM.-LM  p 

RnM  D  = - '■ -  operating  Rn  of  the  Pod 

!  uf\v 


CM 


fs. 

I 


0.075 

log^RnM  -  2\2 


Coefficient  of  Frictional  Resistance  based  on  Strut  Length 
and  using  the  1957  ITTC  line. 


0.075 

Jlog^RnMpJ-2'i2 


Coefficient  of  Frictional  Resistance  based  on  Pod  Length 
and  using  the  1957  ITTC  line. 


215 


j  /  2 

RMfs  =  CM  fs  •—  p  fa ■  ^VM.j  -4-SMg  Frictional  Resistance  of  four  Struts 


RM^J  '--p  ^VM.j  -4-SMp  Frictional  Resistance  of  four  Pods 


RMr  :  =  RMtl  -RMfs  -RMfp 


Residual  Resistance  for  the  Model  as  measured. 


T'Pfw  (VM.j  SM  t 


143046 

288438 

431485 

481903 

524113 

553426 

578049 

612051 

641364 

699990 

771513 

964978 

1205343 

1448053 

1688418 

1928783 


230865 

465514 

696378 

777749 

845873 

893181 

932920 

987798 

1035106 

1129723 


0.101 

0.346 

0.704 

0.857 

0.995 

1.096 

1.185 

1.312 

1.426 

1.667 

1.983 

2.96 

4.413 

6.139 

8.097 

10.297 


0.025 
0.172 
3.5 
5.589 
8.37 
9.194 
8.126 
5.953 
4.619 
3.641 
3.425 
2.981 
3.843 
5.638 
8.51 
1 1.556 
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Note:  these  values  for  Rn  are  Low  for  ITTC  1957  use.  Should  be  on  order  of  10A7 


Scaled  Resistance  Calculations 


CS  r  -  CM  r  Residual  Resistance  Coefficient  remains  the  same  in  the  saling. 

i  i 

VS.-LMs-Jl 

RnS  s  - -  operating  Rn  of  the  Strut  in  the  new  scale 


RnS 


Pi 


VS.LMn-)i 
i  P 


u 


s\v 


operating  Rn  of  the  Pod  in  the  new  scale 


CS 


0.075 


fs. 


I'logjRnS  s  j  -  2j2 


Coefficient  of  Frictional  Resistance  based  on  scaled  Strut  Length 
and  using  the  1957  ITTC  line. 


CS 


fPi 


_ 0075 _ 

log/  RnS  p  J  -  l\2 


Coefficient  of  Frictional  Resistance  based  on  scaled  Pod  Length 
and  using  the  1957  ITTC  line. 


CS„  :=  0.0004 

a 


Correlation  Allowance 


RSft.  =tsfs“PSw(VSi)2-«SMsX2 

^fp.  fp.'~,f)sw'  (^i) 

I  1  z 

RS*:=cs.~Psw-(vSi)J-SMt-12 

i  2 


Frictional  Resistance  of  the  Struts 


Frictional  Resistance  of  the  Pods 


Scaled  Residual  Resistance 


Correlation 


RS  t  =RS  fs  -  RSfp  +  RS  r  +  RS  a 

i  i  *i  i  i 
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Scaled 

Speed 

(fps) 

Scaled 

Speed 

(knots) 

Total  M-2 

Resistance 

at  Full 

Scale 

Residual 

Component 

Pod  Frictional 
Component 

StrutFrictional 

Component 

Reynolds  # 

of  the  Strut 

Reynolds  # 
of  the  Pod 

VS. 

VSkts. 

RSt 

i 

RS  r 

i 

RSfPi 

RSfs. 

1 

R”Ss. 

RnSp 

2.952 

5.952 

8.904 

9.944 

10.815 

11.42 

1 1.928 

12.63 

13.235 

14.444 

15.92 

19.912 

24.872 

29.881 

34.841 

39.801 

10129794 

20425651 

30555445 

34125783 

37114902 

39190680 

40934333 

43342235 

45418013 

49569568 

54634465 

68334597 

85355972 

102543410 

119564786 

136586162 

87.1 

314.68 

659.99 

809.02 

944.46 

1044.18 

1 131.5 

1257.39 

1370.81 

1611.16 

1928.51 

2917.3 

4404.6 

6189.62 

8231.34 

10540.35 

4.98 

10.05 

15.03 

16.79 

18.26 

19.28 

20.13 

21.32 

22.34 

24.38 

26.87 

33.61 

41.98 

50.44 

58.81 

67.18 

16348607 

32965225 

49313832 

55076046 

59900225 

63250349 

66064454 

69950598 

73300723 

80000972 

88175275 

1 10286097 

137757117 

165496147 

192967168 

220438188 

200.31 

727.02 

1528.5 

1874.83 

2189.74 

2421.67 

2624.82 

2917.76 

3181.78 

3741.42 

4480.68 

6785.95 

10257.03 

14426.78 

19199.71 

24600.85 

446.65 

2025.26 

19224.96 

29761.5 

43537.84 

47855.66 

43117.98 

33268.15 

27363.12 

23664.57 

23900 

25749.59 

35850.64 

51632.94 

73622.44 

97565.13 

118.56 

818.18 

16666.37 

26616 

39857.57 

43780.95 

38697.41 

28348.32 

21992.81 

17337.95 

16307.54 

14195.24 

18300.87 

26848.18 

40524.35 

55028.49 

WRITE(m2fulfx)  :  =  RSt 


WRITE(m2fulv)  :=  VS  kns 

i 


Writes  the  scaled  data  to  a  file  for  future  use. 


RS 


RSft.  =RSfp.+  RSfs. 

I  r  1  I 


ft. 


FF. 


RS, 


RS  r 

RF.  = - -  Percent  of  total  resistance 

t  attributable  to  friction/wave. 


RSa. 


I 


percent  of  total  resistance  that  is  attributable  to  correlation  allowance 


WRITE(  m2fulfp )  -  FF.  WRITE(  m2fttlrp )  .  =  RF.  WRITE(  m2fttlcap )  .  =  C A. 


EHP 


ship. 


RS.  VS. 

1.  I 

I 

550 
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vs  las.  EHP  ship. 


4.046 

36.996 

525.35 

908.305 

1445.139 

1677.3 

1578.486 

1289.539 

1111.446 

1049.079 

1167.775 

1573.641 

2736.688 

4735.1 

7872.416 

11917.794 

2.952 

5.952 

8.904 

10.815 

11.42 

11.928 

12.63 

ME 

29.881 

34.841 

39.801 

4.983 

10.047 

15.03 

16.786 

18.256 

19.277 

20.135 

21.319 

22.34 

24.382 

26.873 

41.985 
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This  Calculation  Scales  the  Test  Data  obtained  for  Model  M-3  to  the  scale  that  Model  M-1 
is  built  to.  This  enables  a  direct  comparison  of  the  Data  obtained  from  the  two  models. 


NOTE:  The  System  of  Subscripts  in  this  Routine  can  become  rather  bulky...  “R",  “C“,  "M",  "S"  are  Resistance,  Coefficient 
of  Resistance,  Measured  and  Scaled  respectively...  where  measured  refers  to  the  measured  model  (Model  M-3  Data 
unsealed)  and  Scaled  refers  to  the  Model  M-3  data  scaled  to  the  Model  M-1  scale  factor.  Also,  "s",  "p",  Y\  T,  T  as 
subscripts  are  strut,  pod,  residual,  frictional  and  total  respectively.  Other  uses  of  symbols  are  defined  as  they  arise. 

Input  Data  Section 


i  =0..  16 


Sets  the  data  range 


RM  t|  -  READ(m3tlfx)  Total  measured  X  component  force  for  test  #1 :  fixed  in  pitch  and  heave, 
i 


VM.  -  READ(m3tlv)  Speed  increments  used  in  the  test,  (inputted  in  fps) 


15.1 

. -  Scaling  Factor  for  Lengths  M3/M1  scale  factor. 

17.26 


VM  i^c  -  VM.- -  Measured  Velocity  in  Knots  (converts  the  inputs  which  are  in  fps) 

i  1  1.688 


r 

VS  kts  =  kts  "  V ^  Scaled  Velocity  in  Knots 

i  i 


VS.  VM.-^/X  Equivalent  Scaled  Velocity  in  FPS 


.=  1.2260- 10‘ 


usw  =1. 1794-10' 
pfw  =  1.9384 
psw  =1.9890 


ftA2/sec  Kinematic  Viscosity  of  Fresh  Water  at  59  F  (Test  Temp) 


ftA2/sec  Kinematic  Viscosity  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


lb  secA2/ftA4  Density  of  Fresh  Water  at  59  F  (T est  Temp) 


lb  secA2/ftA4  Density  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 
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Geometry  Section 


LMS 

=  1.589 

SMs 

=  1.716 

LMp 

=  5.958 

SMp 

=  6.45 

SM, 

=  4SM , 

feet  Length  of  one  Strut 


ftA2  Wetted  surface  of  one  Strut 


feet  Length  of  one  Demihull 

ftA2  Wetted  Surface  of  one  Demihull 


ftA2  Wetted  surface  of  the  entire  underwater  body 


Resistance  Calculations 


NOTE:  In  this  section  Rn  is  calculated  for  the  pod  and  the  strut  separately  since  their  lengths 
are  different.  Coefficients  of  friction  are  then  calculated  separately  for  each  based  on  the  1957 
ITTC  Line.  These  are  converted  to  Resistance  Values  and  then  subtracted  from  the  Total 
Resistance  value.  Finally  this  is  converted  to  a  Coefficient  of  residuary  resistance  which  is  equal 
for  the  model  and  the  ship. 


VM.LM  s 

RjlM  r -  operating  Rn  of  the  Strut 

N  uf\v 


RnM 

n 


VM.-LM  n 
i  P 

uf\v 


operating  Rn  of  the  Demihull 


0.075 

log  (RnM  s  j  -  2)2 


Coefficient  of  Frictional  Resistance  based  on  Strut  Length 
and  using  the  1957  ITTC  line. 


CM 


fPi 


0.075 

logj'RnM  p  j  -  2)2 


Coefficient  of  Frictional  Resistance  based  on  Demihull  Length 
and  using  the  1957  ITTC  line. 
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RMfs  -  CM  j*s  ■—  p  VM.^2'4  SM  s  Frictional  Resistance  of  four  Struts 
i  i  2  l' 


RM  fp  "  CM  |p  '  ~  ‘  P  fw‘  ^ “2  •  SM  p  Frictional  Resistance  of  two  Demihulls 


RMr  .  =  RMtl  -  RMfs  -  RMf 

i  i  i  *1 


Residual  Resistance  for  the  Model  as  measured. 


CM, 


RM, 


RnM , 


RnM 


Pi 


RM 


tl. 


RMfs  RM 


fPi 


RM, 


CM 


fs. 


CM 


fPi 


CM 


334390 

451038 

502881 

540467 

593607 

619529 

644154 

671372 

773763 

896243 

1008354 

1121113 

1289604 

1401068 

1684910 

1960976 

2240931 

1253804 

1691178 

1885566 

2026498 

2225746 

2322940 

2415274 

2517328 

2901245 

3360487 

3780852 

4203646 

4835408 

5253343 

6317618 

7352736 

8402432 

0.801 

3.17 

4.22 

6.01 

7.14 

6.75 

6.35 

6.08 

6.74 

8.23 

8.93 

9.74 

12.2 

13.865 

19.08 

25.875 

33.185 

0.267 

0.453 

0.548 

0.623 

0.735 

0.793 

0.849 

0.914 

1.176 

1.527 

1.884 

2.276 

2.924 

3.391 

4.72 

6.199 

7.881 

0.372 

0.635 

0.772 

0.879 

1.041 

1.124 

1.206 

1.299 

1.678 

2.188 

2.709 

3.282 

4.231 

4.918 

6.877 

9.064 

11.56 

0.162 

2.082 

2.899 

4.508 

5.364 

4.833 

4.295 

3.867 

3.886 

4.514 

4.337 

4.182 

5.046 

5.556 

7.483 

10.612 

13.744 

0.006 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.003 

0.003 

0.003 

0.003 

0.003 

0.001 

0.009 

0.01 

0.014 

0.013 

0.011 

0.009 

0.008 

0.006 

0.005 

0.004 

0.003 

0.003 

0.002 

0.002 

0.002 

0.002 

Note:  these  values  for  Rn  are  Low  for  ITTC  1957  use.  Should  be  on  order  of  10A7 
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Scaled  Resistance  Calculations 


CSr  =  CM  r 
i  i 


Residual  Resistance  Coefficient  remains  the  same  in  the  saling. 


VS.LM  A 
RnS  = — ! - 

1  ufw 


operating  Rn  of  the  Strut  in  the  new  scale 


VS.LM  nX 
■  P 

ufw 


operating  Rn  of  the  Pod  in  the  new  scale 


0.075 

i'log^RnSs  j  -  2'; 2 


Coefficient  of  Frictional  Resistance  based  on  scaled  Strut  Length 
and  using  the  1957  ITTC  line. 


cs 


0.075 

Jlog^RnSpJ-2':2 


Coefficient  of  Frictional  Resistance  based  on  scaled  Pod  Length 
and  using  the  1957  ITTC  line. 


csa  =0 


Correlation  Allowance 


RS 


fs. 


=  csfs;--pfW-fVS.->2-4-SM  A2 


i  2 


Frictional  Resistance  of  the  Struts 


Frictional  Resistance  of  the  Demi  Hulls 


RS  r  =  CS  r ,  “  Pfvv'  ^VS.J2  SM  t-A2  Scaled  Residual  Resistance 

RS  t.  =  RS  fs.  +  RS  fp.  +  RS  r 

i  l  ri  i 
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Scaled 

Scaled 

Total  M-3 

Speed 

Speed 

Resistance 

(fps) 

(knots) 

at  17.26:1 

Scale 

vs. 

VSkts. 

1 

RS  t. 
i 

0.56 

2.16 

2.87 

4.07 

4.83 

4.58 

4.31 

4.14 

4.59 

5.62 

6.11 

6.67 

8.36 

9.51 

13.08 

17.72 

22.72 

2.41 

3.25 

3.63 

3.9 

4.28 

4.47 

4.65 

4.85 

5.58 

6.47 

7.28 

8.09 

9.31 

10.11 

12.16 

14.15 

16.17 

1.43 

1.928 

2.15 

2.311 

2.538 

2.649 

2.754 

2.87 

3.308 

3.832 

4.311 

4.793 

5.513 

5.99 

7.203 

8.384 

9.581 

Residual  PodFrictional  StrutFrictional  Reynolds  #  Reynolds  # 

Component  Component  Component  of  the  Strut  of  the  Pod 


RSr. 

1 

oil" 

1.39 

1.94 

3.02 

3.59 
3.24 
2.88 

2.59 
2.6 
3.02 
2.9 
2.8 
3.38 
3.72 
5.01 
7.11 
9.2 


RSfps 

0.26 

0.44 

0.54 

0.61 

0.73 

0.78 

0.84 

0.91 

1.17 

1.52 

1.88 

2.28 

2.94 

3.42 

4.78 

6.29 

8.02 


RSfSj 


0.19 

0.32 

0.38 

0.44 

0.52 

0.56 

0.6 

0.64 

0.82 

1.07 

1.32 

1.59 

2.04 

2.37 

3.29 

4.32 

5.49 

RnS  s  RnS  p 

i  *i 


369077 

411500 

442256 

485739 

506951 

527101 

549373 

633158 

733381 

825120 

917390 

1055263 

1 146472 

1378736 

1604636 

1833718 

1025968 

1383864 

1542929 

1658251 

1821293 

1900825 

1976381 

2059890 

2374043 

2749834 

3093812 

3439778 

3956738 

4298728 

5169608 

6016628 

6875579 

WRITE(m3sclfx)  =RSt  WRITE(m3sclv)  ■=  VS 


Writes  the  scaled  data  to  a  file  for  future  use. 


RS  ft.  RS*p.  ;  RS  fs 


i 


RF.  =  1  -  FF. 

i  i 


Percent  of  total  resistance 
attributable  to  friction/wave. 


WRITE(m3ff)  =FF. 


WRITE(m3rf)  :=RF. 
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This  Calculation  Scales  the  Test  Data  obtained  for  Model  M-3  to  full  scale. 


NOTE:  The  System  of  Subscripts  in  this  Routine  can  become  rather  bulky...  "R“,  "C",  "M",  "S"  are  Resistance,  Coefficient 
of  Resistance,  Measured  and  Scaled  respectively...  where  measured  refers  to  the  measured  model  (Model  M-3  Data 
unsealed)  and  Scaled  refers  to  the  Model  M-3  data  scaled  to  the  Model  M-1  scale  factor.  Also,  "s",  "p",  V,  T,  "t"  as 
subscripts  are  strut,  pod,  residual,  frictional  and  total  respectively.  Other  uses  of  symbols  are  defined  as  they  arise. 

Input  Data  Section 

i  -0..  16  Sets  the  data  range 

RM  •  -  READ(  m3tlfx)  Total  measured  X  component  force  for  test  #1:  fixed  in  pitch  and  heave, 
i 


VM.  -  READ(m3tl  v)  Speed  increments  used  in  the  test,  (inputted  in  fps) 


X  15.1  Scaling  Factor  for  Lengths  M3/M1  scale  factor. 


VM  t,tc  -  VM. - Measured  Velocity  in  Knots  (converts  the  inputs  which  are  in  fps) 

>  1  1 .688 

VS  ktS  =  VM  •  V  X  Scaled  Velocity  in  Knots 

i  i 

r 

VS.  r  VM. A  Equivalent  Scaled  Velocity  in  FPS 


l)  ^  -  1 .2260- 1  0  J  ftA2/sec  Kinematic  Viscosity  of  Fresh  Water  at  59  F  (Test  Temp) 

u  SW  =  1  •  1 794-  1 0  5  ftA2/sec  Kinematic  Viscosity  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 


P  f\v  =  1  9384  lb  secA2/ftA4  Density  of  Fresh  Water  at  59  F  (Test  Temp) 


p  sw  :=  1 .9890  lb  secA2/ftA4  Density  of  Salt  Water  at  65  F  (Assumed  Operating  Temp) 
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Geometry  Section 


LMS 

.=  1.589 

feet 

Length  of  one  Strut 

SMs 

-  1.716 

ftA2 

Wetted  surface  of  one  Strut 

LMp 

=  5.958 

feet 

Length  of  one  Demihull 

SMp 

=  6.45 

ftA2 

Wetted  Surface  of  one  Demihull 

SMt 

=  4-SM$-  2  SM p 

ftA2 

Wetted  surface  of  the  entire  underwater  body 

Resistance  Calculations 


NOTE:  In  this  section  Rn  is  calculated  for  the  pod  and  the  strut  separately  since  their  lengths 
are  different.  Coefficients  of  friction  are  then  calculated  separately  for  each  based  on  the  1 957 
ITTC  Line.  These  are  converted  to  Resistance  Values  and  then  subtracted  from  the  Total 
Resistance  value.  Finally  this  is  converted  to  a  Coefficient  of  residuary  resistance  which  is  equal 
for  the  model  and  the  ship. 


VM.LM  s 

RnM  s  = -  operating  Rn  of  the  Strut 

!  uf\v 

VM.LM 

RnM  _  - - - -  operating  Rn  of  the  Demihull 


CM 


fs. 


0.075 

I  log /RnM  s  j  -  2 

'•  \  i  / 


Coefficient  of  Frictional  Resistance  based  on  Strut  Length 
and  using  the  1957  ITTC  line. 


CM 


fPi 


0.075 

log:  RnM  p  .  -  2' 2 


Coefficient  of  Frictional  Resistance  based  on  Demihull  Length 
and  using  the  1957  ITTC  line. 
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Frictional  Resistance  of  four  Struts 


RM 


RM 


RM  r 


CM  r 


RnM  s 

i 

334390 

451038 

502881 

540467 

593607 

619529 

644154 

671372 

773763 

896243 

1008354 

1121113 

1289604 

1401068 

1684910 

1960976 

2240931 


=  CM 


fs.|pfv(VMi}’4SM 

i  2 


s 


:  =  CM 


•2-SM , 


Frictional  Resistance  of  two  Demihulls 


=  RMtl  -  RMfs  -  RMfp 

i  i  *i 


Residual  Resistance  for  the  Model  as  measured. 


RM, 


RnM  D 

1  i 


1253804 

1691178 

1885566 

2026498 

2225746 

2322940 

2415274 

2517328 

2901245 

3360487 

3780852 

4203646 

4835408 

5253343 

6317618 

7352736 

8402432 

RM  ^  j  RM 

i  i 


0.801 

0.267 

3.17 

0.453 

4.22 

0.548 

6.01 

0.623 

7.14 

0.735 

6.75 

0.793 

6.35 

0.849 

6.08 

0.914 

6.74 

1.176 

8.23 

1.527 

8.93 

1.884 

9.74 

2.276 

12.2 

2.924 

13.865 

3.391 

19.08 

4.72 

25.875 

6.199 

33.185 

7.881 

RMf  RMr 

i 


0.162 

2.082 

2.899 

4.508 

5.364 

4.833 

4.295 

3.867 

3.886 

4.514 

4.337 

4.182 

5.046 

5.556 

7.483 

10.612 

13.744 

0.372 

0.635 

0.772 

0.879 

1.041 

1.124 

1.206 

1.299 

1.678 

2.188 

2.709 

3.282 

4.231 

4.918 

6.877 

9.064 

11.56 

CMfs  CMf  CMr 

1  ri  1 


0.001 

0.009 

0.01 

0.014 

0.013 

0.011 

0.009 

0.008 

0.006 

0.005 

0.004 

0.003 

0.003 

0.002 

0.002 

0.002 

0.002 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.004 

0.003 

0.003 

0.003 

0.003 

0.003 

0.006 

0.006 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.005 

0.004 

0.004 

0.004 

0.004 

0.004 

Note:  these  values  for  Rn  are  Low  for  ITTC  1957  use.  Should  be  on  order  of  10A7 


Scaled  Resistance  Calculations 


CS  r  -  CM  r  Residual  Resistance  Coefficient  remains  the  same  in  the  saling. 

i  i 

VS.  LM  s‘l 

RnS  s  - -  operating  Rn  of  the  Strut  in  the  new  scale 


VS.LMnX 

RnSp 

*  r 

*  usw 

operating  Rn  of  the  Pod  in  the  new  scale 

CSfs. 

1 

0.075 

(log(RnSs.j-2)2 

Coefficient  of  Frictional  Resistance  based  on  scaled  Strut  Length 
and  using  the  1957  ITTC  line. 

CSfp 

0.075 

Coefficient  of  Frictional  Resistance  based  on  scaled  Pod  Length 

Ipi 

(log(RnSpji-2j2 

and  using  the  1957  ITTC  line. 

CSa 

=  0.0004 

Correlation  Allowance 

RSfs. 

1 

i  2 

Frictional  Resistance  of  the  Struts 

RSfPj 

=  CSlp.PPSw(VS022SMPi2 

Frictional  Resistance  of  the  Demi  Hulls 

RS  r 

i 

=  csrPPsw(VSi)2'SMt‘2 

Scaled  Residual  Resistance 

RS  a. 

i 

Correlation 

Rst. 

i 

=  RSfs+RSf  +RSr  +  RSa 

1  r  1  1  1 
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Scaled 

Scaled 

Total  M-3 

Speed 

Speed 

Resistance 

(fps) 

(knots) 

at  Full 

Scale 

vs. 

VSkts. 

1 

RS  t 

i 

1805.61 

9515.06 

12890.08 

18956.26 

22563.11 

20988.62 

19383.41 

18211.4 

19668.69 

23779.04 

25097.62 

26706.78 

33361.84 

37792.85 

52161.65 

71750.23 

92642.78 

10.03 

13.52 

15.08 

16.2 

17.8 

18.57 

19.31 

20.13 

23.2 

26.87 

30.23 

33.61 

38.66 

42.01 

50.52 

58.79 

67.19 

5.939 

8.011 

8.932 

9.6 

10.543 

11.004 

11.441 

11.925 

13.743 

15.919 

17.91 

19.913 

22.905 

24.885 

29.927 

34.83 

39.803 

Residual  PodFrictional  StrutFrictional  Reynolds  #  Reynolds  # 

Component  Component  Component  of  the  Strut  of  the  Pod 


RS  r.  RSfp.  RSfs. 

1  1 


416.2 

721.46 

881.47 

1006.64 

1196.59 

1294.71 

1391.23 

1501.61 

1951.46 

2560.23 

3183.76 

3873.76 

5020.44 

5854.09 

8242.17 

10923.55 

13996.93 

637.02 

1109.41 

1357.69 

1552.15 

1847.58 

2000.34 

2150.67 

2322.7 

3024.61 

3976.3 

4952.65 

6034.61 

7835.39 

9146.31 

12907.73 

17139.31 

21997.35 

572.22 

7356.37 

10243.43 

15926.77 

18951.13 

17075.09 

15172.89 

13660.76 

13727.86 

15948.15 

15322.77 

14773.05 

17826.12 

19629.28 

26437.11 

37490.84 

48556.41 

RnSs 

I 


20396136 

27511067 

30673259 

32965848 

36207094 

37788190 

39290231 

40950382 

47195710 

54666388 

61504627 

68382394 

78659517 

85458229 

102771228 

1 19609899 

136685734 

RnS 

r  1 


76475883 

10315351 

11501024 

12360637 

13575951 

14168787 

14731982 

15354460 

17696163 

20497315 

23061332 

25640170 

29493606 

32042802 

38534359 

44848066 

51250698 

WRITE(m3fulix)  :=RSt 

i 

WRITE(m3fi)lv)  :=VSkls 

i 

Writes  the  scaled  data  to  a  file  for  future  use. 

Rs  ft.  =RSfp  .  +  Rsfs. 

i  r  i  i 

RS  ft. 

FF  -  1  PF 

RS  r. 

1 

Percent  of  total  resistance 

attributable  to  friction/wave. 

II.-  X  VI  . 

1  RS , 

Rs  t. 

1  I 


1 


percent  of  total  resistance  that  is  attributable  to  correlation  allowance 
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WRITE(m3fulfp)  =FF.  WRITE(m3fulrp)  -RF. 


EHP 


ship. 


RS.VS. 
1.  1 
l 

550 


VS. 


VS 


kts. 


EHP 


ship. 


10.026 

13.523 

15.077 

16.204 

17.797 

18.574 

19.313 

20.129 

23.199 

26.871 

30.232 

33.613 

38.664 

42.006 

50.516 

58.793 

67.187 

5.939 

8.011 

8.932 

9.6 

10.543 

11.004 

11.441 

11.925 

13.743 

15.919 

17.91 

19.913 

22.905 

24.885 

29.927 

34.83 

39.803 

32.913 

233.946 

353.357 

558.489 

730.113 

708.823 

680.632 

666.498 

829.613 

1 161.749 

1379.551 

1632.163 

2345.303 

2886.43 

4790.937 

7669.871 

11317.027 

WRITE(m3fulcap)  :=CA. 


M-3  in  full  scale 
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Appendix  H 


Approximation  of  Impingement  Force 
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This  Worksheet  approximates  the  impingement  force  created  by 
the  rooster  tail  on  the  towing  frame. 

,,  lb 

p  =o4— -  density  of  water 

ft3 


V 


XT 


sec 


estimate  of  the  maximum  relative  x  component 
of  velocity  between  the  frame  and  water  before  impact. 


xa 


sec 


water  x  component  of  velocity  after  the  impact. 


an  approximation  for  the  time  of  impact 
of  the  water  on  the  frame 


acceleration  of  the  water 


a  x  =40*  ft*  sec 


vol  -  1 .25- in-6- in-  ^  t ;  an  approximation  for  the  volume  of  water  that  is 

impacted  during  the  impact  period. 


vol  =0.13*  ft3 


m  w  =  pvol 


m  w  =  8.333*lb 


mass  of  the  water  that  is  impacted  during  the  impact 
period. 


F  impact  mw'ax 


Approximated  Impact  Force  on  the  frame. 


F  impact  =  '0  36 -lbf 
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Appendix  I 


Numerical  Codes 
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********** 


Q  m********** 


%  This  routine  solves  for  the  total  resistance  of  an  8  bodied  SLICE  type 
%  vessel.  Solution  to  the  friction  resistance  is  based  on  the  ITTC  1957 
%  friction  line.  Solution  to  the  wave  resistance  portion  is  based  on  the 
%  solution  to  Mitchell's  Integral  proposed  bv  Lunde  in  1951  and  expanded 
%  by  Wilson. 


%  GENERAL  INPUTS  SECTION 

ls=  input('Input  the  slowest  integer  speed  as  x.xx  (f/s)  ='): 
hs=input('The  number  of  increments- ); 
lani=input('The  incrementation  (decimal)^); 

g=32.2; 


water=input('For  Fresh  Water  enter  L  Salt  Water  enter  2:  '); 
if  water==l ; 
rho=  1.9367; 
nu=  1.0804E-5; 
el  seif  water==2; 
rho=  1 .9912: 
nu=  1.3343E-5: 
else 

disp('  ') 

disp('  ') 

disp('*  *******************************  **********') 
disp('A  valid  entry  was  not  made  for  water  type') 
disp('The  program  may  not  provide  valid  results’) 
disp('  ******************************************') 
end 


end 


disp('  ') 
disp('  ’) 
disp('  ') 

disp('The  body  numbering  is:  l=fwd  stbd  strut') 
disp('  2=fwd  port  strut') 
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disp(' 

disp(' 

disp(' 

disp(' 

dispC 

dispC 

dispC  ’) 

dispC  ’) 

dispC  ’) 

disp('The 

dispC 

dispC 

dispC 

dispC 


3=aft  stbd  strut') 
4=aft  port  strut') 
5=fwd  stbd  pod') 
6=fwd  port  pod') 
7=aft  stbd  pod’) 
8=aft  port  pod') 


global  coordinate  system  is:') 

x  positive  in  the  direction  of  motion') 
y  positive  to  port') 
z  positive  upward') 

The  origin  is  at  the  undisturbed  free  surface') 


TS=  inputCThe  height  of  the  struts  (ft)-); 

LS=  inputCThe  length  of  the  struts  (ft)='): 

0/oWSS=mput('The  wetted  Surface  of  one  strut  (ftA2)-); 

cs6=input('The  6th  power  coefficient  for  the  strut  potential  fn='); 
cs5=input('The  5th  power  coefficient  for  the  strut  potential  fn='); 
cs4=input('The  4th  power  coefficient  for  the  strut  potential  fn='); 
cs3=input('The  3rd  power  coefficient  for  the  strut  potential  fn— ); 
cs2=input('The  2nd  power  coefficient  for  the  strut  potential  fn- ); 
csl=input('The  1st  power  coefficient  for  the  strut  potential  fn-); 
csO=input('The  Oth  power  coefficient  for  the  strut  potential  fn='); 

TP=  input('The  max  diameter  of  the  pods  ( ft)- ); 

LP=  input('The  length  of  the  pods  (ft)-); 

%WSP=input('The  wetted  surface  of  one  pod  (ftA2)- ); 

cp6=inputCThe  6th  power  coefficient  for  the  pod  potential  fn- ); 
cp5=input('The  5th  power  coefficient  for  the  pod  potential  fn-); 
cp4=input('The  4th  power  coefficient  for  the  pod  potential  fn-); 
cp3=input('The  3rd  power  coefficient  for  the  pod  potential  fn- ); 
cp2=input('The  2nd  power  coefficient  for  the  pod  potential  fn- ); 
cpl=input('The  1st  power  coefficient  for  the  pod  potential  fn-); 
cpO=input('The  Oth  power  coefficient  for  the  pod  potential  fn- ); 


sbl=input('The  x  setback  of  body  1  (+  or  -  ft)-); 
osl=input('The  y  offset  of  body  1  (+  or  -  ft)-); 
subl=input('The  z  submergence  of  bodyl  (-  or  -  ft)-); 
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sb2=input(’The  x  setback  of  body  2  (-  or  -  ft)=’); 
os2=input('The  y  offset  of  body  2  (+  or  -  ft)='); 
sub2=input('The  z  submergence  of  body  2  (+  or  -  ft)-); 


sb3=input('The  x  setback  of  body  3  (-  or  -  ft)-  ); 
os3=input('The  y  offset  of  body  3  (+  or  -  ft)='); 
sub3=input('The  z  submergence  of  body 3  (+  or  -  ft)='); 


sb4=input(’The  x  setback  of  body  4  (-  or  -  ft)-); 
os4=input('The  v  offset  of  body  4  (~  or  -  ft)='); 
sub4=input('The  z  submergence  of  body4  (+  or  -  ft)='); 


sb5=input('The  x  setback  of  body  5  (-  or  -  ft)=’); 
os5=input('The  v  offset  of  body  5  (-  or  -  ft)='); 
sub5=input('The  z  submergence  of  bodv5  (+  or  -  ft)='); 


sb6=input(’The  x  setback  of  body  6  (-  or  -  ft)-); 
os6=input('The  y  offset  of  body  6  (-r  or  -  ft)='); 
sub6=input('The  z  submergence  of  bodv6  (+  or  -  ft)='); 


sb7=input('The  x  setback  of  body  7  (-  or  -  ft)=’); 
os7=input('The  v  offset  of  body  7  (+  or  -  ft)=’); 
sub7=input(’The  z  submergence  of  body 7  (+  or  -  ft)='); 


sb8=input('The  x  setback  of  body  8  (-  or  -  ft)='); 
os8=input('The  y  offset  of  body  8  (+  or  -  ft)='); 
sub8=input('The  z  submergence  of  bodv8  (+  or  -  ft)='); 


tseg=input('The  number  of  Theta  Segments- ); 
zseg=  input('The  number  of  segments  in  Z- ); 
xseg=  input('The  number  of  segments  in  X— ); 
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°o  This  section  calls  the  routines  that  computes  the  Rw  value 
%  for  each  of  the  body  and  interface  terms  and  the  Rf  value  for 
%  the  body  terms  only. 

Rwb=0: 

Rwl  1=0: 

Rw22=0: 

Rw33=0: 

Rw44=0: 

Rw55=0: 

Rw66=0: 

Rw77=0: 

Rw88=0: 

Rws=0: 

Rwl  5=0: 

Rw26=0: 

Rw37=0: 

Rw48=0: 

R\vl=0; 

Rwl  3=0: 

Rwl  7=0; 

Rw53=0: 

Rw57=0: 

Rw24=0: 

Rw28=0: 

Rw64=0: 

Rw68=0: 

Rwd=0: 

Rwl4=0; 

Rwl  8=0: 

Rw54=0: 

Rw58=0: 

Rw23=0: 

Rw27=0: 

Rw63=0; 

Rw67=0: 

Rwt=0; 

Rwl  2=0: 

Rwl  6=0: 

Rw52=0; 

Rw56=0: 

Rw34=0; 

Rw38=0: 

Rw74=0: 
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Rw78=0; 

Rf=0; 

Rfs=0; 

Rfp=0; 

Rt=0; 


o/0***********  body  TERMS********* 

oneone 

twotwo 

thrthr 

forfor 

fivfiv 

sixsix 

sevsev 

ateate 

%*******SAME  CORNER  INTERACTION  TERMS****** 

onefiv 

twosix 

thrsev 

forate 

%*****  *  LONGITUDINAL  INTERACTION  TERMS********* 

onethr 

onesev 

fivthr 

fivsev 

twofor 

twoate 

sixfor 

sixate 

%*****DIAGONAL  INTERACTION  TERMS********* 

onefor 

oneate 

fivfor 

fivate 

twothr 

twosev 

sixthr 

sixsev 

%******TRANSVERSE  DIAGONAL  TERMS*********** 
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onetwo 

onesix 

flvtwo 

fivsix 

thrfor 

thrate 

sevfor 

sevate 


Rwb=Rwl  1+Rw22-Rvv33-Rw44+Rw55+Rw66+Rw77+Rw88; 
Rws=Rwl  5-rRw26^R\v37-rRw48; 

Rwl=Rw]3-t-Rwl7-Rw53-Rw57+Rw24+Rw28+Rw64+Rw68; 

Rwd=Rwl4+Rwl8-Rw54-Rw58+Rw23+Rw27+Rw63+Rw67; 

Rwt=Rwl2-LRwl6-Rw52^Rw56+Rw34+Rw38+Rw74+Rw78; 


Rw-Rwb^-Rws+Rwl-Rwd^Rwt 


Rf=Rfs+Rfp 

Rt=Rw-i-Rf 


plot(V,Rt,'-\V,Rw/-.') 
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%  *******************podgeo************************** 

% 

%  This  file  creates  the  offsets  for  the  pod  by  obtaining  nose 
%  tail  and  parallel  midbody  information.  The  nose  is  formed 
%  using  an  ellipse  and  the  tail  is  formed  using  a  parabola. 

%  This  is  then  fit  to  a  seventh  order  polynomial  and  the  results  compared. 

seg=input('The  min  number  of  segments  (must  be  even)  for  the  pod  =  '); 
D=input(’The  max.  diameter(ft)  of  the  pod  =  '); 

L2D=input('The  desired  Legnth  to  Diameter  ratio  of  the  pod  =  '); 
Lpmb=input('  The  amount  of  Parallel  Midbody(ft)  in  the  pod  =  '); 

Pn=input('  The  decimal  amount  of  the  pod  that  is  nose  section  =  '); 

Lpod=  D*L2D; 

Ln=  Pn*(Lpod-Lpmb); 

Lt=  Lpod-Ln-Lpmb; 


npn— input)'  The  nose  shape  exponent. (bigger  is  fuller)  =  '); 
npt=input('  The  tail  shape  exponent. (bigger  is  fuller)  =  '); 

ppod=zeros(  1 ,8); 
xn=zeros(  l,seg*2); 
offn=zeros(  1  ,seg*2); 
xpmb=zeros(  1  ,seg/2); 
offpmb=zeros(  1  ,seg/2); 
xt=zeros(l,seg*2); 
offt=zeros(  1  ,seg*2); 
xpod=zeros(  1  ,(seg*4+seg/2)); 
offpod=zeros(  1  ,(seg*4+seg/2)); 

Xpod=zeros(  1  ,seg*6); 

OFFpod=zeros(  1  ,seg*6); 

xn=  (linspace(0,-Ln,seg*2))’; 
offn=  D/2*(l-((Ln-abs(xn))./Ln)/npn).A(l/npn); 


xpmb=  (linspace(-Ln,-(Ln+Lpmb),seg/2))'; 
offpmb=  D/2*xpmb./xpmb; 


xt  =  (linspace(-(Ln+Lpmb),-Lpod,seg*2))'; 
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offt=  D/2*(l-((abs(xt)-Ln-Lpmb)./Lt)  Anpt): 


xpod=[xn:xpmb;xt]: 

offpod=[offn;offpmb;offt]; 


ppod=polyfit(xpod,offpod,7)' 


Xpod=(  linspace(0,-Lpod,seg*6)); 

OFFpod=(ppod(l).*Xpod  A7)+(ppod(2).*Xpod  A6)+(ppod(3).*Xpod  A5)+... 
(ppod(4).*Xpod  A4)+(ppod(5).*Xpod.  3)+(ppod(6).*Xpod  A2)+... 
(ppod(7)*Xpod.Al)-t-ppod(8); 


plot(Xpod.OFFpod,'-^Xpod.-OFFpod,,:,) 
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%  ***************  ****£trutgeo*  **************  *********** 

% 

%  This  file  creates  the  offsets  for  the  strut  by  obtaining  nose 
%  tail  and  parallel  midbodv  information.  The  nose  and  tail  are  formed 
%  using  a  parabola.  This  is  then  fit  to  a  seventh  order  polynomial 
%  and  the  results  compared. 


segs  input(  The  min  number  of  segments  (must  be  even)  for  the  strut  =  '); 
B=input('The  max.  beam(ft)  of  the  strut  =  '); 

L2Ds=input('The  desired  Legnth  to  Diameter  ratio  of  the  strut  =  '); 
Lpmbs=input('  The  amount  of  Parallel  Midbody(ft)  in  the  strut  =  ’); 
Pns=inputf  The  decimal  amount  of  the  strut  that  is  nose  section  =  ’); 

Lstrut=  B*L2Ds; 

Lns=  Pns*(Lstrut-Lpmbs); 

Lts=  Lstrut-Lns-Lpmbs; 


nsn=inputf  The  nose  shape  exponent,(bigger  is  fuller)  =  '); 
nst=inputf  The  tail  shape  exponent,(bigger  is  fuller)  =  '); 


xns=zeros(  1  ,segs*2); 
offhs=zeros(  1 , segs*  2); 
xpmbs=zeros(l  ,segs/2); 
offpmbs=zeros(  1  ,segs/2); 
xts=zeros(l,segs*2); 
offts=zeros(  1  ,segs*2); 
xstrut=zeros(  1  ,(segs*4+segs/2)); 
offstrut=zeros(  1  ,(segs*4+segs/2)); 
Xstrut=zeros(  1  ,segs*6); 

OFFstrut=zeros(  1  ,segs*6); 

xns=  (linspace(0,-Lns,segs*2))'; 
offns=  B/2*(  ]-((Lns-abs(xns))./Lns).Ansn); 


xpmbs=  (linspace(-Lns,-(Lns+Lpmbs),segs/2))'; 
offpmbs=  B/2*xpmbs./xpmbs; 


xts  =  (linspace(-(Lns+Lpmbs),-Lstrut,segs*2))'; 
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offts=  B/2*(l-((abs(xts)-Lns-Lpmbs)./Lts).Anst); 


xstrut=[xns;xpmbs:xts]; 

offstrut=[offns;offpmbs:offts]; 


ps=polyfit(xstrut,offstrut.7)' 


Xstrut=(linspace(0,-Lstrut,segs*6))'; 

OF  F  strut=(  ps(  1 ) .  *  X  strut . '  7 )+( ps(  2 ) .  *  X  strut .  A6)+(ps(  3 ) .  *  X  strut  A5 )+ . . . 
(ps(4).*Xstrut  A4)-(ps(5).*Xstrut  A3)+(ps(6).*Xstrut  A2V-... 

(ps(  7 ) .  *  Xstrut  )-(ps(  8)); 


plot(Xstrut,OFFstrut.'--',xstrut,offstrut.'-'); 


%*  *****************  ONEONE 


theta  1  =zeros(  tseg+ 1,1); 
z  1  =zeros(  zseg+ 1,1); 
xl=zeros(xseg+l,l); 
p  1  x=zeros(  xseg+ 1,1); 
qlx=zeros(xseg+l,l); 
p  1  z=zeros(  zseg+ 1 , 1 ) ; 
qlz=zeros(zseg+l,l); 
plt=zeros(tseg+l,l); 
qlt=zeros(tseg-i-l,l); 
resl  l=zeros(tseg+l,l); 
q  1  At=zeros(  tseg+ 1,1); 
p  1  At=zeros(  tseg+ 1,1); 
p  1  Az=zeros(  zseg+ 1,1); 
q  1  Az=zeros<  zseg+ 1,1); 
zl  A=zeros(zseg+l,l); 
p  1  Ax=zeros(  xseg+ 1,1); 
q  1  Ax=zeros(  xseg+ 1,1); 
x  1  A=zeros(  xseg+ 1 , 1 ); 
sara  1  =zeros(  xseg+ 1,1); 
katep  1  =zeros(xseg+ 1,1); 
kateq  1  =zeros(  xseg+ 1,1); 
greg  1  =zero  s(  xseg+ 1 , 1 ); 
sara  1  A=zer  o  s(  x  seg+ 1 , 1 ) ; 
katep  1  A=zeros(xseg+ 1,1); 
kateq  1  A-zeros(xseg+ 1,1); 
greg  1  A=zeros(xseg+ 1,1); 
V=zeros(hs- 1  -Is,  1 ); 

Rwl  l=zeros(hs+l-ls,l); 


for  iv=l  :hs; 
V(iv)=ls-(iv*lani); 
k=g/(V(ivr2); 
w=-V(iv)/(2*pi); 


for  it=l:tseg+l ; 
theta  1  (it)=(pi/2)*((it- 1  )/tseg); 


for  izl  =  1  :zseg-*-l ; 
zl(izl)  =  -TS*((izl-l  )/zseg); 

for  ixl=l:xseg+l; 
xl(ixl)  =-LS*((ixl-l)/xseg): 


sara  1  ( i  x  1  )=w*  (( cs6*  x  1  ( ix  1  )A6)+(cs5  *  x  1  (ix  1 )''  5 )+( cs4  *  x  1  ( ix  1  )A4)+(cs3  *  x  1  (ix  1  )A3  )+(cs2 
*x  1  (ix  1  )A2)+(cs  1  *x  1  (ix  1  ))+csO); 

katepl(ixl)=cos(k*(sbl*cos(thetal(it)hosl*sin(thetal(it)))*(sec(thetal(it)))A2); 
kateq  1  (ixl  )=sin(k*(sbl  *cos(  theta  1  (it))— osl  *sin(thetal  (it)))*(sec(thetal  (it)))A2); 
gregKixl  )=exp(k*subl*(sec(thetal(it)))A2); 

plx(ixl  )=saral(ixl  )*katepl(ix1  )*gregl(ixl ): 
q  1  x(  ix  1  )=saral  (ix  1  )*kateq  1  (ix  1  )*greg  l(ix  1 ); 

end  0 othis  end  closes  the  ****  for  ixl=l:xseg-l  loop  **** 

p  1  z(  iz  1  )=sum(p  1  x)/'xseg; 
qlz(izl)=sum(qlx)/xseg; 

end  0 othis  end  closes  the  ****  for  izl  =  l:zseg-l  loop  **** 

plt(it)=sum(plz)/zseg; 
q  1 1(  it)=sum(q  1  z)/zseg; 

%  *******This  starts  the  next  body  calculations  inside  the  theta  loop***** 


for  izl  A=  1  :zseg-l; 
zlA(izlA)  =  -TS*((izlA-l)/zseg); 

for  ixl A=1  :xseg+l; 
xlA(ixlA)  =-LS*((ixl  A-l  )/xseg); 

saral  A(  ix  1  A)=w*((cs6*x  1  A(ix  1  A)A6)+(  cs5*x  1  A(ix  1  A)A5)-(  cs4*x  1  A(ix  1  A)A4)+(cs3*x 
1  A(ix  1  A)A3)+(cs2*x  1  A(ix  1  A)A2)+(csl  *x  1  A(ix  1  A))+csO); 

katepl  A(ixl  A)=cos(k*(sbl  *cos(thetal(it))+osl*sin(thetal(it)))*(sec(thetal(it)))A2); 
kateql  A(ixl  A)=sin(k*(sbl  *cos(thetal(it))+osl  *sin(thetal(it)))*(sec(thetal(it)))A2); 
gregl  A(  ix  1  A)=exp(k*subl  *(sec(thetal(it)))A2); 

p  1  Ax(ix  1  A)=sara  1  A(ix  1  A)*katep  1  A(  ix  1  A)*greg  1  A(ix  1  A); 
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qlAx(ixlA)=saralA(ixlA)*kateqlA(ixlA)*greglA(ixlA); 

end  %this  end  closes  the  ****  for  ixlA=l:xseg+l  loop  **** 

p  1  Az(iz  1  A)=sum(p  1  Ax)/xseg; 
q  1  Az(iz  1  A)=sum(q  1  Ax)/xseg; 

end  %this  end  closes  the  ****  for  izlA=l:zseg+l  loop  **** 

p  1  At(it)=sum(pl  Az)/zseg; 
q  1  At(it)=sum(q  1  Az)/zseg; 

resl  l(it)=16*pi*rho*kA2*((plt(it)*pl  At(it))+(qlt(it)*qlAt(it)))*(sec(thetal(it))A3); 

end  %this  end  closes  the  ****  for  it=  l:tseg+l  loop  **** 

Rwll(iv)=  sum(resll)/tseg 
end  %  this  closes  the  ****  for  iv  =  ...  loop  ***** 
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%*  *****  **  **  *  **********  SIXSIX 


Rw66=zeros(hs~  1  -1  s,  1 ); 
res66=zeros(tseg^  1,1); 
theta  1  =zeros(tseg- 1,1); 
zl=zeros(zseg+l,l); 
xl=zeros(xseg-i-l,l); 
p  1  x=zeros(xseg^- 1.1); 
q  1  x=zeros(xseg- 1,1); 
plz=zeros(zseg*l,l); 
qlz=zeros(zseg-!-l,l); 
plt=zeros(tseg+l,l); 
qlt=zeros(tseg-*-l.l); 
q  1  At=zeros(tseg^  1,1); 
p  1  At=zeros(tseg- 1,1); 
p  1  Az=zeros(zseg- 1,1); 
q  1  Az=zeros(zseg- 1.1): 
z  1  A=zeros(zseg-^  1.1); 
p  1  Ax=zeros(xseg- 1,1); 
q  1  Ax=zeros(xseg- 1,1); 
x  1  A=zeros(xseg- 1,1); 
sara  1  =zeros(xseg- 1,1); 
katep  1  =zeros(xseg- 1,1); 
kateq  1  =zeros(xseg- 1,1): 
greg  1  =zeros(xseg- 1,1); 
saral  A=zeros(xseg-l  ,1 ); 
katep  1  A=zeros(xsegj- 1,1); 
kateq  1  A=zeros(xseg^l  ,1 ); 
greg  1  A=zeros(xseg- 1,1 ); 
V=zeros(hs+ 1  -Is.  1 ): 


for  iv=l  :hs; 
V(iv)=ls+(iv*lani); 

k=g/(V(iv)A2); 

w=-V(iv)/(2*pi): 


for  it=l:tseg+l: 

thetal(it)=(pi/2)*((it-l)/tseg); 

for  izl=  l:zseg*l; 
zl(izl)  =  -TP*((izl-l)/zseg); 


for  ix  1 =1  :xseg+l : 
xl(ixl)  =-LP*((ixl-l)/xseg); 


saral(ixl)-w*((cp6*xl(ixl)A6)+(cp5*xl(ixl)A5)+(cp4*xl(ixl)A4)+(cp3*xl(ixl)A3)+(cD 

2*xl(ixl)A2)+(cpl*xl(ixl))-rcpO); 

katepl(ixl)=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it)))A2); 
kateql(ixl)=sin(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it)))A2);, 
greg  1  (ix  1  )=exp(k*sub6*(sec(thetal  (it)))A2); 

p  1  x(ix  1  )=sara1  (ix  1  )*katep  1  (ix  1  )*greg  1  (ix  1 ); 
q  1  x(ix  1  )=sara  1  (ix  1  )*kateq  1  (ix  1  )*gregl  (ix  1 ); 

end  %this  end  closes  the  ****  for  ixl  =  l:xseg+l  loop  **** 

p  1  z(iz  1  )=sum(p  1  x)/xseg; 
q  1  z(iz  1  )=sum(q  1  x)/xseg; 

end  %this  end  closes  the  ****  for  izl  =  l  :zseg+l  loop  **** 

plt(it)=sum(plz)/zseg; 

qlt(it)=sum(qlz)/zseg; 

%  *******This  starts  the  next  body  calculations  inside  the  theta  loop***** 


for  izlA=  l:zseg+l: 
zlA(izlA)  =  -TP*((izl  A-l  )/zseg); 

for  ixl  A=1  :xseg+l ; 
xlA(ixlA)  =-LP*((ixlA-l)/xseg); 

saral  A(ixlA)=w*((cp6*xl  A(ixl  A)A6)+(cp5*xlA(ixlA)A5)+(cp4*xl  A(ixlA)A4)+(cp3* 

x  1  A(ix  1  A)A3)+(cp2*x  1  A(ix  1  A)A2)+(cp  1  *x  1  A(ix  1  A))+cpO); 

kateplA(ixlA)=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it)))A2); 

kateq  1  A(ix  1  A)=sin(k*(sb6*cos(theta  1  (it))+os6*sin(theta  1  (it)))*(sec(theta  1  (it)))A2); 
gregl  A(ix  1  A)=exp(k*sub6*(sec(thetal  (it)))A2); 

p  1  Ax(ix  1  A)=saral  A(ix  1  A)*katep  1  A(ix  1  A)*greg  1  A(ix  1  A); 
q  1  Ax(ix  1  A)=saralA(ix  1  A)*kateq  1  A(ix  1  A)*greg  1  A(ix  1  A); 

end  %this  end  closes  the  ****  for  ixl  A=1  ixseg-^l  loop  **** 
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p  1  Az(iz  1  A)=sum(p  1  Ax)/xseg; 
ql  Az(izl  A)=sum(ql  Ax)/xseg; 


end  %this  end  closes  the  ****  for  izlA=l:zseg+l  loop  **** 

p  1  At(it)=sum(p  1  Az)/zseg; 
q  1  At(it)=sum(q  1  Az)/zseg; 

res66(it)=16*pi*rho*kA2*((plt(it)*plAt(it)Hqlt(it)*qlAt(it)))*(sec(thetal(it))A3); 

end  %this  end  closes  the  ****  for  it=  l:tseg+l  loop  **** 

R\v66(iv)=  sum(res66)/tseg 
end  °o  this  closes  the  ****  for  iv  =  ...  loop  ***** 
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%*********************ONESIX.m************************** 


Rw  1 6=zeros(hs+ 1  -Is,  1 ); 
res  1 6=zeros(tseg+ 1,1); 
theta  1  =zeros(tseg+ 1,1); 
zl=zeros(zseg+l,l); 
xl=zeros(xseg+l,l); 
p  1  x=zeros(xseg+ 1,1); 
qlx=zeros(xseg+l,l); 
p  1  z=zeros(zseg+ 1,1); 
q  1  z=zeros(zseg+ 1,1); 
plt=zeros(tseg+l,l); 
q  1  t=zeros(tseg+ 1,1); 
q  1  At=zeros(tseg+ 1,1); 
p  1  At=zeros(tseg+ 1,1); 
p  1  Az=zeros(zseg+ 1,1); 
q  1  Az=zeros(zseg+ 1,1); 
zl  A=zeros(zseg+ 1,1); 
plAx=zeros(xseg+l,l); 
q  1  Ax=zeros(xseg+ 1,1); 
x  1  A=zeros(  xseg+ 1,1); 
sara  1  =zeros(xseg+ 1,1); 
katep  1  =zeros(xseg+ 1,1); 
kateq  1  =zeros(xseg+ 1,1); 
greg  1  =zeros(xseg+ 1,1); 
sara  1  A=zeros(xseg+ 1,1); 
katep  1  A=zeros(xseg+ 1,1); 
kateq  1  A=zeros(xseg+ 1,1); 
greg  1  A=zeros(xseg+ 1,1); 
V=zeros(hs-t-l-ls,l); 


for  iv=l  :hs; 
V(iv)=ls+(iv*lani); 

k=g/(V(iv)A2); 

w=-V(iv)/(2*pi); 


for  it=l:tseg+l; 
thetal  (it)=(pi/2)*((it- 1  )/tseg); 


for  izl=  l:zseg+l; 
zl(izl)  =  -TP*((izl-l)/zseg); 
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for  ixl=l  :xseg+l; 
xl(ixl)  =-LP*((ixl-l)/xseg); 


saral(ixl)=w*((cs6*xl(ixl)A6)+(cs5*xl(ixl)A5)+(cs4*xl(ixl)A4)+(cs3*xl(ixl)A3)+(cs2 

*x  1  (ix  1  )A2)+(csl  *x  1  (ix  1  ))+csO); 

katepl(ixl)=cos(k*(sbl*cos(thetal(it))+osl*sin(thetal(it)))*(sec(thetal(it)))A2); 
kateq  1  (ix  1  )=sin(k*(sb  1  *cos(thetal  (it))+osl  *  sin(thetal  (it)))*(sec(theta  1  (it)))A2); 
gregl(ixl)=exp(k*subl*(sec(thetal(it)))A2); 

plx(ixl)=saral(ixl)*katepl(ixl)*gregl(ixl); 
ql  x(ix  1  )=saral  (ix  1  )*kateql  (ix  1  )*greg  1  (ix  1 ); 

end  %this  end  closes  the  ****  for  ixl=l  :xseg+l  loop  **** 

p  1  z(iz  1  )=sum  ( p  1  x)/'xseg; 
qlz(izl)=sum(qlx)/xseg; 

end  %this  end  closes  the  ****  for  izl  =  l:zseg+l  loop  **** 

p  1  t(it)=sum(p  1  z)/zseg; 
q  1  t(it)=sum(q  1  z)/zseg; 

%  **»****This  starts  the  next  body  calculations  inside  the  theta  loop***** 


for  izl  A=  1  :zseg+l ; 
zlA(izlA)  =  -TP*((izlA-l)/zseg); 

for  ixl  A=1  :xseg+l ; 
xlA(ixlA)  =-LP*((ixl  A-l  )/xseg); 

saralA(ixlA)=w*((cp6*xlA(ixlA)A6)+(cp5*xlA(ixlA)A5)-i-(cp4*xlA(ixlA)A4)+(cp3* 
x  1  A(ix  1  A)A3)+(cp2*x  1  A(ix  1  A)A2)+(cp  1  *x  1  A(ix  1  A))+cpO); 
kateplA(ixlA)=cos(k*(sb6*cos(thetal(it))+os6*sin(thetal(it)))*(sec(thetal(it)))A2); 
kateq  1  A(ix  1  A)=sin(k*(sb6*cos(  theta  1  (it))+os6*sin(thetal  (it)))*(sec(  theta  1  (it)))A2); 
greg  1  A(ix  1  A)=exp(k*  sub6*(sec(theta  1  (it)))A2); 

p  1  Ax(ix  1  A)=saral  A(ixl  A)*katep  1  A(ix  1  A)*greg  1  A(ix  1  A); 
q  1  Ax(ix  1  A)=saral  A(ix  1  A)*kateq  1  A(ix  1  A)*greg  1  A(ix  1  A); 

end  %this  end  closes  the  ****  for  ixl A=l:xseg+1  loop  **** 
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p  1  Az(izl  A)=sum(p  1  Ax)/xseg; 
q  1  Az(iz  1  A)=sum(q  1  Ax)/xseg; 

end  %this  end  closes  the  ****  for  izlA=l:zseg+l  loop  **** 

p  1  At(it)=sum(p  1  Az)/zseg; 
q  1  At(it)=sum(q  1  Az)/zseg; 

resl6(it)=32*pi*rho*kA2*((plt(it)*plAt(it))+(qlt(it)*qlAt(it)))*(sec(thetal(it))A3); 

end  %this  end  closes  the  ****  for  it=  l:tseg+l  loop  **** 

Rwl6(iv)=  sum(resl6)/tseg 
end  %  this  closes  the  ****  for  iv  =  ...  loop  ***** 
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